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On the fraction of star formation occurring in bound stellar clusters 
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ABSTRACT 

We present a theoretical framework in which bound stellar clusters arise naturally at the high- 
density end of the hierarchy of the interstellar medium (ISM). Due to short free-fall times, 
these high-density regions achieve high local star formation efficiencies, enabling them to 
form bound clusters. Star-forming regions of lower density remain substructured and gas- 
rich, ending up unbound when the residual gas is expelled. Additionally, the tidal perturbation 
of star-forming regions by nearby, dense giant molecular clouds imposes a minimum density 
contrast required for the collapse to a bound cluster. The fraction of all star formation that oc- 
curs in bound stellar clusters (the cluster formation efficiency or CFE) follows by integration 
of these local clustering and survival properties over the full density spectrum of the ISM, and 
hence is set by galaxy-scale physics. We derive the CFE as a function of observable galaxy 
properties, and find that it increases with the gas surface density, from T ~ 1% in low-density 
galaxies to a peak value of T ~ 70% at densities of E g ~ 10 3 M Q pc~ 2 . This explains 
the observation that the CFE increases with the star formation rate density in nearby dwarf, 
spiral, and starburst galaxies. Indeed, comparing our model results with observed galaxies 
yields excellent agreement. The model is applied further by calculating the spatial variation 
of the CFE within single galaxies. We also consider the variation of the CFE with cosmic time 
and show that it increases with redshift, peaking in high-redshift, gas-rich disc galaxies. It is 
estimated that up to 30-35% of all stars in the Universe once formed in bound stellar clus- 
ters. We discuss how our theory can be verified with Gaia and ALMA, and provide possible 
implementations for theoretical work and for simulations of galaxy formation and evolution. 

Key words: galaxies: ISM — galaxies: starburst — galaxies: star clusters — galaxies: stellar 
content — stars: formation — stellar dynamics 



1 INTRODUCTION 



The formation of stars from a turbulent interstellar medium 



ing giant molecular clouds (e.g. 


LarsorJ 1981: Elme 


jreen & Scald 


2004: Mac Low & Klessedl2004 


; McKee & Ostriker 2007). While 



this by necessity implies that stars generally do not form 
alone, it has been known for a long time that the star for- 
mation process does not exclusi v ely produce stars in bound 



stellar clusters (Elmegreerj 1 19831 : lLadd [l987). Instead, some 
fraction of stars is born in unbound stellar associations 
over a bro a d ran g e of physical scales 
Clarke etai] l200d: iMegeath et ail [20041: iPortegies Zwart et 
20ld : lGieles & Portegies Zwartll201ll ; lBastian et alj|2012h . and in 
rare cases individual star s may even form in relative isolation 



9641; 
iUd] 



IWhite & Reesl 1 19781 ; ISearle & Zinnl 1 19781 : Iwhite & Frenkl Il99lh 
several decades ago, the concept of star formation occurring in 
quantized, gravitationally bound systems has remained remarkably 
popular. 

It was pointed out by lLada & Ladal d2003h that in the Milky 
Way, embedded stellar groups constitute the vast majority of star 
formation and that these groups are about 10-20 times more numer- 
ous than gas-rid star clusters, suggesting that only a small fraction 
is capable of surviving the gas-embedded phase. With their review, 
the concept of infant mortality became firmly established - in this 
scenario, most (if not all) stars form in clusters, but the popula- 
tion of embedded clusters is t hen decimated due to the expulsion of 
residual gas by feedback (e.g. Tutukov 1978; Hills 1980; Lada et al. 
19841: lAdamsl I2OO0I; iGever & Burkertl 1200 ll; iBoilv & Kroun 



jParker & Goodwinl 120071 : iKrumholz etal] 120091 : iBressert et all 
l2012ah . Surprisingly though, star clusters are often still con- 
sidered to be a fundamental unit of star formation (e.g. 



Pflamm-Altenburg et alj 120071 : IPfalznej 12009] ; lAssmann etal] 
201 ll) . Whereas galaxy formation studies augmented t he scenario 
of mo nolithic collapse (see e.g. the seminal paper by Egg en et ail 
1 19621) with the current picture of hierarchical galaxy formation (e.g. 



2003allbl feastian & Goodwid 120061: iGoodwin & Bastiad l200d : 
Baumgardt & Kroupdl2007l : |Parmentier et all2008fc 

However, a rece nt study by Bress ert et al ] boid. also see 
IParker & Meveill2012[) of the spatial distribution of young stellar 
objects (YSOs) in the solar neighbourhood has shown that star for- 
mation follows a continuous spectrum of number densities, indi- 
cating that there is no separate or critical density scale for star 
cluster formation imprinted in the star formation process. Their 
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results imply that surface den sity thresholds used for the identifi- 
cation of stellar clusters (e.g . lLada & Ladall2003l ; |j0rgensen et al.l 
120081 : [Gutermut h~etai]|2009h are arbitrary, and do not correspond 
to a physical scale. It is thus impossible to conclude that the ma- 
jority of stars form in clusters by dividing the density spectrum of 
star formation at a certain density. In other work, it has been shown 
that Y SOs are strongly c orrelated with the hierarchically structured 
ISM jTestietalJ|200Ct I Allen et al J [20071 ; iGutermuth et"alll201ll ; 
iBressert et al.l l2012ch , and that the impact of gas expulsion on 
bound stellar structure in a dynamical envi ronment resulting from 
such a hierarchy is modest, if present at a ll jKruiissen et all 2012a; 
iGirichidis et al.l2012l ; ICottaar et alj2012h . Rather than the scenario 
of star cluster formation in which infant mortality (i.e. gas expul- 
sion) plays an important role, these results favour a framework in 
which the observed scarcity of bound, gas-rid clusters is a sim- 
ple result of the star formation process: only some fraction of star 
formation reaches the densities required to attain high star forma- 
tion efficiencies and result in bound stellar clusters, while the rest 
forms in a more dispersed fash ion throughout the natal cloud (e.g. 
Elmegreen & Elmegreen 2001). 

The fraction of star formation occurring in bound stellar clus- 
ters is ofte n quantified a s the cluster formation efficiency (CFE or 
r. see e.g.lBastiadl2008l:lGoddard et al]|201Cl : lAdamo et alj|20 111 ; 
ISilva- Villa & Larsenll201 1). The CFE is a crucial quantity in many 
respects. Not only can it ena ble a better understanding of the star 
formation process itself (e.g. lElmegreenll2002h . but it is also a key 
ingredient for work that aims to trace the (star) format ion histories 
of galaxies using their star cluster populations (e.g. Larsen et al 



200l |: |Bastian et al]2003:ISmith et al.l2007l : lKonstantopoulos et al 
2009; Fedotov et al. 201 1), or studies of star cluster disruption (e.g. 
Gieles et al.l 2005 ). On galactic scales, the CFE may be used to infer 



whether the most massive stars inject their feedback energy within 
the same bu bbles and hence how efficiently feedback coup les to 
the ISM (cf. IStrickland & Stevendfl999l ; lKrause et alj|2012h . Ad- 
ditionally, the CFE is an important parameter in numerical simula- 
tions on galactic or cosmological scales th at aim to model the as- 
sembly history of star cluster systems (e.g. lPrieto &Gnedin| [2008; 
iKruiissen et al.|[20ll I. l2012bl) . Thusfar, such simulations have had 
to assume that the CFE is constant throughout cosmic time. 

Observational studies of extragalactic cluster populations at 
ages older than the embedded phase suggest that the CFE increases 
with the star formation rate surface density of the host galaxy 
dLarsen & Richtledl200d:lGoddard et alj|20 id ; lAdamo et al.ll20Tll ; 
ISilva- Villa & Larsenl 201 ll) . but at present there is no theoretical 
understanding of the variation of the CFE with the galactic envi- 
ronment. It is notoriously hard to model the formation of stellar 
clusters in numerical simulations due to the large dynamic range 
that needs to be covered to resolve the necessary physics. Cur- 
rent state-of-the art si mulations can model systems up to masses 



of a few 10 3 M a (e.g . iTillev & Pudritzll 20071 ; iBonnell et al]|2008l ; 



iKrumholz et al.|[2012r) . but this is still insufficient by several orders 
of magnitude to address the formation of cluster populations on the 
desired galactic scales. 

The aim of this paper is to derive and apply a theory of clus- 
ter formation that is based on simple anal ytical consideratio ns. In 
part, it will follow the recent approach by Elmegreen (2008), who 
showed that the CFE can be related to the density spectrum of the 
ISM, but did so by defining a critical density (or pressure) for clus- 
ter formation. While not suitable to do quantitative predictions due 
to an arbitrarily defined density threshold, the work is a big step 
forward from the idealized picture in which centrally concentrated 
gas is expelled from spherically symmetric, bound clusters. In the 



present paper, the avenue suggested by lEImegreenl {2008) is ex- 
plored further and expanded to provide a self-consistent framework 
for star cluster formation in galaxies. 

The paper is organized as follows. In ij2] the required physi- 
cal mechanisms are discussed that should be included in a theory 
of cluster formation. The model is derived and presented in ^2.2\ - 
12.81 In Sj3]the parameter space is explored and it is addressed how 
the CFE varies among different galaxy types. The calculated CFEs 
are compared to observed cluster populations in Sj4] giving excel- 
lent agreement. In Sj5]the spatial variation of the CFE within single 
galaxies is addressed. In [J6]the variation of the CFE with cosmic 
time is discussed, accounting for the evolution of galaxy proper- 
ties in a cosmological context. A discussion of possible caveats is 
presented in SjT] together with an outlook to the observational ver- 
ification of the model, and potential applications in theoretical and 
numerical work. The conclusions of this paper are given in ij8] We 
provide Fortran and IDL routines for calculating the CFE with the 
model of this paper at http://www.mpa-garching.mpg.de/cfe (see 
Appendix [A}. 



2 FORMATION OF BOUND STELLAR STRUCTURE IN 
STAR-FORMING REGIONS 

This section discusses the physical mechanisms that should be in- 
cluded in a theory of the CFE, and presents the derivation of a 
model in which these mechanisms are accounted for. In the remain- 
der of the paper, the model is applied, verified and discussed. 



2.1 Relevant physical mechanisms 

The fraction of stellar structure that eventually ends up being bound 
after star formation has ceased is the result of several physical 
mechanisms that act on a range of spatial scales. These mecha- 
nisms postulate a set of requirements for a theory designed to de- 
scribe and explain the CFE. Observations, simulations and the- 
ory show that star formatio n is likely a universal process (e.g . 
l Elmegreen & Efremovl 1 19971: iKennicuttl 1 19981: lEImegreenl |2002|; 
IKrumholz & McKeell2005l ; iMcKee & Ostrikedl2007h . in the sense 
that the physics of how a collapsing cloud converts its gas to stars 
seem to hold everywhere in the Universe. This implies that the 
blueprint for the initial stellar population is to a large degree al- 
ready present in th e characteristics of the ISM at the onset of star 
formation (see e.g. |Gutermuth et al ] |201ll ; lLongmore et al]|2012bl : 
IBressert et al.ll2012br) . It is therefore essential that a theory of the 
CFE is based on the structure and properties of the ISM that govern 
the process of star and cluster formation (see point (i) below). 

In a hierarchical ISM, the t raditional picture of mono- 
lithic star cluster formation (e.g. iBaumgardt & Kroupal 120071) 
breaks down, as is illustrated by star-forming regions such as 
M16, M17, Taurus, Ophiuchus or even in the region around 
the Orion Nebula Cluster, which is the local archetype of clus- 
tered star formation. These regions may contain local concen- 
trations of stars (which are generally unbound after subtract- 
ing the gravitational potential of the gas, see e .g. Ide Zeeuw et al.l 
Il999l : IScallv etal.ll2005l ; iHuff & Stahlen 120061) . but on a global 
scale they are co nstituted by a stellar populat i on that is hierar- 
chical in nature dEfremov~ & Elmegreen] 1 19981 ; iTesti et al] l200d : 
IClarke et alj|200d :l Brcsscrt ct al . 2012c). Per consequence, obser- 
vational and theoretical studies do not find any evidence for dis- 



crete modes of star formati on in the density spectrum of stel- 
lar structure in such regions ( IBressert et al.ll20l"ol ; IGutermuth et al.l 
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1201 it iGieles & Portegies Zwarlll201 it iKruiissen et alj|2012ah . In- 
stead, dispersed and clustered stellar structure represent opposite 
extremes of a continuous density spectrum of star formation. Al- 
though it is tempting to define a critical density that separates both 
extremes, it should be noted that this is not possible, especially 
for a hierarchically structured distribution - at i ntermediate den- 
sities, only a certain part of the structure is bound dGoodwinlll997t 
Ide Zeeuw etai1ll999t IKruiissen et alj|2012al) and no hard separa- 
tion exists. This calls for a continuous description of cluster for- 
mation, in which the (partial) boundedness of stellar structure is 
evaluated locally (see point (ii) below). This will enable the formu- 
lation of the CFE across the stellar density spectrum, which in turn 
allows the CFE to be connected to the density spectrum of the ISM 
if one adopts a local prescription for star formation (e.g. Elmegreer] 



l2002t iKrumholz & McKeelf2005t IPadoan & Nordlun J201 ll) ~ 

It is irrelevant to the problem at hand whether unbound stars 
result from gas expulsion or from a possibly intrinsically unbound 
state of the (stellar) substructure. The hierarchical nature of star- 
forming regions implies that there can exist spatially concentrated 
group s of stars that each ind i vidually may or may n ot be unbound 
(see IKruiissen et al.l l2012al : iGirichidis et al.1 l2012h . while other, 
more distributed sequences of stellar groups could have ended up 
merging and becoming a star cluster had the gas between them not 
been expelled (de pending on the virial state of the cloud, see e.g. 
Dobbs et al. 2011b). This does not mean that gas expulsion is not 
important, but it does indicat e that the classical picture of infant 
mortality dLada & Ladall 20031) due to gas expulsion from centrally 
concentrated embedded clusters does not hold in a picture of hier- 
archical star formation. In such a dynamic environment, a singular 
observational definition of a cluster in terms of current or future po- 
tential groups is not possible, and therefore the CFE requires spec- 
ifying an age t at which it is evaluated (see point (iii) below). 

While the relative amounts of clustered and distributed star 
formation may already be largely set by the local properties of the 
ISM, it has also been shown that environmental effects can play an 
important role during the aftermath of star and cluster formation. 
This is to be expected - given a universal star formation process, 
the characteristics of the ISM would only contain the complete 
blueprint for the state of the new-born stellar structure if star forma- 
tion were an instantaneous mechanism. But because the star forma- 
tion process takes time, environmental effects come into play as the 
conversion of gas into stars proceeds. The truncation of star forma- 
tion by feedback and the potential dispersal of the stellar structure 
may appear to be mostly an internal process, but theoretical and 
numerical work have shown that the effectiveness of gas expulsion 
is largely determined by how the feedback couples to the ISM (e.g . 



is largely determined by now the feedback couples to the IsM (e.g. 
Elmegreen & Efremovi ll997l : IPelupessv & P ortegies ZwarJ 12012 ). 



This implies that the time-scale on which the inflow of gas onto a 
star forming region can be halted depends on a combination of the 
porosity and pressure of the ambient ISM (see point (iv) below). 

Another environmental eff ect is entirely external in nature. 
In lKruiissen et al.l j201ll , l2012d) we identified the cruel cradle ef- 
fect, which refers to the tidal disruption of star-forming regions or 
young stellar clusters by encounters with giant molecular clouds 
(GMCs) or other substructure in the gas of th e dense, natal envi- 
ronment (also see Elmegreen & Hunter 2010). If the density con- 
trast between the star-forming region and surrounding GMCs is low 
enough, tidal shock£]are capable of destroying the new-born stellar 



1 The similarity in nomenclature with hydrodynamic shocks is unfortunate, 
especially in the context of this paper. However, tidal shocks are completely 



structure before the primordial gas has been cleared. The cumula- 
tive effect of this mechanism increases as star formation proceeds. 
The cruel cradle effect thus needs to be included in a theory of the 
CFE (see point (v) below). 

It is clear that the fraction of star formation producing bound 
stellar clusters is the result of a combination of factors, some of 
which are already imprinted in the ISM before the onset of star 
formation, while others are environmental effects that act on the 
stellar structure throughout. The requirements for a theory of the 
cluster formation efficiency are summarized as follows. 

(i) The theory must account for the hierarchical structure of the 
ISM, and do so in a continuous manner, i.e. without invoking any 
arbitrary thresholds for the formation of stars or stellar clusters. 

(ii) The dispersed and unbound fractions of star formation 
should follow from a local criterion that can be applied to the entire 
hierarchy of the ISM. 

(iii) The theory should enable predictions that are resolved in 
time, to account for the continuous nature of cluster formation. 

(iv) The truncation of star formation by feedback has to be ac- 
counted for, as well as its dependence on the ambient ISM (unless 
star formation is so efficient that the gas is depleted before then). 

(v) The disruption of substructure by tidal shocks acting within 
the star formation time-scale should also be included. 

A model that satisfies the above requirements enables the predic- 
tion of the CFE as a function of the galactic environment. This can 
be between individual galaxies as well as locally resolved, both in 
space and time. 



2.2 Outline of the model 

The cluster formation efficiency T can be formulated as the product 
of two fractions: 



100% 



/"bound fa 



(1) 



The first fraction, /bound, accounts for the naturally bound part of 
star formation and includes the effects of gas expulsion by feed- 
back, thus incorporating the second, third and fourth points of the 
set of requirements postulated in the previous section. The second 
fraction, / cce , covers the fifth of these points and indicates the frac- 
tion of potentially clustered star formation that is left after applying 
the tidal disruption due to the cruel cradle effect. 

Expressions for /bound and / ccc are derived in detail below, 
but the theory can be summarized as follows. We translate the den- 
sity spectrum of the hierarchical ISM into a local clustering of star 
formation. This can then be integrated to determine which part of 
star formation occurs in bound stellar clusters. The model consists 
of the following steps. 

(i) The starting point is the overdensity probability distribution 
function (PDF) of the ISM. This PDF reflects the distribution of 
density contrasts with respect to the mean density. (*j23J 

(ii) By assuming that star formation occurs in a gas disc that 



unrelated to fluid discontinuities, and refer to the transient injection of en- 
ergy into a gravitational system by the passage of a massive object such 
as for instance a GMC (see e g. ISpitzei|[l987r. iKundic & OstrikeJ 1 19951 ; 
iGnedin et al]| 19991; IGieles et alj|2006hl; IKruiissen et alj|201ll) . Throughout 
the paper, both types of shocks are distinguished using the adjectives 'tidal' 
and 'supersonic'. 
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obeys hydrostatic equilibrium, the mid-plane density can be de- 
rived, allowing the overdensity PDF to be written as a PDF of the 
absolute density. ( i]2.4t 

(iii) Given a certain (absolute) density, the local free-fall time 
is known. By adopting a specific star formation rate per free- 
fall timeQ and determining when the feedback pressure is strong 
enough to prevent the further inflow of cold gas and hence halt star 
formation, it is possible to determine the final star formation effi- 
ciency at each density. ( i]2.5l > 

(iv) For a local star formation efficiency, it is possible to formu- 
late a local fraction of stars that remains bound upon instantaneous 
gas expulsion. To establish this relation, we use a numerical simu- 
lation of star formation in a turbulent ISM. This then provides the 
naturally bound fraction of star formation at each density. ( i!2.6l > 

(v) It follows from the classical theory of tidal shocks in stellar 
clusters that there is a certain overdensity above which regions can 
survive the tidal perturbations by their environment. The theory is 
adjusted slightly to be applicable to gas-rich star-forming regions. 

(vi) The final expression for the CFE is obtained by integration 
of the naturally bound fraction of star formation over the density 
range of the PDF where it survives the cruel cradle effect, and di- 
viding it by the integral of the star formation efficiency over the 
entire density range of the PDF. ( i]2.8t 



The reader who is not so much interested in the details of the 
derivation below is referred to i j2.8l where the theory is summa- 
rized graphically in Figure Q] 



2.3 The density spectrum of the interstellar medium 

The hierarchical structure of the ISM is driven by a combina- 
tion of turbulence and local gravitational contraction. It has been 
known for well over a decade now that the PDF of the mass 
density in isothermal, supersonically turbulent clouds and larger 
structures of cold gas is well-described by a log-normal function 
( Vazquez-Semadeni 1994; Padoan et al. 1997; Padoan & Nordlund 
20021; IFederrath et all I2008L l20ld ; IPadoan & Nordlund 1201 ll ; 
Hill et alJl2oT3T When formulated as the PDF of the m'erdensity 
with respect to the mean density in the turbulent region x = 
Pg/ Pism, it can be written independently of the physical scale as: 



dp 
dx 



■ exp 



(In a; — Inx) 2 



2ai 



(2) 



where the logarithmic mean In x is related to the standard deviation 



i p 

lnx = - T , 

and 



(3) 



a 2 = ln(l + b 2 M 2 m ) =m(l + 3& 2 M 2 ), (4) 

with b « 0.5 (e.g. IPadoan & Nordlund! 12002). which is consis- 
tent with a mix of sol enoidal and compressive turbulence forcing 
( IFederrath etal.ll20ld) . and A4:iu and M the three-dimensional 
and linear Mach numbers, respectively. Equations (O-© hold at 
scales smaller than the scale at which the turbulence is driven, 
which in the case of disc-like ga s distributions has an upper l imi t 
equal to the disc scale height (e.g. Ma c Low & Klessenll2004l) . The 

2 This is the fraction of the gas that is converted into stars per free-fall time. 



smallest scale on which these equations hold is set by the ambipo- 
lar diffusion length, i.e. a few hundreths of a parsec, or the length- 
scale of protostellar cores {Ossenkopf & Mac Low 2002). It was re- 
cently shown that the form of the overdensity PDF is very similar in 
the presence of strong magnetic fields, requiring only a correction 
factor in the second term of equation ([4} of /3n/(/3o + 1), where 
Po = -Pth/f mag is the ratio of the thermal pressure to the mag- 
netic pressure (Padoan & Nordlund 2011; Molina et al. 201 3). For 
simplicity, magnetic fields are omitted in the present model for the 
CFE, although it should be noted that the adopted value of b = 0.5 
is roughly consistent with recent results of magnetohydrodynam- 
ical simulations with weak magnetic field strengths (in which the 
Al fvenic velocity is co mparable to or smaller than the sound speed, 
see lMolina et alj2012l) . 



2.4 Gas discs in hydrostatic equilibrium 

The above set of equations shows that the density PDF solely de- 
pends on the Mach number of the g as M and the mid-plane den sity 
of the galaxy disc pism. Following Krumholz & McKee (2005) we 
assume that the gas disc of the galaxy is in hydrostatic equilib- 
rium[f] which allows both quantities to be expressed in terms of 
the me an surface density E g , angular velocity f2, and the iToomrel 
dl964l) stability parameter Q. This implies that the overdensity PDF 
becomes a PDF of the absolute density, which is set by those three 
variables. In this formulation, the mid-plane density is given by 



PISM 



tvGQ 2 



= 2.8 x 10~ 21 Q~ 2 fio g cm" 3 , 



(5) 



where <f>p — 3 is a constant to account for the gravity of the stars 
and fin = fl/Myr -1 is a convenient scaling of the angular veloc- 
ity. By deriving an expression for the velocity dispersion in GMCs 
and comparing it to the typical sound speed in the cold ISM, the 
Mach number of the star-forming regions in a galaxy disc is ap- 
proximated as 

M = 2.82<j>^ s QQo ^,2, (6) 

with E gi 2 = Eg/10 2 Mq pc -2 . In this equation, (f>-p is the ratio of 
the mean GMC pressure to the mid-plane disc pressure and can be 
expressed as 

<j>p = 10 - 8/gmc « 10 - 8 (l + 0.025E" 2 ,) ~' , (7) 

where /gmc is the mass fraction of the ISM in GMCs. The lToornrel 
(1964)Q parameter is defined as 



Q 



7rGE E 



7tGE s 



(8) 



in which k is the epicyclic frequency, <7 g the one-dimensional ve- 
locity dispersion of the gas, and fi the angular velocity within the 
galaxy. The second equality assumes a flat galaxy rotation curve. 
Gas discs with Q < 1 are considered unstable to gravitational col- 
lapse, whereas Q > 1 indicates stability by kinetic support. As a 
fiducial value Q = 1.5 is assumed, but it is known to vary between 
0.5 and 6 (e.g. lKennicuttlll989l : lMartin & Kennicuttll2001l) . Finally, 
the angular velocity can be related to the surface density as 



3 Perhaps surprisingly, this assumption is not too inaccurate for (e.g. 
merger-induced) starburst galaxies, in which the rapid cooling of star- 
forming gas implies that a substanti al fraction of the st ar-forming regions 
follows a disc-like morphology (see Hopkins et al. 2009 and i]7.U . 
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fi = 0.058E 



g,2 ! 



(9) 



where the Mach number is estimated from E E 



and Q, as 



with an intrinsic scatter of about 0.5 dex in the range E gi 2 = 10 -2 - 
10 3 . Note that this is a fit to observed, nearby galaxies in their en- 
tirety, implying that Q. should still be treated as an independent vari- 
able when the spatially resolved CFE within a model galaxy is cal- 
culated, or when high-redshift galaxies are considered. Moreover, 
if the rotation curve is not flat, Q. should be replaced by 



V2 #gc 



1 + 



din v c 
d In R e , 



(10) 



in all model equations. In this definition, v c indicates the circular 
velocity, and R gc represents the galactocentric radius. The quan- 
tity k/V2 increases from Q. for flat rotation curves to y/2£l for 
solid body rotation. Neglecting this variation thus introduces an er- 
ror of up to a factor \/2 in the angular velocity, which translates to 
a smaller uncertainty on the CFE (see [j3]for the variation of F with 
Q). 

For details regarding the d erivation of equations (6) , (TJ, 
and (9]l the reader is referred to iKrumholz & McKeel (20051) . In 
combination with equations (2)-(4), this set of equations defines 
the absolute density PDF as a function of the gas surface density 
Eg, the Toomre Q parameter, and the angular velocity fi. 



before, and ow is the virial parameter of a typical GMC (see 
lBertoldi&McKeelll992h . The virial parameter is given by 



5aLR 
GM ' 



(13) 



in which M is the cloud mass, R is the cloud radius, and <Ji nt is the 
internal velocity dispersion, which is assumed to equal the turbu- 
lent velocity dispersion because cold GMCs are highly supersonic. 
GMCs with Q v i r ~ 1 are in virial equilibrium, while those with 
Qvir > 1 are not self-gravitatinjjf] and those for which a < 1 are 
contracting or supported by magn etic pressure. A fiducial value of 
Qvir = 1-3 (McKee & Tan 2003) is adopted here, but it should be 
noted that values in the range ow = 0.2-10 are found in observa 



tions (Solomon et al. 1987; Heyer et al 



2009) and numerical work 



dTasker & Tanll2009l : lDobbs et al.ll2011bl ; lHopkins et alj|2012l) . For 
typical values of a v ir = 1.3 and M = 100 (this is appropriate 
for E g , 2 ~ 1), equation (12) gives sSFR| M05 = 0.014, very 
similar to sSFRf 02 = 0.012. In our fiducial model, the empirical 
star formation law of equation (TTJ will be used; the influence of 
the adopted star formation law (and, per consequence, of the GMC 
virial parameter) on the CFE is discussed in more detail in Sj3] 



2.5 The local star formation efficiency 

Given the density PDF of the turbulent gas, we can compute the 
fraction of the mass that ends up in stars (the star formation effi- 
ciency, SFE or e) as a function of density and thereby establish how 
much each density contributes to the total amount of star formation 
in a galactic region. 



2.5.1 The specific star formation rate per free-fall time 

Assuming that st ar formation proceeds on a free-fall time ts (e.g. 
Elmegreen 2000), the star formation efficiency can be expressed in 
terms of the specific star formation rate per free-fall time sSFRg, 
i.e. the mass fraction of a GMC that is converted into stars each 
free-fall time. For this, it is needed to assume a star formation law, 
whi ch has been the topic of extensive debate in the recent literature 
(e.g.lElmegreej2002l;lKrurnholz & McKej2005l;lKrumholz & far] 
l2007l ; lElmegreenll2007l ; |Padoan & Nordlundll201 ll) . A star forma- 
tion law specifies sSFRgp, which interestingly exhibits little varia- 
tion over a large dynamic range. Obs ervational results on galactic 
scales suggest that sSFRg ~ 0.01 (Kennicutt 19981; lElmegreenl 
l2002h . whereas in the Cores to Disks (c2d) Spitzer Legacy project, 
a value sSFRg ~ 0.04 is found within nearby star-forming regions 
jEvans et al]|2009h . Considering the difference in physical scales, 
this relative consistency is remarkable. In this paper, the CFE is de- 
rived using two different star formation laws. The first is based on 
the empir ically motivated a ssumption that sSFRff is approximately 
constant (Elmegreen 2002): 



sSFR| 02 



0.012. 



(11) 



The second st ar formation law that i s used in this work is the pre- 
scription from IKrumho lz & McKee (2005). They estimate sSFRg 
from the part of the overdensity range that can be interpreted as 
protostellar cores: 



sSFR KM05 



0.13 



, oi -ln0.68a 2 ir X 4 
1 + erf 1 



2 3/2 



(12) 



2.5.2 The end point of the star formation process 

If the star formation process freely continues until it is halted by 
feedbackjf] one can formulate a feedback time-scale tft, to char- 
acterise the duration of star formation. The SFE then becomes a 
simple function of sSFRg , iff, and tib'. 



sSFR ff 
e = — £fb- 

iff 



(14) 



In this equation, the free-fall time tg is defined for each density p g 
as 



iff 



3tt 



32G/9 g 



(15) 



The feedback time, i.e. the time it takes to halt star formation, 
is taken to depend on when pressure equilibrium is attained be- 
tween feedback and the surrounding ISM. There is some debate in 
the literature as to which feedback mechanisms halts star formation 
on the spatial scale of entire GMCs. There are indications that su- 
pernova feedback is the most imp ortant agent for truncating the star 
forma t ion process on thes e scales ( Larson 1974; McKee & O strikeJ 
19771: iKorpi et all 1 19991 : Ijoung & Mac Low! 120061 : iDobbs et al.1 
201 lal) . However, it has recently also been shown that supernova 



4 Note that this does not imply zero star formation, since part of the GMC 
will be gravitationally unstable. 

5 The Galactic Center provides an exception to such a scenario due to its 
high density and angular velocity: the Arches cluster is thought to have 
stopped forming stars due to a cloud-cloud collision, which would explain 
why it is no longer associated with a GMC at its present age of ~ 2-4 Myr 
(Wang et alj |2006). For regions where the time-scale for cloud-cloud col- 
lisions is much shorter than the feedback time-scale, i.e. i co n <C tf^, the 
feedback time-scale in equation 1141 can be replaced by the cloud -cloud 
collision time-scale t coii = 2QEQMC / '4> pf2£g (see e.g. lSillJl997l) . with 
Eqmc the typical surface density of GMCs. However, on a global scale the 
formation of stellar clusters in galaxy discs is generally unaffected, because 
the cloud-cloud collision time-scale increases approximately linearly with 
galactocentric radius, and greatly exceeds tft, outside the central 100 pc of 
the Milky Way. 
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feedback could be ineffective in high-density galaxies due to effi- 
cient coolin g and instead radi ative feedback would be the dominant 
agent (e.g. [Thompson et al.l I200I iFall et alj I20T0I : iMurrav et all 
l20ldlDaleetalJl2012h . Throughout this paper, we will assume 
that supernovae truncate the star formation process on the largest 
spatial scales, but in Appendix |Cl we also consider the (additional) 
effect of radiative feedback. The relevance of supernovae is sup- 
ported to some degree by numerical simulations of embedded clus- 
ters - those clusters that span the largest spatial scales and have 
the weakest coupling between feedback and the ISM, reflecting a 
high degree of inhomogeneity, have dynamical histories that are 
most strong ly impacted by supernovae rather th an winds from mas- 
sive stars jPelupessv & Portegies Zwartll2012h . Crucially though, 
the description of feedback in our model serves the purpose of trun- 
cating star formation on a time-scale that is broadly consistent with 
observations, rather than reflecting the details of the feedback pro- 
cess. Additionally, Appendix [C] shows that varying the feedback 
prescription has only a modest influence on the predicted CFE. 

We write the feedback time as the sum of the time since the 
onset of star formation until the first supernova t sn and the sub- 
sequent time until pressure equilibrium between feedback and the 
surrounding ISM t eq : 



tfb — tsn te 



(16) 



At times t < tfb, star formation still is still ongoing because the 
ambient ISM pressure Pism is higher than the outward pressure 
provided by feedback Pfb, and gas is continuously fed into the star- 
forming region. The ambient pressure that needs to be overcome to 
halt star formation in GMCs is dominated by turbulent motion and 
is therefore given by 



2 7T G Q E g piSM 
P[SM = PlSMOg = ^2 : 



(17) 



where equation (f8]( was used in the second equality. The condition 
for pressure equilibrium depends on how the feedback couples to 
the ISM - since the structure of the ISM is hierarchical, its por osity 
cause s the pressure support by feedback to be inefficient (e.g.|SilkJ 
1 19971) . Taking this into account, equilibrium can then be formulated 
as 

D Pfb / , , , 7r 2 G 2 Q 2 pisMSg 

Pfb = — = <PfbP S Ioq = (PfbCPglcq = ^2 ' 

where Pfb is the feedback energy providing the outward pres- 
sure, V is the volume of the region, p s is the stellar density, 
and </)fb is a constant that represents the rate at which feed- 
back injects energy into the ISM per unit stellar mass. It thus 
includes a certain efficiency factor smaller than unity and in- 
cludes the porosity of the ISM as well as additional energy loss 
due to cooling. While it is known that the feedback-ISM cou- 
pling is crucial in s haping the properties of galaxies in numeri- 
cal simulations (e.g. Efstathiou 2000; Springel & Hernquist 2003; 



2001 ; 
2008; 



Abadi et al. 2003; Roberts on et al.H2005l ; iDubois & Tevssier 2008 
Dalla V ecchia & Schave 2008), it s characteristics are still uncertain 
dSilkll 19971; iMac Low & Ferrarall 19991 ; iNavarro & Steinmet j|2000l ; 
Dib et a IT l2006l) . and therefore the value of 0fb has not been de- 
termined conclusively. In Appendix [B] <j>fb is chosen by consid- 
ering order-of-magnitude estimates from the literature, and it is 
shown that the influence of this choice on the CFE is minor. In 
the remainder of the paper, we adopt <^>fb ~ 0.16 cm 2 s -3 = 
i2 erg s _1 Mq -1 . From equation (18K the equilibrium 



3.2 x 10 



_ tt 2 G 2 Q 2 E 2 
* cq ~~ 2(f> fb exQ. 2 ' 



(19) 



where x = p g /pisM as in equation ©. 

Equations J14t . l |16t , and l |19t can be solved for the SFE e and 
the feedback time tfb. The solution for tfb reads 



tfb - 2 I 1 + Y </>fbsSFRfft 2 n fi 2 x 



and the SFE is given by 

sSFRff t s 



£fb = 



1 + + 



27r 2 G 2 tffQ 2 E 2 



b sSFRfft 2 n fi 2 :r / ' 



(20) 



(21) 



where the subscript 'fb' indicates that this is the SFE after star for- 
mation is halted by feedback. For a solar neighbourhood-like pa- 
rameter set, an overdensity of x = 10 1 , and t en = 3 Myr, the 
feedback time-scale is typically tfb ~ 5 Myr (with etb ~ 0.01), 
and it quickly approaches to, = t En for x 10 2 . 

In practice, star formation need not be halted by feedback as 
the region may either be observed when the star formation process 
is still ongoing, or the density may have been so high that all of 
the available gas was consumed before tfb. If star formation is still 
ongoing, the SFE at time t follows from a slight modification of 
equation dl4b as 

sSFRff 

fine = — 1, (22) 

£ff 

where the subscript 'inc' indicates that the conversion of gas to stars 
is still incomplete. Of course, ej Ilc depends on the moment of eval- 
uating the CFE, and a standard value of t — 10 Myr is assumed 
here. However, by varying t one can obtain the time-evolution of 
the CFE, which is shown in £13.31 If star formation proceeds at such 
a high density that it is optimally efficient before the onset of feed- 
back, the SFE equals the core star formation efficiency e corc : 

(23) 

where the subscript 'max' indicates that this is the maximum SFE. 
The value of e corc follows from the mass fraction of protostellar 
cores that ends up in the actua l stars after account i ng for the mass 
loss in outflows. According to lMatzner & McKed d200oT) . e corc = 
0.25-0.75, and a constant value of e corc = 0.5 is therefore adopted 
in this paper. Finally, this implies that the actual SFE is given by 



e = min {e max , eib , ftnc } • 



(24) 



The above treatment of the feedback-induced truncation of 
star formation assumes that the specific star formation rate per free- 
fall time sSFRff is constant until star formation is halted. This might 
not actually hold once the first supernovae have started to inject 
energy into the ISM, as the heating might eit her prevent or stim- 
ulate some regions from undergoing collap se dElmegreen & Lad3 
1 19771; bale et al.ll2005l ; IWunsch et alj 2008 ) . Due to the ambiguous 
effects of feedback, it is not clear a priori whether a gradual trun- 
cation of star formation would be appropriate. A simple correction 
is explored in Appendix [B] which shows that the change of the re- 
sulting CFE is very marginal either way. This justifies the use of a 
constant sSFRff . 

Another word of caution is in order regarding the possible 
variation of the onset and strength of feedback, covered in the con- 
stants t sn and (j>fb. These may depend weakly on metallicity (see 
the discussion in Appendix IB1. but they are affected much more by 
the stochasticity of star formation at the high-mass end of the mass 
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range. The sampling of the masses of the few most massive stars 
determines whether a continuous (and in that sense statistical) ap- 
proach as in equation J21t is suitable. It should therefore be kept in 
mind that the treatment of feedback is idealised in the sense that it 
may hold for an ensemble of star-forming GMCs, but not necessar- 
ily for individual examples. 



2.6 The naturally bound part of star formation 

It has long been known that the CFE depends on the SFE on the 
local scale of young stellar clusters. The less efficient star for- 
mation is locally, the smaller the part of the stell ar structure that 
remains bound after the gas has been expelled dTutukovl 1 19781; 
Hilldll98d:lAdamdl200d:lGever & BurkertlboOlklBoilv & Kroupal 
2003allbl: iGoodwin & Bastianl 120061: iBaumgardt & Kroupal 120071 ; 
Parmentier et al.l l 2008; Smi th et alj|201ll) . However, the details of 
the relation between the local CFE 7 and the SFE e have so far 
always been determined for a static background potentials and/or 
some degree of equilibrium between the gas and stars. It is an im- 
portant property of star-forming environments in nature that the 
gas and stars can evolve independently of eac h other. The decou- 
pled kinematics of both components (see e.g. lOffner et alj|2009l) 
and the accretion of the gas in the direct vicinity of the stars im- 
ply that stellar structure can evolve to a gas-poor state while re- 
maining embedded in an eva cuated cocoon of gas dKruiissen et al.l 
2012ahlGirichidis et alJl2012T) . As such, the assumption of dynam- 
ical equilibrium between gas and stars does not hold. For the pur- 
pose of this paper it is therefore not possible to follow the 'classical' 
equilibrium-type results to obtain 7(e). Instead, numerical simula- 
tions of turbulent fragmentation can be used to establish what the 
fraction of stars is that is gravitationally bound on some scale, for 

a given gas-to-stellar mass ratio or SFE. 

In Appendix IE] the analy s is of t he lBonnell et all d2008h sim - 
ulations by Maschber ger et alj d2010h andytruiissen e t al. 1 2012al) 
is used to address the local CFE as a function of the local SFE in 
a hierarchical star-forming region. Since the simulation does not 
include feedback, the effect of gas expulsion is accounted for by 
ignoring the gravitational potential of the gas when calculating the 
boundedness of the stellar structure. This reflects the most violent 
form of feedback, as the gas is removed instantaneously. Obviously 
this is an idea lized approach. A more realistic treatme nt of feed- 
back (see e.g. lWang et"ai]l2010l : lKrumholz et alj|2012h may work 
in two ways. Either it could accelerate the evolution of stellar sub- 
clusters to a gas-poor state, in which case the approach is still valid, 
or it could slow down star formation to such a degree that the stars 
and gas retain some degree of dynamical equilibrium. It is shown 
in Appendix ID1 that the extreme case of adopting the results from 
iV-body simulations that assume complete dynamical equilibrium 
between gas and stars does not substantially change the resulting 
CFE. Another implication of th e absence of f e edbac k is that star 
formation is quite rapid in the Bonnell et al. (2008) simulation, 
with sSFRft ~ 0.3. As a result, the collapse into stars proceeds so 
rapidly that turbulence may not have reached statistical equilibrium 
yet, which could drive the density PDF away from the characteris- 
tic log-normal shape if it is generated by the turbulence. However, 
it has been shown th at turbulence is no t required to produce a log- 
normal density PDF dTassis et alj201Ct) and indeed we find that the 



gas in the central 5 pc of the simulation is roughly consistent with 
the log-normal of equation $2^ for A4 — 3-10. 

The analysis in Appendix 151 shows that in the boundedness 
of stellar structure in the iBonnell et alj 120081) simulation can be 
described with 30% accuracy by the simple relation 7(e) = e. 
This does not account for protostellar outflows, i.e. the maximum 
SFE in the analysis is 100%, while in the model of this section, 
a maximum SFE equal to e corc = 0.5 is assumed. The ques- 
tion thus rises whether the local CFE could still be optimal at a 
maximum SFE that is smaller than unity due to outflows. Recent 
observational and numerical results on protostellar outflows show 
that they should be capable of driving the turbul ence on scales of 
~ 1 pc, but not on the scales of entire GMCs (e.g.|Arce et alj2O10t 
IWang et alj20l3 : lHansen et alj2012t iBuckle et al.l2012h . As such, 
their influence may n ot reach very far, contrary to earlier theoreti- 
cal expectations (e.g. Matzner & McKee 2000). The typical proto- 
st ellar outflow velocities of ~ 0.5 km s _1 in the numerical work 
of iNakamura & Lil d2007h are consistent with this picture, and for 
protoclusters with radii of 1 pc they imply crossing times of 2 Myr. 
Since t his is longer than th e typical free-fall time of such protoclus- 
ters (cf. lEvans et al. 2009), the stellar structure can respond adiabat- 
ically to any mass loss in protostellar outflows. The mere existence 
of bound, young stellar clusters without a massiv e halo of escaping 
stars (e.g. Roch atTet al.l2010llCottaar et al.l20 12) evidences that in 
at least some cases protostellar outflows alone are not sufficient to 
unbind stellar clusters. The local CFE should thus not be strongly 
affected by the outflows. This leads to the following adopted rela- 
tion for the local CFE: 

7 = e Acorc, (25) 

which increases linearly from to 1 as the SFE increases from to 

^corc • 

The naturally bound part of star formation /bound from equa- 
tion (0 can now be obtained by integration of the overdensity PDF 
of the ISM: 



l(x)e(x)x(dp/dx)di 
e(x)x(dp/dx)dx 



(26) 



where the numerator represents the part of star formation that re- 
sults in bound structure, and the denominator denotes the total 
amount of star formation. The quantities that set /bound are the 
gas surface density E g , the Toomre Q parameter, and the angular 
velocity 



2.7 The cruel cradle effect 

Star-forming regions are characterized by their enhanced densities 
relative to the mean ambient density in a galaxy disc. Over the 
course of a star formation time-scale, it is possible that these re- 
gions encounter other density peaks or GMCs with which they in- 
teract gravitationally. Such interactions result in tidal shocks, dur- 
ing which the perturbation injects energy into the star-forming re- 
gion. This may prevent collapse on a global scale, but because the 
susceptibility of an object to tidal shocks is set by its density, local 
condensations into protostellar cores are not inhibited. Tidal shocks 
thus decrease the CFE, but not the SFE. 



B Except for a normalization of the gas potential that decreases with time 
when the gas is expelled. 



7 As mentioned previously, if the galaxy rotation curve is not flat, f2 should 
be replaced by equation UOt in all equations. 
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2.7.1 The Spitzer theory of tidal shocks 

The disruption time-scale of a hierarchically structured star- 
forming region due to tidal shocks t CC c can be derived from 
the energy consid e rations that also hold for stellar clusters (see 
e.g. ISpitzerl Il987l; iKundic & Ostriken Il995l: iGnedin et alj Il999l : 
IGieles et al.ll2006bl : IPrieto & Gnedinll2008l; iKruiissen et aljhoilh . 
The only modifications lie in the adopted proportionality constants, 
which account for the physical and structural differences between 
star clusters and hierarchical star-forming regions. The time-scale 
is written as 



^cce — 



E (dEy 1 

1 V&t) 



(27) 



where / is a constant that accounts for the fraction of the energy 
gain that is used to unbind the star-forming region, E is the in- 
ternal energy per unit mass of the region, and the time derivative 
indicates the energy gain per unit mass due to tidal perturbations. 
It is assumed that the relevant regions experience more than one 
encounter before t = 10 Myr, which is verified in i]2.7.2l below. 

In order to calculate the internal energy of a star-forming re- 
gion, it is convenient to assume a certain radial density profile. 
While such spherical symmetry in principle violates the philoso- 
phy of this work to treat the ISM as a hierarchy, it is a reasonable 
approach since stars and stellar clusters are born in local density 
peaks. This assumption is particularly appropriate in the context 
of tidal perturbations: the details of the density profile of the per- 
turbed region transl ate into secon d-order changes in the proportion- 
ality constants (e.g. Spitzer 1987). For mathematical sim plicity, the 
regions are thus considered to follow a lPlummeil dl91 ll) potential: 



$ = 



GM 



i — 

Rl 



-1/2 



(28) 



where $ is the potential energy, M the mass of the region, Rq the 
characteristic (Plummer) radius, and R the radius. For the above 
potential, the internal energy per unit mass is given by 

3tt GM 



E 



64 Ro 



(29) 



The rate of energy change per unit mass dE/dt of the star- 
forming region can be derived acc ording to the formalis m for the 
perturbation of stellar clusters of ISpitzed dl958l Il987l) . but ap- 
pli ed to the gas-rich, s t ar-fo rming region as a whole. As laid out 
bv lBinnev & Tremaind dl987l) . dE/dt is given by the phase-space 
integration of the product of the encounter rate and th e energy in- 
jected b y each encounter. Generalizing their result as in Gie les et all 
(2006b) and assuming a Maxwellian velocity distribution, one ob- 
tains 



dE 4ir 3/2 gG 2 (f> sh 
dt 3cr g 

4V2^gG<l> sh {l 



3E g Q 



£GMCPlSM-R 2 <^>ad 
EGMCPlSM-R 2 0ad, 



(30) 



where g ~ 1.5 is a correction factor to account for the ex- 
tended nature of the perturbing clouds (s ee below), </> s h ~ 2.8 ac- 
counts for the higher-order energy gairFl dKundic & OstrikeJ 19951: 
Gnedin & Ostriken [l997l ; IPrieto & Gnedinl 120081 : IKruiissen et al.l 
201 ll) . Sgmc is the surface density of GMCs, and equation l[8j 



8 Calculated for iPlummed 1 1 9 1 ll) parameters using equation (23) of 
IKruiissen etai]d201ll) . 



was u sed in the second equality. The factor a d is the ISpitzerl 
dl987l) adiabatic correction, which accounts for the dampening 
of the energy injection by the adiabatic expansion of the re- 
gion. The expression for a d is derived below in H2J.2\ For the 
Milky Way and galaxies in the Local Group , observations show 
Sp,mg ~ 100 M a pc -2 dSolomon et alj 1 9871 ; iBolatto et alj2008l 



iHever et "all 2009), but pressure equilibrium implies that for galax- 
ies with Eg > 100 Mq pc -2 the GMC surfa ce density is set by th e 
disc surface density Egmc ~ E g (also see lHopkins et al. I l2012l) . 
We therefore adopt 



Egmc = max {E GMC , E g }. 



(31) 



Contrary to the distant encounters to which star clusters are 
often subjected, in the case of star-forming regions it is cru- 
cial to account for the spatially extended nature of the perturber. 
There are several corrections for tidal sho cking by extended mass 
distributions available in the literature jA guilar & White 1985; 
iBinnev & Tremai ne 1987; Gnedin et alj 19991 : IGieles et alj2 006b). 
Here, the approach of iGieles et alj d200 6bi) is adopted for mathe- 
matical simplicity and consistency - both in terms of the general 
derivation and the use of Plummer profiles. Their analysis shows 
that if the perturber is about an order of magnitude more massive 
than the perturbed object, their extended nature can be accounted 
for by adopting g ~ 1.5 in equation i30\ . Such a mass ratio is 
plausible because the GMC mass spectrum covers about two or- 
ders of magnitude, d istributed according to a power law with an in- 
dex of about —1 .7 ( Solom on et al.|[l987l ; lElmegreen & Falg arone 
1 19961 : 1 Williams & McK ee 1997). This indicates that most of the 
tidal disruption will typically come from those structures that are 
more massive than the region itself, but likely not by more than an 
order of magnitude. 

In TV-body simulations of the disruption of star clusters by 
tidal shocks, a fractio n / = 0.25 of the injected energy is used 
to unbind the cluster (Gieles et al. 2006b). This low efficiency is 
due to the excess energy that is carried away by escaping stars. 
However, in a dissipative medium such as the ISM the energy is 
not removed by a small number of high-velocity escapers, but in- 
stead it can be radiated away by the dissipation of turbulent en- 
ergy(f| The influence of turbulent energy decay on the value of / 
is addressed in Appendix lEl for a large set of Monte-Carlo exper- 
iments, which are used to model the stochasticity of the strength 
and time-separation of tidal perturbations. It is found that for those 
regions where t ccc ^ t, the fraction of the tidally injected energy 
that is used to unbind a region is about / ~ 0.7, independently of 
the environmental conditions. This relatively large fraction arises 
because regions are only disrupted if they happen to encounter a 
rapid sequence of tidal perturbations, thereby not allowing the en- 
ergy to be dissipated. Hence, a large fraction of the energy can be 
used to overcome the gravitational potential, and those regions that 
do not survive until time t are thus characterized by efficient dis- 
ruption. It is verified in Appendix |E]that this holds for the complete 
parameter range that is relevant for the presented theory (see £|2. 8 1 
and Table [T]below). 



9 This assumes that the kinematics of the star-forming region are domi- 
nated by the gas, i.e. the SFE is low. As will be shown in Figure[TJbelow, at 
the ages of interest this is indeed the case for those overdensities where the 
cruel cradle effect is important. 
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2. 7.2 The adiabatic correction 

In the harmonic approximation (cf. the adiabatic 

correction of equation ( 130b is given by ^ a d = exp (— 2<f>t), in 
which 4> t = tcnc/^diss, with t CYlc the typical time interval between 
encounters and tdiss the time-scale for the dissipation of turbu- 
lent energy. This assumes that the duration of the tidal shock is 
of the same order as the time interval between subsequent per- 
turbations, as should be expected due to the continuous nature 
of the ISM. The dissipati on time-scale is approximately given by 
jMcKee & Ostrikej|2007h 
R 



(Tint 



which in combination with equations i |13l > and {T3J becomes 



40 



-t s « 1.5 



/ Qyir \ 
V 1.3 ) 



-1/2 



te, 



(32) 



(33) 



implying that tdiss falls in the range 0.5£ff-3.8tfr for the val- 
ues of avir found in observations and numerical work. The 
encounter time-scale is d efined from simple kinetic theory as 
jBinnev & Tremaindl 19871) : 



= (4^7TnGMC0-g&max) 



(34) 



where tigmc is the number density of GMCs and 6 max = 
(48/7rnGMc) 1,/3 is the mean cloud separation. The GMC num- 
ber density can be specified by writing tigmc = Pism/Mgmc, 
with Mgmo the typical GMC mass. If one assumes that the char- 
acteristic GMC mass is set by the Jeans mass, it can be written as 
dKrumholz & McKeell2005h : 



Mgmc = 



4ft 4 



(35) 



Substitution of equations O, <[8l» and i35i into equation d34t then 
yields 



768 2 / 3 ^/G^ 



0.025^p /6 (G*pism)" 



1/2 



(36) 



Since at the solar galactocentric radius fio ~ 0.026 Myr -1 (as- 
suming a circular velocity of 220 km s" 1 ), equations (O and l !36t 
imply that t cnc ~ 2.5 Myr in the solar neighbourhood, which is a 
relatively low-density environment. The number of encounters in- 
creases with angular frequency, so at the radii where the bulk of the 
Galactic star forma tion occurs (between galac tocentric radii of 2 
and 7 kpc, see e.g. iMcKee & Williamd [19971) always t cnc < t. 
This justifies the earlier assumption that a star-forming region typ- 
ically suffers multiple encounters before t = 10 Myr. Finally, the 
combination of the final expressions for tdiss and t cnc now provides 
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This expression for <f> t as a function of a v ir and x now allows the 
specification of a a, which becomes 

1/2 ( x V /2 " 

Vio 4 J ' 

For gaseous regions, this ISpitzeJ dl987h formulatio n of the adia- 
batic correction is preferred over the lWeinbergMl994h correction, of 
which the power law tail at large 4> t is caused by stellar oscillations 
that do not necessarily apply to a turbulent medium. Since the pre- 
sented the ory approximate s gas-rich, star-forming aggregates lo- 
cally with Plum met J 1 9 1 lh potentials, which are parabolic in the 
centre, the harmonic approach is more appropriate. 



2.7.3 A critical overdensity for surviving the cruel cradle effect 

Substitution of equations J29l > and OOt into equation J27b now gives 
the disruption time-scale of star-forming regions due to the cruel 
cradle effect: 
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where R n is the half-mass radius of the star-forming region, and 
the second equali ty uses R n = 1.3-R p and R^/R 2 w 0.25 for a 
Plummer profile jGieles et alJl2006bl) . By writing t CC c = t, it is 
possible to derive a threshold overdensity a; CC c above which star- 
forming regions will survive the cruel cradle effect at least until 
time t, which is given by: 

124/^ sh f2E G Mc£ 
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where a d is a monotonically decreasing function of a; ccc , which 
indicates that equation J401 > has a unique solution. Because the re- 
lation is implicit, .t ccc has to be obtained by solving equation d40t > 
numerically. Regions with overdensities x < x CC c cannot form 
bound stellar clusters because they are disrupted by tidal perturba- 
tions due to the primordial environment before time t F°l A typical 
value for solar neighbourhood-type parameters is x CC c ~ 10 2 (see 
the discussion of Figure Q]in cj2.8| below). Also note that the pro- 
portionality x cco oc r2Q _1 E g " 1 implies x CC c oc M~ 1 . Because for 
Eg > E G mc we nave Egmc ~ Eg, high-surface density galaxies 

follow Xccc OC QQ^ 1 OC /Ojg^j. 

The influence of the cruel cradle effect on the CFE can now be 
quantified in two ways. Firstly, the fraction of the naturally bound 
part of star formation that survives the cruel cradle effect, i.e. / ccc 
from equation (TJl, is obtained by integration of the overdensity PDF 
of the ISM: 



/ccc — 



j(x)e(x)x(dp/dx)dx 
j{x)e(x)x(dp/dx)dx ' 



(41) 



where the numerator represents the bound structure that survives 
the cruel cradle effect, and the denominator denotes the part of 
star formation that results in bound structure. This quantity does 
not reflect the absolute contribution of the cruel cradle effect as 
equation J26b does for the naturally bound part of star formation, 
because the parts of the density range that are dispersed by both 
mechanisms overlap. As a result, / cce only gives the additional de- 
crease of the CFE due to tidal perturbations. To directly compare 
the individual contributions to dispersed star formation of naturally 
unbound star formation or the cruel cradle effect, it is thus useful 
to define the fraction of all star formation that survives tidal pertur- 
bations by the dense primordial environment: 



f 

J cc 



e(x)x(dp/dx)dx 
e(x)x(dp/ dx)dx ' 



(42) 



in which the denominator now represents all star formation. This 
definition allows one to determine the CFE if all star formation 



10 Note that this does not introduce a certain threshold density that can be 
used to define bound clusters observationally. Firstly because 3?ccc varies in 
space and time, but also because most regions that survive the cruel cradle 
effect (i.e. x > x CC c) are often intrinsically unbound, i.e. /bound < / CC c 
(see e.g. Fig.fjjand Table[3j- 
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Figure 1. The determination of the cluster formation efficiency (CFE). The 
solid curve shows the overdensity probability distribution function of the 
interstellar medium. Multiplication with the star formation efficiency gives 
the part that is converted into stars, which is represented by the dashed 
curve. The subset of stars that are formed in naturally bound stellar clusters 
is shown as a dotted line, of which the part that has an overdensity larger 
than a certain environmentally dependent, critical value (thick vertical line) 
survives tidal disruption by the dense star-forming giant molecular clouds. 
The CFE is given by the ratio of the grey-shaded region to the integral of 
the dashed curve. The figure shows the distributions for a parameter set that 
is characteristic for the Milky Way (see text). 

were intrinsically bound, giving an absolute measure for the dis- 
ruptive potential of the cruel cradle effect. As for /bound in equa- 
tion i26\ . the quantities that set / ccc and / ccc are the gas surface 
density E g , the Toomre Q parameter, and the angular velocity fl. 

2.8 The cluster formation efficiency 

The final expression for the CFE is obtained by combination of 
equations d26t and J41b , which provide the naturally bound fraction 
of star formation and the fraction thereof that survives the cruel 
cradle effect: 

p J°° -y(x)e(x)x(dp/dx)dx 

100% = - f "- f ^ = 7i e(x)x(dp/dx)dx ■ (43> 

This integral is to be evaluated numerically, and is illustrated in Fig- 
ure!]] It shows the overdensity PDF of the ISM xdp/dx before and 
after multiplication by e(x) and an additional factor 7(2;), respec- 
tively, as well as the environmentally dependent threshold overden- 
sity that ensures survival after the cruel cradle effect. The param- 
eter set is chosen to reflect the properties of the solar neighbour- 
hood, with E g = 12 M pc -2 , Q = 1.5, and Q = 0.026 Myr" 1 
being substituted in equation $6^ to obtain the Mach number and 
set the overdensity PDF. While higher Mach numbers decrease the 
value of the peak overdensity somewhat, they also lead to a much 
broader PDF. The curves that include the factors e and 7 nicely il- 
lustrate the three regimes of star formation given by equation | |24I >. 
In this example, star formation is still ongoing at age t for overden- 
sities x Ss 10°, whereas for x <; 10 5 it has been so efficient that 
an optimal SFE of e corc is reached. At intermediate overdensities, 
star formation is halted by feedback. Due to the cruel cradle effect, 
bound stellar clusters are only formed at high densities, i.e. either in 
the regime where star formation is halted by feedback or has been 
optimally efficient. The CFE equals the ratio of the grey-shaded 



Table 1. List of variables and their typical values. 



Variable 


Minimum 


Typical 


Maximum 


Quiescent 


Starburst 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


s s,o 


10° 


10 1 


10 5 


10 0.7 


10 3 


Q 


0.5 


1.5 


6 


3 


0.5 


Oo 


10 -2.5 


10- 17 


10 0.5 


10 -1.9 


10 -0.7 



Eg and are in units of Mq pc 2 and Myr , respectively. 



Table 2. List of parameters. 



Parameter Typical value 

(1) (2) 

<t>P 3 

Cfvir 1.3 

ten 3 Myr 

t 10 Myr 

(f>f h 0.16 cm 2 s -3 

€corc 0.5 

/ 0.7 

a 1.5 

4>nh 2.8 

s gmc 100 M © P c " 2 



region in Figure[T]to the integral of the dashed curve, which repre- 
sents the product e(x)xdp/dx and reflects the total amount of star 
formation. The figure illustrates that the cruel cradle effect and the 
naturally bound part of star formation both enable a fairly similar 
part of the overdensity range to remain bound. As such, the cruel 
cradle effect accelerates mainly the dispersal of unbound star for- 
mation, whereas its effect on the bound part of star formation is 
smaller. 

The presented theory of the cluster formation efficiency satis- 
fies the five criteria listed in i]2.1| 

(i) The fraction of star formation that is locally bound is formu- 
lated in a scale-free and continuous manner, and is based on the 
hierarchical nature of the ISM and the star formation process. 

(ii) The global bound fraction of star formation is obtained from 
a local criterion that is integrated over the entire overdensity spec- 
trum of the ISM. 

(iii) The formation of bound stellar clusters is considered as a 
continuous process, which allows the naturally bound fraction of 
star formation, the cruel cradle effect, and hence the cluster forma- 
tion efficiency to be calculated at any time t. 

(iv) By addressing the CFE at a specific time t, the theory ac- 
counts for three possible outcomes of the star formation process: it 
can occur on such a short time-scale that the gas is consumed before 
the onset of feedback, it can be truncated by large-scale supernova 
feedback, or it can still be ongoing at the moment of observation. 
This order reflects a sequence of decreasing density. 

(v) The effect of the tidal disruption of star-forming regions by 
the primordial environment is included and quantified at the time 
of observation t. 

The model depends on two or three main variables: (1) the 
gas surface density E g , (2) the Toomre Q parameter, (3) the angu- 
lar velocity Q or epicyclic frequency which are summarized 

If the CFE is determined globally for a low-redshift galaxy, Q and E g 
are roughly related according to equation j9}- 
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Figure 2. Influence of the adopted star formation law on the cluster formation efficiency. From left to right, the panels show the cluster formation efficiency T, 
naturally bound fraction of star formation /bound, fraction of bound star formation surviving the cruel cradle effect / CC c, and the fraction of all star formation 
surviving the cruel cradle effect f' CCB as a function of surface density and angular velocity. The displayed range is indicated in the legends at the top of each 
column, and the dotted contours denote fixed values as given by their labels. The solid lines mark the relation between S g and f2 from equation {5J along 
which most nearby galaxies are located, with dashed lines representing a scatter of a factor of three. The dash-dotted lines in the second column separate the 
subsonic (top left) and supersonic regimes. Top panels: using the e mpirically motivated star fo rmation law of equation jilt and Elmegreer] 120021) . Bottom 
panels: using the theoretical star formation law of equation [\2\ and Krumholz & McKee 1 2005). All other parameters are as in Tabled with Q = 1.5. 



in Table Q] together with their ranges. The 'typical' set of variable 
values reminisces the conditions in the solar neighbourhood and is 
broadly characteristic of disc galaxies. The 'quiescent' and 'star- 
burst' variable set|3 are chosen to reflect the typical extremes of 
star-forming galaxies, whereas the absolute extremes of the values 
themsel ves are listed in the 'minimum' and 'maximum' columns 
(see e.g. iKennicuttll 19981 : iBigiel e t al. 2008 for discussions of the 
observed parameter range). The small set of other relevant parame- 
ters for which constant, fiducial values can be adopted are listed in 
Table [2] These tables will be referred back to in the remainder of 
the paper when discussing examples of cluster-forming galaxies. 



3 THE VARIATION OF THE CLUSTER FORMATION 
EFFICIENCY 

Before making specific predictions for the CFE in different galax- 
ies and comparing them to observed values, it is relevant to ex- 
plore the parameter space and reach an understanding of CFE vari- 
ations across a range of galactic environments. In this section, we 
discuss the variation of the CFE as a function of gas surface den- 
sity and angular velocity, while each time varying one of the other 
model compone nts. In order o f appearance, these are the star for- 
mation law, the lToomrel d 19641) Q parameter, the time-evolution of 

12 Throughout the paper, 'starburst' is used to refer to galaxies with short 
(tdcpl 5s 300 Myr) gas depletion time-scales. The 'starburst' parameter set 
of Table |TJis chosen to be representative of high-density starburst galaxies, 
but does not rely on any criterion for gas depletion. 



the CFE, and some miscellaneous parameters that have a weak in- 
fluence on the CFE. The main results of the analysis in this section 
are summarized in ij3,5l 



3.1 General behavior and the influence of the star formation 
law 

We show the CFE as a function of the gas surface density E g and 
the angular velocity il in Figure [2] for the two diff erent star for- 
mation laws (Elmegreen 2002; Krumholz & McKee 2005) that are 
adopted in this paper. Because the behavior of the CFE is com- 
parable in both cases, we first focus on common features between 
both star formation laws, before pointing out the more subtle dif- 
ferences. The panels in Figure |2]do not only show the variation of 
the CFE, but also of the naturally bound fraction of star formation 
/bound, the fraction of bound star formation that survives the cruel 
cradle effect / ccc , and the fraction of all star formation that survives 
the cruel cradle effect / cco . We reiterate that the CFE is given by 
the product of /bound and / ccc , which of course also applies to the 
panels in Figure|2] 

The naturally bound fraction of star formation exhibits a pro- 
nounced increase with the gas surface density. This can be under- 
stood in terms of the Mach number A4 as defined in equation ^3). 
We mentioned in the discussion of Figure Q] that the Mach num- 
ber sets the shape of the overdensity PDF of the ISM. Because 
X a S g , the width of the overdensity PDF increases with the 
surface density, which enables a larger fraction of star formation 
to reside in the high-density range that leads to high local forma- 
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tion efficiencies of stars and clusters. This effect is large enough to 
overcome the decrease of the peak overdensity with increasing M . 
The value of /bound is remarkably insensitive to the angular veloc- 
ity for most of the parameter space, because in an equilibrium disc 
the decrease of the Mach number with £1 and the corresponding 
narrowing of the overdensity PDF is compensated by an increase 
of the average mid-plane density. This implies that while the over- 
densities may be smaller, the absolute density remains similar, and 
/bound only weakly changes with Q. The only exception occurs 
in the top left of the parameter space, at low surface densities and 
high angular velocities. In this region, /bound is insensitive to the 
surface density, but instead becomes a function of the angular fre- 
quency. The reason is that above the line Qo = 4.6E g ,2 the gas 
becomes subsonic and the overdensity PDF freezes. This removes 
the dependence on the Mach number, and implies that /bound only 
depends on the angular frequency through the mid-plane density of 
equation l(5j- 

The fraction of naturally bound (or total) star formation that 
survives the cruel cradle effect / ccc (/ ccc ) almost exclusively de- 
pends on the Mach number. This occurs because the critical over- 
density that is required to survive external perturbations scales as 
x ccc oc <7g oc M^ 1 due to the effect of gravitational focusing, 
while the PDF of the overdensity x also is a function of the Mach 
number only. Both dependences result in larger survival fractions 
for higher Mach numbers, and as a result the contours of constant 
/ ccc in Figure [2] mostly follow lines of constant Mach number, 
and mark a rather sudden transition between strong tidal disrup- 
tion and complete survival. Only at high surface densities a 'knee' 
appears in the contours, of which the location corresponds to the 
point where E g = Egmc = 100 Mq pc -2 . Towards higher densi- 
ties, the tidal perturbations grow because Egmc = E g , and hence 
the dependence of a; C cc on the surface density vanishes, implying 
Tccc oc Q.Q^ 1 oc Pjg^j. The upturn of the contours at high surface 
densities and angular velocities marks the point where the criti- 
cal overdensity for survival a: C cc coincides with the overdensity at 
which the free-fall time is so short that the SFE is optimal and all 
available gas is turned into stars (e = e C orc). In this regime, the tidal 
perturbations start to unbind structure that is entirely gravitation- 
ally bound. Finally, the difference between / cce and /c CC is mostly 
a horizontal stretch in which a smaller fraction of all star formation 
survives the cruel cradle effect than of the naturally bound part - as 
anticipated in the discussion of Figure[T] 

The total cluster formation efficiency T is given by the product 
of the second (/bound) and third (/ cco ) panels in each row of Fig- 
ure^ The cruel cradle effect provides a clear threshold for cluster 
formation at the low-density end of the parameter space, but when 
focusing on the region where most observed galaxies lie (solid and 
dashed lines), the CFE is mostly set by the naturally bound frac- 
tion of star formation. Only at very low and very high gas surface 
densities (cf. Table [TJ, the cruel cradle effect notably suppresses 
the CFE. The combination of all effects leads to a CFE that in- 
creases with the gas surface density, from a few per cent in qui- 
escent, low-density galaxies to a maximum of about 60% in high- 
density environments. It decreases again for the most extreme den- 
sities (Eg J; 10 3 Mq pc~ 2 ) due to the cruel cradle effect. 

The differences between t he two star formation laws are mi- 
nor. The iKrumholz & McKeej d2005h prescription leads to CFEs 
that can be up to 0.2 dex higher than for the em pirical star forma - 
tion law, but this falls within the error margins jKennicuttll 19981) . 
The main other di fference is a stronger dependence on the angu- 
lar velocity for the Krumholz & McKe^ d2005h star formation pre- 
scription than for the empirically motivated one. This is caused by 



an explicit dependence of their specific star formation rate per free- 
fall time on the Mach number, which causes the rate of star for- 
mation to increase with decreasing Mach numbers. That aside, the 
overall variation of the CFE for both star formation laws is very 
similar. In the following, we therefore assum e equation jilt unless 
it is cl arifying to also show the results for the Krumho lz & McKeel 
(2005) prescription of star formation. 

3.2 The influence of the Toomre Q parameter 

The variation of the CFE as a function of E g and Q is again shown 
in Figure[3] this time addressing the influence of the lToomrel i 19641) 
Q parameter, which indicates the stability of the gas disc. When E g 
and Q are specified, this parameter reflects the velocity dispersion 
of the gas as Q oc a s . As for Figure [2] the physics of the CFE 
are best understood by separately considering the naturally bound 
fraction of star formation and the fraction thereof that survives the 
cruel cradle effect. Again, the Mach number plays a central role. 

The Mach number scales as Mi oc cr g oc Q at fixed E g and Q, 
which leads to a slight decrease of /bound with increasing Toomre 
Q. This is easily understood in physical terms: at fixed E g and Q, 
a more stable gas disc implies a higher velocity dispersion, which 
lowers the mid-plane density as pism oc a^ 2 and causes the frac- 
tion of star formation occurring in high-density peaks to become 
smaller with increasing Q. The effect is largely offset by the broad- 
ening of the overdensity PDF due to the growth of the Mach num- 
ber, but this is insufficient to overcome the global density decrease 
entirely. 

Interestingly, the fraction of naturally bound star formation 
that survives the cruel cradle effect / ccc shows the opposite be- 
havior. As mentioned above, the overdensity PDF broadens as Q 
increases, while the critical overdensity for survival x ccc oc Q~ oc 
cr~ due to gravitational focusing. Both effects imply that the cruel 
cradle effect is more disruptive in unstable (low-Q) discs. When 
considering unstable, Q — 0.5 galaxies along the typical E g -Q 
relation (solid line in Figure[3}, the cruel cradle effect plays an im- 
portant (> 0.3 dex) role in setting the CFE for E g < 10 M pc -2 . 
It is also non-negligible for E g > 10 3 Mq pc -2 , where it limits 
the CFE by disrupting high-density stellar structure that otherwise 
would have remained bound due to the high local SFE. 

Combining the variation of /bound and / CC e with the surface 
density, angular velocity, and Q parameter, we see that the stabil- 
ity of the gas disc is instrumental mainly in setting the fraction of 
bound star formation that survives the cruel cradle effect. In partic- 
ular, unstable galaxies with high angular velocities (upper dashed 
line in Figure [3} have reduced CFEs due to tidal perturbations by 
the star-forming environment. The naturally bound part of star for- 
mation is much less affected by the disc stability and follows the 
same trends as discussed in ij3.ll Figure[3]also indicates the 'quies- 
cent', 'typical', and 'starburst' galaxies from TableQ] with CFEs of 
r = {4, 7, 59}%, respectively. As by the above discussion should 
be expected, the figure clearly shows that the CFE for the high- 
density set of variables is mainly affected by the cruel cradle effect 
(i.e. /ccc < /bound), whereas the 'quiescent' and 'typical' model 
galaxies reside in the part of the parameter space where a low nat- 
urally bound fraction of star formation determines the CFE (i.e. 

/bound ^ /ccc)- 

3.3 The time-evolution of the cluster formation efficiency 

By varying the time t at which the CFE is determined, it is possi- 
ble to address the time-evolution of the CFE across the ensemble 
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Figure 3. Influence of the Toomre 1 1964) Q parameter on the cluster form ation efficiency. Panels and parameters are the same as in Figure|2] this time adopting 
the empirical star formation law of equation jilt and Elmegreen 12002). From top to bottom, Q varies as indicated in the top left corner of each row. The 
diamonds, triangles, and squares indicate the 'starburst', 'typical', and 'quiescent' variable sets from Table[T] respectively. 




Figure 4. Time-evolution of the cluster formation efficiency. The panels indicate the 'quiescent', 'typical', and 'starburst' parameter sets from Table[T] which 
are also marked with symbols in Figure[3] The solid, dotted, dashed, and dash-dotted lines indicate the cluster formation efficiency, the naturally bound fraction 
of star formation, the fraction of bound star surviving the cruel cradle effect, and the fraction of all star formation surviving the cruel cradle effect, respectively. 
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of all star-forming regions in a galaxy, for the hypothetical situa- 
tion in which they are all coeval. This is a useful thought experi- 
ment for two reasons. Firstly, it can provide indications of how star 
and cluster formation proceeds. Secondly, it should be verified how 
sensitive the obtained CFE is to the time at which it is determined. 
Evidently, a high sensitivity could obstruct a reliable comparison to 
observations. 

Figure [4] shows the time-evolution of the CFE for the 'qui- 
escent', 'typical', and 'starburst' parameter sets from Table[TJand 
Figure[3] We see that the CFE increases with time until the onset of 
feedback at t sn — 3 Myr, which is a natural result of the ongoing 
collapse of gas into stars and bound structure while no gas is yet be- 
expelled^ When the first supernovae explode, star formation 
is halted on a global scale in regions with high SFE and density, 
while in low-SFE regions star formation can still continue because 
the feedback pressure is smaller than the ISM pressure. Dense stel- 
lar clusters thus collapse and complete the star formation process 
first, while the remaining star formation is only possible in those 
regions with increasingly low SFEs. Most of the residual, ongo- 
ing star formation is unbound, and further star formation therefore 
contributes very little to the CFE. This causes /bound to slowly de- 
crease for t > t sn , but otherwise it hardly exhibits any variation. 

The slight decline of the CFE is aided further by the criti- 
cal overdensity above which structure can survive the cruel cradle 
effect, which increases with time as is shown explicitly in equa- 
tion J40t . This has an obvious explanation: over a longer times- 
pan, a star-forming region is subjected to a larger number of per- 
turbations. The resulting decrease of / cco provides the largest time- 
evolution for ages t > t sn , and most prominently so in high-density 
environments (see the third panel of Figure|4j. 

Combining the above considerations, we conclude that the 
CFE is approximately time-independent - as long as t is chosen 
such that the densest regions that dominate the clustered part of 
star formation have been able to collapse and complete the star for- 
mation process (i.e. t > t S n)E3 This also explains why the CFE 
is insensitive to the adopted feedback energy input </>fb (see Ap- 
pendix|B). In the high-density regions where bound stellar clusters 
are born, either star formation is optimally efficient (e = e C orc) 
or the SFE is so high that feedback is sufficiently strong to halt star 
formation on a short time-scale. The details of how star formation is 
stopped elsewhere do not notably shift the balance between bound 
and unbound stellar structure. The only important time-evolution 
may occur in high-density (E g <; 10 3 M@ pc -2 ) galaxies, where 
the cruel cradle effect can lead to a decrease of the CFE with time. 



3.4 The influence of miscellaneous parameters 

The influence of some key parameters from Table [2] on the CFE is 
addressed in Figure|5] The first panel shows the impact of varying 
the characteristic time-scale t sn on which star formation is halted 
by feedback. The CFE increases with surface density due to the 
growth of the naturally bound fraction of star formation, until it 
reaches a peak and decreases at the highest densities due to the cruel 
cradle effect. The choice of t sn — 1.5-5 Myr most strongly affects 



13 When including radiative feedback, the CFE during this phase can be up 
to 0.5 dex lower. However, this occurs only in a small part of the parameter 
space, as is shown in Appendix [C] 

14 This is a natural choice since the CFE is established observationally 
by considering unembedded clusters, implying either a high SFE or some 
influence of feedback. 



the CFE at the lowest gas surface densities, where the variation 
can amount up to 0.6 dex. However, this spread decreases towards 
higher densities. Over the entire density range, shorter supernova 
time-scales yield lower CFEs. This is easily understood in terms of 
the presented model. A short t sn implies a lower SFE, which in turn 
implies that a smaller fraction of star formation will be bound after 
gas dispersal. 

The second panel of Figure [5] shows the impact of the GMC 
virial parameter a v ir on the CFE. For the empirically motivated 
star formation law, this influence is negligible because the virial 
parameter only enters as a second-order effect in the cruel cradle 
effect^ which in itself does not dominate the CFE over most of the 
densit y range. However, when adopting the Krumh olz & McKeel 
(2005) star formation law, the range of virial parameters a v i r = 
0.5-3 produces a variation that is comparable in magnitude to that 
brought about by setting t sn = 1.5-5 Myr in the previous panel. 
This is caused by the specific star formation rate per free-fall time 
of equation d!2t . which decreases with a v ir and thus leads to a 
lower SFE and CFE at higher virial parameters. 

The third panel of Figure[5]addresses the variation of the CFE 
with the maximum SFE e C orc- This also behaves in a predictable 
way - higher maximum SFEs increase the star formation rate in 
the high-density tail of the star formation spectrum, but at lower 
densities they also decrease the local bound fraction 7(e) = e /e COIC 
by requiring a larger SFE to attain the same 7. The latter effect ends 
up dominating the influence of e cor c on the CFE, but its magnitude 
is not as large as for the parameters in the previous panels because 
both effects compete. It is important to realize that the influence of 
this parameter is purely driven by the prior postulate that a region 
reaching e = e coro is bound in its entirety. As such, there is no real 
physics present in this panel. 

The fourth and last panel of Figure [5] shows the influence of 
the GMC surface density E G mc on the CFE. This quantity only af- 
fects the disruptive power of tidal perturbations, and only at surface 
densities E g < Eqmc because at higher densities Egmc = E g . 
Therefore, Eq^c affects the CFE only the low-density regime, and 
only if the cruel cradle effect plays a non-negligible role. For the pa- 
rameter set shown here (with Q — 1.5) this does not occur, but the 
small variations that are visible are indicative of the dependence in 
unstable discs (i.e. Q = 0.5). At low (E g < 10 Mq pc -2 ) surface 
densities, the CFE decreases with increasing Egmc due to the cor- 
respondingly stronger tidal perturbations. For stable galaxies this 
effect is entirely negligible. 

In summary, we see that the variation induced by changing 
the most important, 'fixed' model parameters is modest, and in all 
cases retains the characteristic increase of the CFE with the gas 
surface density. 

3.5 Summary of main model dependences 

The main conclusions of the parameter study in this section are as 
follows. 

(i) The CFE increases with the gas surface density E g of the 
host galaxy. This increase is driven by a growing fraction of all star 
formation that takes place in high-density regions that are capable 
of reaching high SFEs, thus evacuating the surrounding gas and 
becoming bound. At high surface densities, the CFE is inhibited 
by the tidal perturbation of star-forming regions due to the cruel 

lo Specifically, it sets the adiabatic correction of equation {38). 
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Figure 5. Influence of miscellaneous parameters on the cluster formation efficiency (CFE), adopting the parameters from Tableland taking Q = 1.5. The 
CFE is shown as a function of surface density S g only, obtained by writing the angular velocity as a function of S g (see equation (9) and the solid lines of 
Figures|2]and[3j. From left to right, the panels show the effects of the time of the first supernova t Bn , the GMC virial parameter ct v h-, the core SFE e C ore, and 
the giant molecular cloud (GMC) surface density Sgmc- The dotted, solid, and dashed lines indicate values of tsn = {1.5, 3, 5} Myr, cv v j r = {0.5, 1.3, 3}, 
fcoro = {0.25, 0.5, 0.75}, and Sqmc = {40, 100, 25 0} Mq pc -2 . In the second panel, the dash-dotted, long-dashed, and dash-triple-dotted lines indicate 
QVir = {0.5, 1.3, 3} for the Kramholz & McKee (2005) star formation law, which depends on the adopted value of c* v ir- 



cradle effect. This causes a peak CFE V ~ 50-70% at a density of 



10 M pc , beyond which the CFE settles at F 



40- 



60%. 

(ii) In the region of the E g -f2 plane where most galaxies reside, 
increasing the angular velocity SI very weakly affects the CFE, ei- 
ther yielding a modest increase (at intermediate surface densities of 
10 1,5 M Q pc -2 < E g < 10 2,5 M Q pc" 2 ) or decrease (at low sur- 
face densities of E g < 10 Mq pc -2 ). The only exception arises 
at high surface densities (E g > 10 2 ' 5 Mq pc -2 ) where the CFE is 
limited by the cruel cradle effect, of which the strong dependence 
on the ambient, mid-plane density pism oc Q 2 Q~ 2 implies that an 
increase of the angular velocity decreases the CFE. 

(iii) At fixed surfa ce density and angular velocity, the variation 
of the lToornrel d 19641) Q parameter reflects the variation of the tur- 
bulent velocity dispersion <r g . The most important consequence of 
increasing Q for the CFE is the corresponding decrease of the mid- 
plane density pism oc Q, 2 Q~ 2 , which leads to a smaller fraction 
of all star formation in high-density peaks and a weakened cruel 
cradle effect. The CFE thus grows with increasing disc stability. 
This also implies that the cruel cradle effect is most important in 
unstable discs. 

(iv) The CFE grows with time until the onset of feedback at 
t = t sn , after which it hardly changes. Dense clusters collapse and 
complete the star formation process first, after which the remaining 
star formation in lower-density regions contributes weakly to the 
global CFE. 

These characteristics persist when varying the miscellaneous pa- 
rameters from Table[2]that are kept fixed throughout the rest of this 
paper. 



4 COMPARISON TO OBSERVED GALAXIES 

Having addressed the behavior of the CFE across the parameter 
space, it is possible to compare modelled CFEs to observations and 
interpret why potential differences might arise. We first do this for 
the broad sample of galaxies for which CFEs have been determined, 
and then do the comparison on a galaxy-by-galaxy basis. 



4.1 The empirical relation between the cluster formation 
efficiency and the star formation rate density 

The CF E has been determined observationally in several recent 
studies jGod dard et al.l 1201(1 ISilva-Villa & Larsenl |2010|. 1201 ll; 
lAdamo et al]|20ldl2fJl~illCook et al.ll2012l) . which together cover 
a galaxy sample containing dwarf, spiral, and starburst galaxies. In 
most cases, these studies compare the age distributions of bound 
clusters to the star formation rates or histories of the host galaxies. 
The ratio of both quantities then gives the CFE. The observational 
results typically cover a range of V = 1-50%, which is comparable 
to the outcome of our parameter survey in !j3] 

Particular attention has been given to how the CFE varies 
with the star formation rate density of the host galaxy Ss f r- Th e 
analyses by iGoddard et al.l OQlol) ; ISilva- Villa & Larserj OOlll) ; 
lAdamo et alj q20 1 lh have shown that the CFE increases with Esfr 
in a way that is similar to the correlation with the gas surface den- 
sity E gaa th at results from o ur model. Both den sities are related 
through the Kennicutt (1998) formulation of the ISchrriid t ( 1959) 
star formation law as 



Espr,o = 2.5 x 10~ 4 E 1 ' 4 



(44) 



g,o> 

where E S pr,o = E S fr/M yr" 1 kpc -2 and E g , = 
Eg/Mg pc -2 . This relation enables us to convert the model gas 
densities to star formation rate densities and directly contrast theory 
and observations. As a first comparison, the CFE can be calculated 
for a typical parameter set. 

The observed CFEs are shown in Fi gure [6] for the com- 
bined g alaxy sample ofrGoddard et al](l2010h . lSilva- Villa & Larsenl 
j201 lh . lAdamo et alj fcoi ll) . and lCook et alj 00121) . together with 
the modelled T — Esfr relation for Q — 1.5 and with specified 
as in equation l|9j. The agreement between theory and observations 
is striking, and gains significance considering the modest variation 
of the model result with the parameters in Table|2](see [J3]and |7JJ , 
which means that the overall trend of the CFEs in Figure[6]is a solid 
outcome of the theory presented in this paper. We thus obtain a sat- 
isfactory explanation of the observations. As shown in !j3] the CFE 
increases with the surface density because at higher densities (and 
Mach numbers) a larger fraction of the density spectrum of star for- 
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Figure 6. Cluster formation efficiency (CFE) as a function of the star 
formation rate density. Symbols d enote observed galaxi es with lcr error 
bars and indicate th e samples from lGoddard et a l] J2010L blue diamon ds), 
lAdamo et all 1201 ll green triangles), and ISilva- Villa &Larsenl l201lL red 
squares indicate their Tmdd)- The black cross indicates the integrated CFE 
of all dwarf galaxies from the sample of lCook et al.l 120121) , with a surface 
density range indicated by the horizontal error bar. The solid curve repre- 
sents the modelled relation for the 'typical' parameter set from Table[T] but 
can vary for different galaxy and ISM properties. The dashed curve shows 
the fit of equation 1451 , and the dotted line represents the original fit by 
iGoddardet all fcOld) . 

mation is pushed into the regime with short free-fall times, enabling 
it to reach high local SFEs and allowing it to become gravitationally 
bound. 

For practical applications, the dashed line in Figure [6] repre- 
sents a good fit to the model, and is given by 

T = (1.15 + 0.6Egp R 4 + O.OSEgpVo)' 1 x 100%. (45) 

This fit is for one particular, 'typical' parameter set and should 
therefore mainly be used for rough estimates. As explained be- 
low, this parameter set is not accurate for high-density galaxies 
(Esfr,o ii 1) and may overestimate their CFE by up to a factor 
of two. It also does not account for differing properties of indi- 
vidual galaxies, which can cause a similar variation. Recall that the 
star formation rate density can easily be replaced by the gas surface 
density using the Kennicutt-Schmidt law of equation J44b . 

One word of caution is in place. The observational CFEs are 
uncertain up to a factor of two or thre e due to assumptions in their 
derivations (see iGoddard et al.ll201ol for an extensive discussion). 
For instance, they rely on the extrapolation of the observed cluster 
mass function down to some minimum cluster mass below the de- 
tection limit, which is affected by uncertainties on the power law 
slope of the cluster mass function and the lower mass limit. The 
estimates are also sensitive to errors in the age and mass fits with 
simple stellar population models, un certainties in the metallicity , 
and possible selection effects. Finallv. lSilva- Villa & Larserj (2011) 
determined the CFE using two different methods. We adopt their 
MDD sample, but the trend only changes mildly when using their 
MID sample instead - in fact, the agreement between theory and 
observations slightly improves when taking the mean of both meth- 



ods. The error margin resulting from all of the above is realistically 
shown for at least a few of the formal error margins in Figure[6] but 
in some evident cases the errors are strongly underestimated. The 
aforementioned factor of two to three is more reasonable. 

Next to the above considerations, there are also some physical 
grounds for the variation between the individual galaxies in Fig- 
ure|6]and the model relation. Firstly, the line only shows the model 
for one particular parameter set, and is based on the scaling relation 
between Q. and E g of equation (|9)- The scatter around that relation 
is about 0.5 dex, implying that it can influence the CFEs of indi- 
vidual galaxies to a statistically significant degree. Secondly, high- 
density galaxies often undergo starbursts and are therefore typically 
less stable than isolated disc galaxies, which means that there is a 
trend of Q with surface density. This implies that at high densities 
the cruel cradle effect is more prominent than shown by the model 
here (cf. §3.2\ . Thirdly, at low densities galaxies deviate from the 
simple power law Schmidt law that we assume to convert our model 
densities to star forma tion rate densities, having a larger spread at a 
given gas density (e.g. Bigiel et al. 2008). S tochasticity should also 
play a larger role at low densities Jda Silva et ail 2012). and hence 
the increased spread of CFEs in the low-density regime is not very 
surprising. We address most of these considerations in i j4,2| below, 
where we estimate the CFE of each galaxy individually. 

A final remark concerns the inclusion of radiative feedback, 
which is explored i n Appendix Icl A simple inclusion of radiat ive 
feedback (following iThompson et al]|2005l: iMurrav et alj|2010l) in 
the theory of this paper predicts a jump of the CFE around E g ~ 
200 Mq pc -2 , due to the increase of the GMC density above such 
densities. This gas surface density corresponds to Espr,o ~ 0.5 in 
Figure[6] where a vague hint of such a jump might be present. It is 
clear that the uncertainties on the observations are far too large for 
any conclusive statements, and unfortunately the absence of galax- 
ies around Espr,o ~ 0.1 is not very helpful either. However, this 
does show that better observational estimates of the CFE could aid 
in establishing which feedback mechanisms drive turbulence and 
star formation on galactic scales. 

4.2 Cluster formation efficiencies in individual galaxies 

The excellent agreement between theory and observations shown in 
Figure[6]warrants a more detailed comparison. Given observed val- 
ues for the gas surface density, the angular velocity, and the Toomre 
Q parameter, it is possible to calculate the theoretical CFE for a 
given galaxy. There are a number of studies in w hich the actual 
CFE was estimated for a small sample of galaxie s (Goddard et al. 
20ldJGielesj|2010l: lAdamo et al]|2010l : ISilva- Villa & LarserJl20lIF 
Adamo et alj|201 II) . This enables a quantitative test of the theory 
presented in q2.3H2.8l 



4.2.1 The observed properties of dwarfs, spirals and starbursts 

The sample of galaxies for which a comparison of the theoretical 
and observed CFEs is carried out is listed in the left-hand columns 
of Table[3] together with the input parameters for each galaxy. The 
gas surface densities are obtained by averaging the tota l gas mass 
within the optical radius R25 to ensure consistency with lKennicuttl 
( 1998). In the cases in which this is not accurate because the gas 
disc is smaller than R25, the listed surface density is averaged over 
the area covered by the gas only. In the cases of the starburst galax- 
ies NGC 3256 and Haro 11, which should exhibit substantial in- 
ternal variation of the surface density (and therefore of the CFE), 
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Table 3. Properties and cluster formation efficiencies of galaxy discs and starbursts. 
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Col. (1): Galaxy identifier. Col. (2): Logarithm of the gas surface density in M© pc — 2 . Col. (3): Angular frequency in (100 Myr)~ 1 . Col. (4): Velocity 
dispersion of the gas disc in km s — 1 . Col. (5): Toomre Q parameter as defined in equation (8). Col. (6): Theoretical bound fraction of star formation. 
Col. (7): Theoretical surviving fraction of bound star formation due to the cruel cradle effect. Col. (8): Theoretical surviving fraction of all star forma- 
tion due to the cruel cradle effect. Col. (9): Theoretical cluster formation efficiency (CFE) in %. Col. (10): Observed CFE in %. Col. (11): Logarithmic 

dif ference between theor etica l and observed CFEs. Col. (1 2): References as follo ws. 

mlChemin et alj feOOrj) (2)ISilva-Villa & Larserj 1201 ll) (3) iKennicunl 1 19981) (4)ICombes & Becquaerll 1 19971) (5) iGielesj feOld). (6 ) |Rvder et alj 
ll995t). (7) IStil & Israej 120021). (8) iGoddard et alj l2010l). (9) ISakamoto et all 120061). (10) ISwaters et alj jl999h, (11) ISchuster et alj 120071), (12 ) 
iLunt 



Lundgr enet aljj2004) q3) IWalsh et aljj2002h (14) lpicaire et aljl2008l) ri5)lAdamo et alj 1201 ll). (16 ) lBergvall et al.ll200(l) (17)|6stlin et alj«199' 



(18)l6stlin et alJfeOOlh. (19)lwilke et aljl2004l). (20)IStanimirovic et aljfe004h. (21)IWong et alj J2009l). ( 22)1 Alves & Nelson! l2000|). (23)IWolfire et alj 
120031) . q4) lHeiles & Troland J2003l) . (25l lChemin et alj J2009h . (26) teraun et all2009h . (27) lYun et alj fl993h . f28) foownes & Solomorj Jl998h ' 

a If no observed velocity dispersion is available, a fiducial velocity dispersion of 6 km s~ _1 teennicutl 19891) is adopted in the calculation of equation (8) 
and the value for Q is shown in parentheses. 

" This is a lower limit because th e molecular gas mass has not been estimated. 

c Average of the values found by Goddard et al. 1 2010> for the nuclear starburst) and Silv a- Villa & Larser] 1201 ll for the spiral arms). 
^ Because no observed s urface de nsity is available , this value is derived from the average of the star formation rate densities listed in 
ISilva- Villa & Larsenl 1201 ll) . using the lKennicutil 1 19981) star formation law. 

e The listed values are exclusive to the circumnuclear starburst. The gas surface density accounts for molecular hydrogen only. 



the surface density is determined for the same area as for which 
the observed CFE has been established. For NGC 3256, most of 
the gas mass is confi ned to the inner 1.7 kp c, but the disc contin- 
ues to at least 3 kpc jSakamoto et alj[2006j) at 20% of the central 
density. A total gas mass of 3.3 x 10 9 Mq within 3 kpc is there- 
fore adopted. For Haro 1 1 , most of the gas mass resides within 6.5" 
or 2.6 kpc, which c ontains about 2 x 10 9 M dOstlin et all 19991 ; 
iBergvall et al]2000h . Whenever only angular sizes are given by the 
literature sources, the distance from the NASA/IPAC Extragalactic 
Databas43 (NED) is used to convert to physical dimensions, unless 
the distance is given by the original source. For the few galaxies 
where R25 is not given by the corresponding literature sources, the 
dimensions from NED are used to calculate the surface densities. 

The angular velocities are evaluated from the galaxy rotation 
curves at 0.5i?25, which corresponds to Q, = 4 7r/t d vn (w i th tdv n 
the dynamical time-scale) for the galaxies from Ken nicuttl Jl998l) . 



In those galaxies that are not taken from the Kennicutt ( 1998j) sam- 
ple and in which less than the full area of the galaxy is covered by 
either the gas disc or the observed CFE, the angular velocity is de- 
termined at half the maximum radius. Given a surface density and 
angular velocity, the Toomre Q parameter is calculated from the 
velocity dispersion of the gas. Like the angular velocity, the veloc- 
ity dispersion (obtained from H2, CO or Hi observations, in that 
order of preference and depending on availability) is evaluated at 
0.5/?25 when a radial profile is provided by the literature source. 
Whenever the velocity dispersion is only available for certain ra- 
dial intervals, an area-weighted average is listed in Tabled If no 
value is prese nt in the literatur e at all, a fiducial velocity dispersion 
of 6 km s~ JKennicuttll 1989T) is adopted in the calculation of Q. 
It should be noted that the Milky Way values are representative of 
the solar neighbourhood instead of the entire disc. Finally, we cor- 
rect for the presence of spiral arms by halving the Q parameters 
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Figure 7. Comparison of observed cluster formation efficiencies r o b s to the 
values calculated with the theoretical model r tn for the galaxies of Table[5] 
The solid line indicates the one-to-one relation, while the dashed and dotted 
lines indicate deviations by factors of two and three, respectively. 



as in Krumholz & McKee (20051) for all galaxies except the Milky 
and the circumnuclear starbursts of M82 and Arp 220. 



4.2.2 Comparison of observed and theoretical CFEs 

Using columns (2)-(5) of Table [3] for each galaxy we calculate 
the naturally bound fraction of star formation /bound, the fraction 
of bound star formation surviving the cruel cradle effect / cce , the 
fraction of all star formation surviving the cruel cradle effect / cce , 
and the resulting cluster formation efficiency r tn - These quantities 
are listed in Table [3] including the logarithmic difference between 
the observed and theoretical CFEs (log rth/r o bs)- The comparison 
is also shown graphically in Figure|7] 

Given the inhomogeneity of the galaxy sample, the assump- 
tions that were made for certain unknown observables, and the trou- 
blesome process of determining the CFE observationally, the agree- 
ment between theory and observations is remarkable. As discussed 
in £|4. 1 1 a spread of at least a factor ~ 2 should be expected con- 
sidering the uncertainties in the comparison. This is approximately 
consistent with the scatter of up to a factor 2.5 around the one-to- 
one relation in Figure [7] This includes the statistical errors on the 
observed CFE and the model input variables {E g , Q, fl} (also see 
flTTland Figure[[0]!. 

Contrasting the listed values of /bound and / ccc , it is clear 
that the contributions of both mechansisms to the CFE are non- 
negligible. However, they favour the survival of comparable over- 
densities (see Figure QJ. This overlap implies that the additional 
decrease of the CFE due to the cruel cradle effect typically does 
not exceed 30%, as is indicated by the values of / ccc . It has been 
mentioned throughout this paper that tidal perturbations are most 



17 The solar neighbourhood is located right between two spiral arms (see 
e.g. lPorte gies Zwart ^Talj2010l) . 



important for high-density galaxies such as NGC 5236, but per- 
haps surprisingly the results in Table [3] show that it is also relevant 
for low-density galaxies like NGC 45, NGC 4395, and NGC 7793. 
The importance for low-density galaxies is caused by their rela- 
tively low Mach numbers, which imply that the influence of tidal 
perturbations is enhanced due to a larger influence of gravitational 
focusing. In galaxies with intermediate surface densities, the CFE 
is dominated by the naturally bound fraction of star formation. 

Table [3] also includes a prediction for four galaxies that are 
obvious or interesting candidates for future observational estimates 
of the CFE. These include M31, which i s currently undergoin; 
an extensive surve y by the PHAT team ( Dalcanton et alj 1201 
I Johnson et afll2012h for the Hubble Space Telescope treasury pro 
gram. The low gas surface density of M3 1 leads it to have the low- 
est CFE (T = 1.7%) of the entire galaxy sample in Table[3] which 
is a prediction that can be tested using the PHAT data. Another clear 
candidate for a future observational determination of the CFE is 
M81, for which we also obtain a relatively low CFE of T = 3.5%. 
The other two listed galaxies are M82 and Arp 220. Establishing 
the CFE in those galaxies will be quite a challenge due to the high 
extinction and/or the large distance. Nevertheless, these are extreme 
starburst environments, and the model predictions are correspond- 
ingly illustrative. In both galaxies, (almost) all stars are formed in 
naturally bound systems due to the extreme density and the result- 
ing efficient star formation. However, the CFE ends up being much 
lower than 100% as it is limited by the cruel cradle effect, which in 
these cases leads to F = 30-40%. 



THE SPATIALLY RESOLVED CLUSTER FORMATION 
EFFICIENCY 



The theoretical framework of i]2.3M2.8l can be used to calculate the 
radial profile of the CFE within a galaxy. This requires the gas sur- 
face density profile, rotation curve and gas velocity dispersion of 
the galaxy in question. For several nearby galaxies such informa- 
tion is available, and some illustrative examples are given in this 
section. 

We consider the radial variation of the CFE in the Milky Way, 
M31, and for a generic exponential gas surface density profile, 
which happens to give a good representation of M51. The input 
variables for the different galaxies are as follows. 



(i) For the Milky Way, we use the IWolfire et ail d2003l) model 
for the ISM, including a factor of 1 . 4 to account for the pres ence of 
helium, and adopt a g = 7 km s _1 dHeiles & Trol and 2003) to cal- 
culate Q. The angular velocity Q follows from a constant circular 
velocity of v c = 220 km s~ . 

(ii) For M31. IChemin et all (2009) provide the Hi surface den- 
sity profile, which is us ed here with a corr ection factor of 1.3 to 
account for neutral gas dNieten et alj 20061). The adopted circula r 
velocity is v c = 250 km s" 1 teraudl 199lt IChemin et af]|2009h. 
and t he gas velocity dispersion is a g = 8 km s _1 (Bra un et al.1 
l2009h , 

(iii) A generic exponential gas surface density profile is given by 
Eg(i? gc ) = S g (0) exp (— Rgc/Ra). This is a reas onable descrip- 
tion o f the ISM in M5 1 when taking Ra = 3 kpc dSchuster et al.1 
120071) . and therefore we set the other parameters to values that 
are appropriate for that galaxy. We consider central surface den- 
sities of E g (0) = {10,30,100} M Q pc" 2 , where the highest 
value is representative of M51. T he circular velocity is taken to 
be v c = 200 km s" 1 dRandl 19931) and is converted to an angular 
velocity profile using the M51 scale radius of Rd = 3 kpc. The 
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Figure 8. Variation of the cluster formation efficiency (CFE) with galactocentric radius i?, gc . Solid lines show the result for the empirically motivated star for- 
mation law of equation {TT} and Elmegreen 1 2002), while dashed lines are based on the theoretical star formation law of equation d!2t and Krumholz & McKee 
120051) . Dotted lines represent the gas surface density profiles, which in each case are normalized to the outermost CFE for the empirical star formation law. 
See text for the other adopted parameters of each galaxy. Left panel: CFE profile for the Milky Way. Middle panel: CFE profile for M3 1 . Right panel: CFE 
profile for a generic exponential gas disc, with the galactocentric radius expressed in units of the scale radius R^. From bottom to top, the curves indicate 
central gas densities of £ g (0) = {10, 30, 100} Mq pc -2 , respectively. The two CFE profiles with E g (0) = 100 Mq pc -2 represent M51, which has 
i?d = 3 kpc. 



adopt ed gas velocity dispersion is a g = 8 km s 1 dSchuster et al.l 
l2007h . 

As in A4.2I we corre ct for the presence of spiral arms by hal ving 
the Q parameters (cf. lBalbus|[l988l ; lKrumholz & McKeell2005l) . 

The resulting CFE profiles are shown in Figure [8] using the 
two different star formation laws used in this paper (cf. fl3.lt . 
As should be expected from our earlier analysis, the CFE profiles 
largely follow the surface density profiles of the galaxies, especially 
at low densities where T oc E g . In the Milky Way, this gives a 
rather flat CFE profile, with a shallow minimum in the solar neigh- 
bourhood (R sc = 8.5 kpc) and modest peaks at 5 and 13 kpc. The 
overall variation across a 12 kpc radius interval is about 0.15 dex. 
The typical CFE throughout the Milky Way di sc is F ~ 6% (empir- 
ical s tar formation law) or F ~ 12% (for the Krumh olz & McKed 
120051 star formation law). M31 exhibits a pronounced variation of 
the CFE, which again traces fluctuations of the gas surface density 
with galactocentric radius. Typical CFEs range from V — 1-5% 
for the empirical star formation law, with the peak values reached 
between 8 and 15 kpc. Note that we did not account for any ra- 
dial variation of the gas velocity dispersion, which could lead to a 
modest (^ 0.2 dex) increase of the CFE towards the centre of the 
galaxy. 

The exponential surface density profiles imply a much 
stronger trend with galactocentric radius than in the previous two 
examples. As shown in the right-hand panel of Figure [8] the CFE 
covers a broad range, hitting T — 30% in the central parts of M51 
to r = 3% at R gc = 10 kpc ~ 3Rd- Like in the preceding cases, 
the CFE profiles shown here assume a flat rotation curve and no 
radial variation of the gas velocity dispersion. Real galaxies exhibit 
solid-body rotation within the central few kpc (i.e. v c oc R sc and 
f2 ~ constant) and have velocity dispersions that rise towards 
their centres. Of course, this can be accounted for in choosing the 
model input variables, but both effects induce only & 0.2 dex devi- 
ations of the CFE in the galaxy centres with respect to the relations 



shown in Figure [8] In general though, it is questionable to what 
extent the model holds in such environments. Cloud-cloud colli- 
sions might be playing an important role in determining whether 
star formation can proceed and structure can remain bound. Even 
more importantly, recent sub-millimeter observations of the central 
molecular zone (CMZ) of the Milky Way suggest that the tracers 
of dense gas are overabundant by u p to a factor of a hundred com- 
pared to the star formation tracers (Longmore et al. 2012a). This 
might indicate that star formation in the Galactic centre does not 
obey the star formation laws that apply to galaxy discs. While this 
does not directly affect the CFE, it seems plausible that the clus- 
ter formation process would be affected as well. Either way, galaxy 
centres clearly provide important tests for current ideas about star 
and cluster formation. 



Figure [8] shows that the CFE directly traces the gas surface 
density for surface densities E g ^ 20 Mq pc~ 2 . This changes 
at higher densities, as is shown by the M51-like CFE profile. The 
reason is that for E g <; 20 Mq pc~ 2 , the increase of the CFE 
with the gas density becomes shallower than linear (see e.g. Fig- 
ure[5]l. Another effect would potentially be that the central regions 
of the galaxies have enhanced angular velocities compared to the 
gas densities and therefore lie above the E-S7 relation from equa- 
tion ^ that holds for galaxies as a whole. They only evolve towards 
that relation at increasing galactocentric radius. However, for disc 
galaxies this occurs in a regime where the CFE is virtually insen- 
sitive to SI (see fl3.lt . and hence the influence of a varying angular 
velocity is negligible. More extreme spatial variations of the CFE 
should be present in starburst galaxies like M82 and Arp 220. The 
central regions of these galaxies have CFEs that are limited by the 
cruel cradle effect, whereas in their outskirts it should be set by the 
naturally bound fraction of star formation. Starburst galaxies are 
therefore ideal targets for exploring the boundaries of the presented 
theory and for understanding the CFE in more detail. 
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6 STAR CLUSTER FORMATION OVER THE COURSE 
OF A HUBBLE TIME 

The conditions under which stellar clusters formed have changed 
notably with cosmic time. As shown in the preceding sections, 
the efficiency of cluster formation is sensitive to the proper- 
ties of the host galaxy. The cosmic star formation rate den- 
sity peaked at r e dshift z ~ 2-3 (e.g . iMadau et all 1 19961. fl998l; 
Lilly et alj|l999l: IPerez-Gonzalez et al.ll2005l ; iHopkins & Beacoml 



20061 : iReddv & Steidell l2009h . reflecting that galaxies at these 



times were denser and more gas-rich than in the nearby Uni- 
verse. Indeed, observations of high-redshift galaxies show that 
the gas densities (10 2 -10 5 M© pc -2 ) and star formation rates 
(20-3000 Mq yr _1 ) commonly reached extreme values in the 
early Universe that are presently rare (e.g. [ Chapman et alj 20041; 
Solomon & Vanden Boutl 120051: iDaddi et al.1 120071 ; iGenzel et al.1 



20101) . Local, extreme starburst galaxies like the merger remnants 



Arp 220 and NGC 624 are the nearest examples of galaxies reach- 
ing similar conditions (Rieke et al. 1985; S coville et al .11 1 99 lh . Re- 
markably though, many of the high-redshift star-forming galaxies 
appear to be 'normal' star-forming (albeit unstable) disc galaxies, 
of which the high star f ormation rates are simply caused by their 
substantial gas reservoir telmegreen & Elme green 2005; Er b et al.l 
20061: IGenzel et al.l2008l ; lDaddi et alj20ld : lElmegreen et al.l2009l : 
Tacconi et al.ll2010l) . With such observational input at hand, it is 



possible to address the CFE under the high-redshift galactic condi- 
tions. 

We use the cluster formation theory of this paper to predict 
the CFEs of some ~ 150 galaxies in the redshi ft range z = 0-3 . 
The sample of nearby disc galaxies is taken from Ken nicuttl jl998h . 
who tabulated their gas surface densities E g and dynamical times 
tdyn- The angular velocities follow as fi = 47r/td yn , and we as- 
sume a constant velocity dispersion a g — 6 km s _1 to calculate 
Q. This suffices to estimate the CFE, which is a function of E g , 
Q, and SI in the presented model (see §2.8\ . For the high redshift 



galaxy population, the combin e d gala x y sample of Tacconi et al.1 
(2008), iForster Schreiber et al.1 J2009I) . flacconi et alj ( l2010h . and 



Genzel et alj feOlOl) is used. Surface densities are derived from the 
gas masses M gas and Ha radii Rna as E g = 0.5M ga s/7r-RH a , 
and the angular velocity is calculat ed from t he cha racteristic veloc- 
ities V as n = 2V/Rn a . The lToomrel Jl964l) Q parameter is 
chosen based on the gas depletion time-scale £dc P i = M gas /SFR. 
If it exceeds 300 Myr, the galaxy is assumed to be stable and we 
adopt Q — 1.5. For fdepi < 300 Myr, the galaxy is considered to 
undergo a starburst and we assume Q = 0.5. While this separation 
is somewhat arbitrary, it does enable a rough separation between 
stable and unstable systems. 

The calculated CFEs are shown in Figure [9] The left-hand 
panel shows the CFE in the Ssfr-T plane of Figure [6] this time 
covering a larger range to account for the extreme star formation 
rate densities. For reference, two model lines for Q — 0.5 and 
Q = 1.5 are includedEl The predicted CFEs of nearby disc galax- 
ies cover a range that is very similar to that of the observations 



18 This velocity can either represent the circular velocity or the veloc- 
ity dispersion, depending on whether the galaxy is rotation- or dispersion- 
dominated. 

19 Note that the symbols and lines are both theoretical predictions. As 
such, the 'scatter' of the predicted CFEs around the theoretical lines is com- 
pletely unrelated to any (dis)agreement between model and observations. It 
is only caused by (1) deviations of individual galaxies from the S-Q re- 
lation of equation {9)> and (2) deviations of individual galaxies from the 



in Figure [6] increasing from Y ~ 1% at low star formation rate 
densities to F ~ 20% at the high-density end. The high-redshift 
star-forming disc galaxies represent a natural continuation of this 
trend to higher densities and CFEs, increasing from V ~ 30% to 
a saturation at F ~ 70%. These discs are at the point of perfect 
balance between a high naturally bound fraction of star formation 
due to their high densities, and a sufficiently low density (and ade- 
quate stability) to have star-forming regions not suffer too strongly 
from tidal perturbations. The high-redshift starburst galaxies are 
in a different regime, where in some cases the CFEs are strongly 
suppressed by the cruel cradle effect. This is similar to the earlier 
prediction for the CFE in Arp 220 (see E|4.2| and Table f3]l. For a 
given galaxy sample, the highest CFE should thus be attained at a 
star formation rate density of Ssfr = 1-10 M yr _1 kpc -2 or a 
gas surface density of S g ~ 10 3 Mq pc 2 . 

The evolution of the CFE with redshift is shown in the right- 
hand panel of Figure [9] The galaxy sample is divided into galaxy 
discs and starburst galaxies according to the above criterion of the 
gas depletion time-scale, and binned in Az = 1 redshift bins to 
show trends of the mean CF E. To obtain a z = data point for 
starburst galaxies, we add the I Kennicuttl jl998h sample of nearby 
circumnuclear starbursts, using the same analysis as for the local 
disc galaxy sample and assuming Q — 0.5. The model predicts 
that the CFE in starburst galaxies has been relatively constant dur- 
ing the past ~ 11 Gyr, whereas the CFE in disc galaxies expe- 
rienced a strong decrease with cosmic time. Both lines cross at 
z ~ 1, with starbursts being the most efficient cluster factories at 
lower redshifts, and disc galaxies having a slight edge at higher red- 
shifts. Note that the error bars do not indicate the statistical signif- 
icance of the plotted mean value, but instead indicate the spread of 
the underlying sample. The differences between the disc and star- 
burst galaxy samples are statistically significant, particularly for the 
z — sample. This has been verified by running a Kolmogorov- 
Smirnoff test. For the hypothesis that the CFEs of both popula- 
tions are drawn from the same parent sample, we obtain p-values 
of Pks = {0.004, 0.148, 1.73 x 10~ 8 } for z = 2-3, z = 1-2, and 
2 = 0-1, respectively. 

The restricted sample of 166 galaxies that is used for the 
right-hand panel of Figure [9] is not necessarily representative of 
the global ga laxy popul a tion. T he local Universe sample relies ex- 
clusively on I Kennicuttl J 1998h . of which the galaxy selection is 
mostly based on the availability of CO and Hi maps, yielding a 
galaxy sample containing a mix of field and Virgo cluster galax- 
ies. The sample should therefore be free of strong selection bi- 
ases. The high-redshift sample contains intermediate-to-high-mass 
(M s tar > 2 x 10 9 M ), star-forming (SFR > 6 M yr" 1 ) galax- 
ies, with a large contribution from the SINS survey dGenzel et al.1 
l2008l : lForster Schreiber et alj2009h . These ranges are chosen to in- 
clude a broad range of main-sequence, star-forming galaxies and 
hence do not impose strong restrictions on the galaxy population 
at z = 1-3 other th an being limited to the bright end of the lumi- 
nosity function (e.g. lForster Schreiber et al] |2009). For our sample 
of nearby main-sequence galaxies (see Table[3}, there is no appre- 
ciable correlation of the CFE with galaxy mass, and hence extrap- 
olating the results for the SINS sample to fainter galaxies seems 
reasonable. Nonetheless, it should be kept in mind that any low- 
density (i.e. ~ low-SFR) galaxies are not included, and as a result 
the characteristic CFE at high redshifts may be lower in reality than 



Kennicutt-Schmidt law of equation )44K which were both used to for deter- 
mining the model lines. 
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Figure 9. Variation of the cluster formation efficiency (CF E) with cosmic time. Left: Predicted CFE as a function of the observed star formation rate density 
for nearby disc galaxies (diamonds, using input data from lKennicutl|[l998h . high-redshift (z = 1-2,5) normal galaxi e s (squ a res), and high-redshift (z = 1 — 
3) starburst galaxies (pl us signs). The input data (E g and f2) for redshifted galaxies are taken from iTacconi et ail (2008), Forst er Schreiber et alj 12009), 
iTacconi et all feOlOt) . and lGenzel etai] feOlOt) . To guide the eye, the curves indicate the model predictions for Q = 0.5 (sol id) and Q = 1 .5 (dashed), defining 
£2 and EgpR as in equations (5) and respectively. Right: Predicted CFE as a function of redshift and lookback time (cf. Wright 2006). Small grey symbols 
are the same as in the left-hand figure. Large solid symbols indicate the mean CFEs for nearby discs (diamond) and starbursts (triangle, using input data from 
Kennicutt 1998). The lines indicate the mean CFEs in Az = 1 redshift bins for discs (dotted) and starbursts (solid). Error bars do not represent standard errors, 
but indicate the 16 th and 84 th percentiles of the underlying distributions. 



is suggested by this sample. We conclude that while the sample 
is certainly not complete, it does highlight a change of the star- 
forming population across the covered redshift range. The trend of 
the CFE with cosmic time of Figure [9] is thus a tentative, but plau- 
sible prediction. 

Current studies of star cluster populations can reach out to 
distances of ~ 100 Mpc, corresponding to z ~ 0.024. It is 
therefore not possible to confirm the predicted evolution of the 
CFE through direct observations. Even for an indirect verifica- 
tion to be possible, it would be required that some products of 
high-redshift cluster formation have survived until the present 
day. It has long been suggested that at least some fraction of 
the old, massive globular clusters that populate the haloes of 
nearby galaxies are the surviving products of the regular high- 
redshift star formation process, and were formed through the same 
mechanisms tha t currently lead to the formation of young stel- 
lar cl usters (e.g.[Elmegreen & Efremov 1997; Kruiiss en & Cooperl 
120121) . The differences between young cluster populations and 
these old globular clusters would then be caused by differences in 



the b i rth environment (e.g . Kravtsov & Gnedinl 2005 ; El megreei] 



1201 d : IShapiro et all I2OI0F iKruiissen et al.l l2012bh and a Hubble 
time of dynamical evolution (e.g. lElmegreen & Efremov! 1 19971 ; 
iFall & Zhandl200ll ; IVesperInlll200ih . 

If some part of present-day globular clusters are indeed the 
relics of high-redshift, regular star cluster f ormationPI then we can 



20 A common argument agai nst this is the prese nce of abundance anoma- 
lies in globular clusters (e.g. iGratton et alJl2004T) , which do not appear in 
nearby, currently forming stellar clusters. However, there are no known 
resolved regions locally where the formation of such high-mass (M > 



use the framework of the present paper to address the conditions of 
their origin in some more detail. Figure [9] predicts that the CFE in 
high-redshift disc galaxies reached values almost one order of mag- 
nitude larger than in present-day discs, and some factor of 2 higher 
than in nearby starburst galaxies. This does not only imply that 
clusters forming in the early Universe were ten times more numer- 
ous at the same star formation rate. Given that young stellar clus- 
ters are formed according t o a power law initial cluste r mass func- 
tion with an index of —2 (Portegies Z wart et al]|2010l) . the size-of- 
sample effect dictates that a 1 dex increase of the CFE would also 
facilitate the formation of clusters ten times more massive at the 
same star formation rate. Considering that the star formation rate it- 
self attained values SFR > 1000 Mq yr _1 , the masses of the most 
massive clusters should have reached even higher. This statistical 2 *! 
argument suggests that the 'normal' cluster formation mechanisms 
of the nearby Universe would lead to the formation of extremely 
massive (<; 10 7 Mq) and dense stellar systems like (metal-rich) 
globular clusters in high-redshift galaxies, and that such systems 
were commonplace. 

It is possible to make a theoretical estimate of the total fraction 
of all stars in the Universe that formed in bound stellar clusters by 
integrating the CFE over the cosmic mass assembly history. Not- 



10 Mq) globular cluster equivalents can be seen in action. Therefore, it 
is unclear whether the formation of massive clusters in gas-rich galaxies 
could perhaps naturally lead to the observed abundance anomalies. 
21 A possibly more relevant (and currently unanswered) question is 
whether there are any physical mechanisms that limit the maximum clus- 
ter mass, and how these would evolve with the galactic environment (e.g. 
iGieles et all2 006a; Bastian 2008; La rsej2009h . 
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Table 4. Influence of the adopted parameters on the CFE. 



Parameter 


Minimum 


Maximum 


Impact on CFE 


(1) 


(2) 


(3) 


(4) 


<f>P 


1 


6 


+ 




1 


2 


+ 


Csn 


2Myr 


5 Myr 


+ 


t 


5 Myr 


20 Myr 




4>ib 


0.032 cm 2 s~ 3 


0.8 cm 2 s" 3 


+ 


e corc 


0.25 


0.75 




/ 


0.5 


0.9 




9 


1 


2 




4>sh 


2 


3 




VLG 

GMC 


30 M pc~ 2 


300 M Q pc~ 2 





ing that secular star formation in galaxy discs contributes a larger 
fraction of the cosmic mass assembly than starbursts for redshifts 
z<>1 (e.g. lSomerville et alj20ollKauffmann et al]2003h . the red- 
shift evolution of the CFE in Figure [9] warrants a rough, two-step 
approximation in which V — 50% for z > 0.7 and F = 10% 
for z < 0.7. As indicated above, the former of these two val- 
ues is an upper limit due to the possible incompleteness of the 
high-redshift sample. Since the (comoving) ste llar mass density in- 
creased by a factor of 1 . 6-2 since z = 0.7 (e.g. lPozzetti et alj2007l : 
iMarchesini e"tai]|2009l ; lllbert et al.ll201oh . this gives a total cosmic 
integrated CFE of r un i v ^ 30-35%. In other words, the model pre- 
dicts that up to one third of all stars in the Universe once formed in 
bound stellar clusters, while the remainder originated in unbound 
associations. 



7 DISCUSSION 

In this section, we discuss the possible sources of uncertainty in 
the presented theoretical framework, as well as future observational 
tests and potential applications in observations, theory, and numer- 
ical simulations. 



7.1 Influence of parameters and model assumptions 

The first test of the uncertainties in the presented model is to as- 
sess how the adopted parameters might influence the result. While 
the model does not rely on a very large number of parameters (see 
Table |2j, and even though we have shown that the variation of the 
CFE due to individual parameters is generally minor (see e.g. Fig- 
ure^, it is worth checking how the results change if all parameters 
would conspire in the same direction. 

The typical variations or uncertainties of the model pa- 
rameters is shown in Table [4] the fourth column of which in- 
dicates whether the CFE increases (+) or decreases (— ) with 
th e parameter in question. Th e ranges are based o n discussions 



in IKrumholz & McKeej <l2005h. iGieles et alj fe006bl), iHever et al] 
d2009|- ' 



(2009), Po rtegies Zwart et alj J2010h . iKruiissen et alj feOl lh . and 
on the appendices of this paper. Using the table, one can easily 
decide which extremes to choose in order to estimate the maxi- 
mum variation of the CFE. Additional variation may come from 
the adopted star formation law (see Figures|2]and[8]( and the inclu- 
sion of radiative feedback (see Appendix ICl. We therefore consider 
two cases with the following changes to the fiducial model. 
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Figure 10. Worst-case-scenario uncertainty on the predicted cluster forma- 
tion efficiency for the fiducial model (solid line). The dashed lines represent 
the input physics and parameters from Table|4] which were chosen to min- 
imize (lower dashed) or maximize (upper dashed) the CFE. Note that these 
parameter sets are highly improbable (<; 2.9<r), implying that the actual la 
uncertainty is comparable to the dotted lines (see text). 



feedback from Appendix [C] and adopting the minimum and maxi- 
mum parameters of Table|4]leading to the lowest CFEs. 

(ii) The high-CFE case, using the IKrumholz & McKed ( 120051) 
star formation law of equation dl2t and adopting the minimum and 
maximum parameters of Table|4]leading to the highest CFEs. 

The resulting extremes of the CFE are shown as a function of the 
gas surface density in Figure QJj] together with the predicted CFE 
for the fiducial parameter set. The 'uncertainty' around the fiducial 
T-Eg relation is of the same order as the scatter of the observa- 
tions around the 'typical' model in Figure [6] While in principle 
this implies that differences between theory and observations could 
be accounted for by model uncertainties, we remind the reader 
that the discrepancies are smaller when choosing appropriate in- 
put variables for each individual galaxy (Figure[7}. This should be 
expected, because the two extreme cases of Figure [10] are highly 
improbable. For illustration, if one interprets the ranges of Table [4] 
conservatively as lcr limits, the chance that all ten parameters con- 
spire to give the extremes of Figure [TTJl is 0.16 10 ~ 10 -8 . Of 
course, this assumes that all parameters are independent and influ- 
ence the CFE to an equal extent. In practice, the uncertainty on the 
CFE is mainly caused by three independent variables, being 4>p, 
t sn and £ C orc. But even in this case, the probability of the extremes 
is 0.16 3 ~ 0.004, which is a 2.9a limit. Assuming a log-normal 
distribution of errors, the actual la logarithmic uncertainty on the 
model predictions (dotted lines in Figure UOt should therefore be 
about 0.35 times the logarithm of the extreme case. This uncer- 
tainty decreases with the gas surface density, from or ~ 0.3 dex at 
E g ~ 2 M pc -2 to or < 0.1 dex at E g ~ 10 3 M pc -2 . 

In the discussion so far, we solely focused on the uncertainty 
of the adopted model parameters. However, other assumptions that 
are physical in nature might also influence the results of this paper. 



(i) The low-CFE case, including the prescription for radiative 



(i) Magnetic fields are not included explicitly, i.e. there are no 
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physical terms for magnetic fields in any of the equations. This 
means that the influence of magnetic fields on cluster formation 
cannot be studied using the theoretical framework of this paper. 
However, their influence is implicitly present to some degree in 
the assumption of an ISM overdensity PDF that is consistent with 
weak-field magnetohydrodynamics simulations, and likewise in the 
adopted specific star formation rate per free-fall time. The model 
predictions should thus not be strongly affected by any (future) ex- 
plicit formulation. A possible, first-order modification of the the- 
ory to allow for strong magnetic fields would be to add a fac- 
tor /3o/(/3o + 1) to the second term of equation (|4j, where /3o = 
fth/ fmag is the ratio of the thermal pressure to th e magnetic pres- 
sure jPadoan & Nordlundll201lMMolina et alj2012h . However, this 
does not account for the influence of a strong magnetic field on the 
physics of gravitational collapse. 

(ii) When considering individual galaxies, the overdensity PDF 
of the ISM might not follow the log-normal shape of equation 

It is known from nume rical work that deviations arise mainly at low 
overdensities (e.g. IWada & Norman] 120071 : iTasker & Tanfl2009l) . 
which is in the range where the formation of stars and stellar clus- 
ters proceeds slowly. As a result, the impact on the absolute star 
formation rate and CFE are generally well within the other uncer- 
tainties. 

(iii) The model assumes that star-forming regions reside in gas 
discs that are in hydrostatic equilibrium. This could pose a problem 
when modeling star formation in galaxy mergers or in dynamic, 
high-redshift environments that might be out of equilibrium. Ac- 
knowledging this, it is important to note that when gas cools to 
form stars, the energy d issipation drives the formation of a disc 
(e.g. iHopkins et al. 2009). This means that while the presence of a 
disc might not be obvious on a global scale, the main loci of star 
formation should follow a disc-like morphology. This means that 
our assumption is satisfied at least locally. 

(iv) In our model, star formation in intermediate-density regions 
is halted by the large-scale obstruction of further gas inflow by 
supernova feedback. The role of different feedback mechanisms 
in truncating the star formation process is traditionally a very ac- 
tive topic in the literature, and hence alternatives such as stellar 
winds and radiative feedback have been discussed extensively (e.g. 
iMurrav et alJl201oT) . It is important to note that the total energy 
inputs from these feedback mechanisms are in the same ballpark, 
which implies that the influence of this approach on the CFE is not 
major (see Appendix|B]and|C]l. The description of feedback mainly 
serves the purpose of truncating star formation on a time-scale that 
is broadly consisten t with observations. With tf b ~ 3 Myr, this is 
indeed the case (e.g. lPortegies Zwart et alj2010l) . 

(v) The cruel cradle effect is modelled by calculating a critical 
overdensity below which tidal perturbations are able to inhibit the 
collapse of a star-forming region. In doing so, we neglect the den- 
sity increase of such regions during collapse and the correspond- 
ingly increased resistance to tidal disruption. This is especially rel- 
evant in view of the stochasticity of tidal perturbations, which are 
unlikely to start precisely at t = and hence do not act on the initial 
ISM density of a star-forming region, but on a (somewhat) higher 
density. The quoted impact of the cruel cradle effect is therefore an 
upper limit. If we consider the extreme case of only a single per- 
turbation per free-fall time(3 then on average regions have the po- 
tential to collapse halfway before experiencing their first tidal per- 



Even lower encounter rates would imply the cruel cradle effect is negli- 
gible due to the adiabatic dampening of equation (38). 



turbation, causing at most a factor-of-eight density increase. This 
is similar in magnitude to the variation of the product fgcj> a h in 
Table[4] which contributes only a small fraction of the total uncer- 
tainty in Figure [10] We conclude that the treatment of tidal pertur- 
bations is adequate for the problem at hand. 

7.2 Observational tests and applications 

Because the presented theory makes a number of clear assumptions 
and predictions, there are several possible ways in which (compo- 
nents of) it can be tested. The model can also be applied as a tool in 
interpreting observational data. Note that Fortran and IDL routines 
for calculating the CFE with our model are publicly available (see 
Appendix lAl>, 



7.2.1 Observational verification with Gaia using kinematics 

On small spatial scales, the model postulates that the formation of 
bound stellar clusters proceeds hierarchically, reaching high SFEs 
in high-density regions. In this picture, the gas and stars decouple 
dynamically as the stellar component accretes and shrinks, allow- 
ing it to attain virial eq uilibrium while embedded in a cocoon of 
gas dKruiissen et all2012al) . The direct implication is that the clas- 
sical picture of early cluster disruption by infant mortality does not 
work as such. Gas expulsion does not unbind the dense, spherically 
symmetric systems that went through hierarchical merging and vi- 
olent relaxation. Such systems should largely survive gas expul- 
sion. Instead, gas removal only affects star-forming regions on the 
scales where they are still gas-rich and substructure is still present. 
This prevents the further merging of individual sub- or protoclus- 
ters. With kinematic data from Gaia, it is possible to derive from 
dispersed structure whether gas expulsion acted as classical infant 
mortality on a spherically symmetric cluster, or whether the stars 
still followed a hierarchical structure at the onset of the dispersal. 
In the former case, the velocity vectors would be radial and have 
a common origin, whereas in the latter case, the velocity vectors 
should be traceable to multiple centres and filamentary structure. 

Gaia will also be able to measure the relative occurrence of 
the dispersion of stellar structure due to being naturally unbound, 
and dispersion due to the cruel cradle effect. Tidal perturbations 
induce structure in velocity space, in that the stellar velocity vectors 
tend towards the plane of the interaction lKrui isser]|201ll) . Such 
kinematic evidence would be visible until the stars reach the tidal 
boundary of the region, and allows one to quantify the fraction of 
early cluster disruption due to the cruel cradle effect. In the Milky 
Way, we estimate that while about half of all stellar structure is 
affected by tidal perturbations, a large fraction of this is already 
naturally unbound. Therefore, the cruel cradle effect is critical only 
for S; 1% of all young stellar structure, and unbinds some 10% of 
all naturally bound structure (see Table[3}. 

Our model can also be considered using the relation between 
the local SFE e (on any scale) and the fraction of stellar mass that 
will remain bound upon gas removal 7 (see Appendix [Dj». The 
model predicts a certain distribution of SFEs, which through the 
adopted 7(e) relation yields a distribution of bound mass fractions. 
With Gaia, it will be possible to identify expanding haloes of es- 
caping stars around young, unembedded clusters and associations. 
Correlating the masses of these haloes with the cluster masses will 
provide a direct measure of 7. The measured distribution of bound 
fractions should be similar to the predicted distribution of SFEs that 
follows from the properties of the ISM. Another prediction of the 
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model is that expanding stellar haloes should be sparse or absent 
around massive and dense stellar clusters! 23 ! 



7.2.2 Observing gas-embedded regions with ALMA 

The theory of this paper causally connects the properties of the ISM 
to the eventual population of bound stellar clusters. This is well- 
suited for verification with the sub-millimeter capacity of ALMA. 
Within the Milky Way, ALMA observations can trace the dynami- 
cal decoupling of gas and stars in protoclusters. It should also en- 
able a statistical analysis of the protocluster population, and estab- 
lish how the density PDF of bound stellar structure follows from 
the density PDF of the ISM. Specifically, the model predicts that 
the high-density tail of the density spectrum should promote bound 
cluster formation (see Figure QJ. With ALMA observations it is 
possible to witness and quantify this process as it takes place. 

In extragalactic systems, the study of populations of proto- 
clusters and their dynamical states will provide a direct measure 
of the dependence of cluster formation on the galactic environ- 
ment. Our theory predicts that the key factor that sets the CFE 
is the gas surface density, with angular velocity and disc stabil- 
ity being the secondary variables. Are gravitationally bound proto- 
clusters more numerous in high-density galaxies than in gas-poor 
galaxies? An obvious nearby example is M31, where we predict 
that only 2% of all star formation eventually produces bound stel- 
lar clusters. Contrasting a sub-millimeter census of the protocluster 
population in this galaxy with a similar survey of systems like M83, 
M82 or NGC 3256 would be very valuable for our understanding 
of how the galactic environment influences the conversion of gas 
into bound stellar clusters. 



7.2.3 Application to observed cluster populations 

Another evident avenue for the application of the presented theory 
is in the study of observed cluster populations. Firstly, our frame- 
work provides a theoretical benchmark to test observational esti- 
mates of the CFE against. The model requires input variables that 
are easily obtained for large samples of galaxie^fj and hence any 
new observed value of the CFE can be compared to the model 
prediction (see Table [3). Secondly, this approach can also be re- 
versed. When the star formation history and cluster age distribu- 
tion are known, having an estimate for the CFE will enable one to 
con strain the cluster disruption law in that particular galaxy (see 
e.g. ISilva- Villa & Larsenll201lh . Thirdly, it is straightforward to 
convert observed cluster age distributions (e.g. iGieles et alj|2005l ; 
ISmith et al. 2007) to star formation rates using the modelled CFE. 
A case in point is the recent study by Ford & Bre gmanl J2012h . who 
use the cluster age distribution to derive extremely low star for- 
mation rate densities Esfr ~ 10~ 5 M© yr _1 kpc -2 in nearby 
elliptical galaxies, but do so assuming that all stars form in clus- 
ters. The model relation of Figure [6] shows that at such low sur- 
face densities only a minor fraction of all stars forms in clusters. 

23 Low-mass, dense clusters could still have a substantial halo bec ause of 
rapid evaporation due to two-body relaxation (Moeckel et al. 20lJ). Also 
note that due to the hierarchical nature of star and cluster formation, it is 
possible that a non-negligible m ass fraction may form dispersed in the pe- 
riphery of a compact cluster (e.g. lBressert et alj2012ah 

24 As mentioned throughout the paper, only the gas surface density is es- 
sential for calculating the CFE. The angular frequency and Toomre Q pa- 
rameter can be approximated if no observational input is available. The de- 
fault f2 follows from equation j9), and Q ~ 1.5 for regular disc galaxies. 



A self-consistent solution is obtained by taking F ~ 0.01, which 
increases the estimated star formation rate density to Esfr ~ 
10 -3 Mq yr _1 kpc -2 . This example illustrates that only when 
used together, the CFE and the cluster age distribution provide an 
accurate, independent tracer of star formation activity. 

The CFE is also a potentially useful quantity for globular clus- 
ter studies. We have shown in Sj6]and Figure [9] that the CFE was 
much higher (r ~ 50%) during the formation epoch of metal- 
rich globular clusters than it is in the local Universe (V ~ 10%). 
When studying globular cluster populations it is relevant to know 
that these objects represent the survivors of a cluster population 
that initially constituted half of all star formation, rather than just a 
minority. Combining this with prescriptions for cluster disruption, 
it is possible to derive the coeval stellar mass of present-day globu- 
lar cluster populations, which refines their use as tracers of galaxy 
assembly. 



7.3 Future applications in theoretical and numerical work 

The model is applicable to a range of theoretical and numerical 
problems. Fortran and IDL routines for calculating the CFE with 
our model are publicly available (see Appendix [At. 



7.3. 1 Models of star and planet formation 

There is some steadily a ccumulating evidence that the stellar initial 
mass function (IMF, e.g. lKroupall200ll ; IChabriej|2003l) is universal 
in the local Universe, and d oes not vary between clustered and dis- 
persed star-forming regions ide Wit et al. 2005; Parker & Goodwin 
l2007l : lBressert et alj|2 012a). However, indirect IMF measurements 
in giant elliptical galaxies have suggested t hat these systems could 
have an excess of low-mas s stars (e.g. [van Dokku m & Conrovl 
|2010| : ICappellari et al.l |2012|). Since giant el l ipticals assembled 
at very high redshift (e.g. Icimatti et alj |2004 iDaddi et all 120051 ; 
iBundv et alj|2006l ; iNaab et alj|2007h and at correspondingly high 
surface densities, the CFE must have been high as well. It is un- 
clear whether a possible variation of the IMF is caused by the vari- 
ation of the large-scale galactic environment, or by an influence of 
high-density clustered environments on the star formation process. 
The absenc e of IMF variations in Galactic stellar clusters favours 
the former (Bas tian et al]|2010h , but the possibility that the high- 
redshift conditions of cluster formation may have affected star for- 
mation cannot be ruled out. In such a case, the CFE would deter- 
mine to what extent the IMF of an entire galaxy changes. 

On the planet formation side, the influence of clus- 
tered star formation is more evident. It has been known for 
some time that circumstellar discs an d planet formation are 
affected by exte rnal photoevaporation ( Scally & Clarke 2001; 
lAdams et alj 120041) and dynamical encounters with other stars 
dOlczak et al. I 120061; iDukes & Kramhoizl 1201 2t IParker & Ouanzl 
|2012| ; dc Juan Ovclar et al .120121) . Because both effects occur on an 
encounter time-scale t cnc oc p _1 , external perturbations of plane- 
tary systems or circumstellar discs are more important in bound 
clusters, and hence in galaxies with high CFEs. Recent observa- 
tional work suggests that environmental effects can truncate cir- 
cumstellar discs above stellar number densities S > 10 3 ' 5 pc -2 
on a 1-Myr time-scale Jde Juan Ovelar et al.ll2012T) . A natural ex- 
tension of this argument would be that in galaxies with gas surface 
densities E g ~ 10 3 Mq pc~ 2 only few planets may form in rela- 
tive isolation because T ~ 70%. When modeling the properties of 
the integrated planet population of such a galaxy or galactic region, 
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it may be necessary to account for the fraction of stars (and planets) 
that was bom in bound stellar clusters. 



7.3.2 The chemical composition of stellar clusters 

A range of recent papers focus on the chemical (in)homogeneity 
of globular and open clusters in th e Milky Way (e.; 



way 

Gratton et al] |2004 l_Randich et alj 120061 : |Pe Silva et ail 120061 : 
Bland-Hawthorn et al .l201Ch . While homogen eity may be desirable 



to identify cluster remains in stellar streams dPe Silva et alj|2007l ; 
IWvlie-de Boer et al.ll201(il; iBland-Hawthorn et alj |2010L although 



Maiewski et al.l2012l for a counterexample), there appears to be 



a natural transition to chemically inhomogene ous clusters at higher 
masses, as evidenced by old globular clusters jGratton et all2004l) . 
This trend seems to be qualitatively consistent with the picture pre- 
sented in our model. For a given density and free-fall time shorter 
than the star formation time-scale, the spatial extent from which 
mass is hierarchically accumulated increases with GMC mass, thus 
allowing for a larger degree of inhomogeneity in the eventual clus- 
ter. Alternatively, it has been argued th at the collapse time-scale in - 
creases weakly with the GMC mass (Bland-Hawthorn et al.ll2O10h . 
This would enable self-enrichment and increasing chemical inho- 
mogeneity at higher cluster masses, which is the prevailing expla- 
nation for the chem ical inhomogeneities in globular clusters (e.g. 
lGrattonetal1l2004l) . However, we should remind the reader that 
this relies on the identification of mass quanta (cf. lTan et all2006h . 
whereas there is no intrinsic mass-scale present in the continuous 
formulation of this paper. Quantized clusters only appear after the 
star formation process has ended, the time of which depends princi- 
pally on the density. As a result, a mass-dependence of the collapse 
time-scale can only be obtained from our model using additional 
assumptions about the relation between mass and density or radius. 
The above considerations imply that high-mass, low-density col- 
lapse should lead to the largest chemical inhomogeneities, provided 
that such inhomogeneity increases with both the spatial extent of 
the region and the star formation time-scale. 



7.3.3 Numerical simulations of star cluster populations 

Analytic models for the evolution of star cluster populations have 
recently been impleme nted in numerical simulations of galaxy for- 
mation and evolution ([Prieto & G nedin 2008; Muratov & Gnedinl 
l20ld :l Kruiissen et alj|201ll . l2012bT) . This is a promising new av- 
enue for investigating the origin of old globular cluster populations, 
and for relating the properties of these populations to the high- 
redshift galactic environments they formed in. So far, these simula- 
tions have only included a physical treatment of the dynamical evo- 
lution and disruption of stellar clusters. For cluster formation, these 
studies have needed to assume that all stars (or a certain arbitrary 
fraction thereof) formed in bound stellar clusters. This assump- 
tion is especially inconvenient because the properties of the surviv- 
ing cluster population are strongly influenc ed by how the cluster 
disruption rate vari es with space and time (Kruiisse n et alJuOllI ; 
Bastia net alJl2012T) . which could be offset (or enhanced) by sim- 
ilar variations of the CFE. Numerical simulations that follow the 
star cluster population with a subgrid, semi-analytic model should 
thus benefit from the theoretical framework that is presented in this 
paper. 

Contrary to observed galaxies, the state of the ISM in numeri- 
cal simulations is ill-described by quantities such as the gas surface 
density, angular velocity and Toomre Q parameter, which all rely 



on a suitable projection or choice of coordinate system. Instead, for 
implementation in numerical simulations it is more useful to con- 
sider strictly local quantities. The best way to do this is to discard 
the assumption of an equilibrium disc (see §2.4i , and to determine 
a subgrid overdensity PDF of the ISM from the local gas density 
pioc, gas velocity dispersion <ji oc , and sound speed of the cold ISM 
c s ,ioc- The local Mach number then follows as Alloc = <xi oc /cs,ioc, 
which means that the PDF of the overdensity a^i oc = p g /pi oc is 
uniquely determined (see §2.3\ . The conversion to a PDF of the ab- 
solute density p s then proceeds simply through multiplication with 
pioc which is determined by integrating over the smoothing kernel 
(for Smoothed Particle Hydrodynamics or SPH codes) or averaging 
over neighbouring cells (for grid codes). Most other steps in our 
model are formulated locally, except for the feedback-limited SFE 
£fb and the critical overdensity for surviving the cruel cradle effect 
x cca . Both can be rewritten to be a function of the local variables 
introduced here. For the SFE we write 



£fb,loc = 



sSFR ff t s 



1 + A 1 + 



^fbSSFRfftln^lo 



(46) 



(47) 



The critical overdensity becomes 

_ 87 '.5^/ttG f 'g(fi Bh T,GMct , , . 

^ccejoc — (pad (^Ccce,loc / 5 

O"loc 

where the GMC surface density can either be taken to be constant 
or be written as Eqmc = max {Sgmo ^g} as m me rest °f tn i s 
paper. In that case, E g can be estimated very roughly by again as- 
suming an equilibrium disc and writing 



2 Pi°cof oc 



, . (48) 

nO(pp 

Note that the aforemention ed models of IPrieto & Gnedinl 

(2008) and Kruii ssen et al.1 d201ll) account for the tidal disruption 
of clusters from the moment they are born. In such a context, the 
cruel cradle effect implies that clusters have high chances of being 
disrupted at young ages, due to still residing in their dense natal en- 
vironment. This is complementary to the manifestation of the cruel 
cradle effect that is treated in this paper, in which it prevents the 
global collapse of star-forming regions to form bound clusters. The 
best ways of connecting both forms in numerical work are either 
(1) to delay the tidal disruption of clusters until the time t at which 
the CFE is evaluated, or (2) to exclude the cruel cradle effect in 
the calculation of the CFE (i.e. r/100% = /bound) without delay- 
ing the tidal disruption of clusters. Because lower-density regions 
are more easily disrupted, the first option underestimates the CFE 
by neglecting the density increase of collapsing regions (although 
this effect is minor, see *j7jj. Conversely, the second option may 
overestimate the CFE if the modelled evolution of the cluster radius 
does not include the collapse of clusters during their formation, and 
their correspondingly lower densities at t — 0. 

The local and global prescriptions for the CFE are both in- 
cluded in publicly available Fortran and IDL routines (see Ap- 
pendix[A]l- This enables a straightforward implementation into ex- 
isting numerical models. 



7.3.4 The feedback efficiency in cosmological simulations 

In numerical simulations of hierarchical galaxy formation and evo- 
lution, the injection of feedback energy into the ISM generally 
includes an efficiency factor Xfb < 1 that indicates the frac- 
tion of the energy that successfully couples to the gas. This feed- 
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back efficiency is often taken to be Xfb ~ 0.05 ( see the dis 



cussion in Appendix iBl and, e.g., Efstathioul 2000l: lAbadi et"al1 



2003; Rob ertson et"alll2005l : IBubois & Tevssierll2008h . but unfor- 
tunately, a detailed understanding of this parameter is still lacking 
dSilkil 19971; iMac Low & Ferraralll999l ; lNavarro & Steinmetj|2000l ; 
iDib et alj |2006). It is not hard to imagine that a systematic varia- 
tion of the feedback efficiency during certain episodes of galaxy 
evolution c ould affect the pr operties of the resulting galaxy popu- 
lation (e.g. Sales et alj |2010). and therefore it is relevant to inves- 
tigate how small-scale physics may influence large-scale galactic 
outflows. 

If feedback sources are clustered, then some fraction of 
the feedback energy is deposited within existing bubbles, which 
can allow the energy to accumulate at a rate higher than the 
cooling rate, and hence increases the feedb ack efficiency (e.g. 
IStrickland & Stevens! 19991 ; I Krause et alj|2012h . Quantitatively, the 
relative change of the feedback efficiency for clustered (xfb.ci) as 
opposed to dispersed (xfb.o) feedback can be parameterized as 
4>d = Yfb.c i /Xfb.o - According to the recent numerical work of 
iKrause et alj d2012f) <j> c i ~ 2, indicating that feedback from clus- 
tered star formation is more effective than if it were dispersed. 
Combining this with the model of the present paper, it is possible to 
formulate a composite feedback efficiency from the clustered and 
dispersed parts of star formation. When simulating the assembly 
of galaxies and their evolution, the CFE can easily be calculated 
for each star-forming particle (in SPH) or cell (in a grid code) as 
detailed in i)7.3.3l The composite feedback efficiency Xfb then fol- 
lows to first order by noting that a fraction T of the particle (or cell) 
mass is clustered and a fraction (1 — T) is dispersed: 



Xtb = (1 - T)Xfb,o + I>clXfb,0 

= [i + (<t>d - i)r] xfb.o. 



(49) 



The maximum variation of Xfb is a factor of cj> c i per definition. 
While for cj> c i ~ 2 this is not substantial, the systematic rise of the 
CFE at high densities (and star formation rates) could imply that 
the galaxy population is affected over the course of a Hubble time. 
With equation J49l > it will be straightforward to investigate thisF^l 



7.3.5 Semi-analytic galaxy modeling 

Galaxy populations are often modelled semi-analytically to alle 
viate the computational ef fort required for the numerical Simula 
tion of their formation (e.g. White & Frenk 1991; Kauffman n et al 



1993; Somerville & Primack 1999; Coleetal. 2000; Bower et al 



20od ; lGuo et al.l201 1). These analyses are covering the c haracteris 



tics o f the modelled galaxies in increasing detail (see e.g. lGuo et al.l 
l20Tlh . It seems only natural that the origin and evolution of the 
star cluster populations hosted by the modelled galaxies will be in- 
cluded in future work. This will enable an understanding of how 
the co-formation and co-evolution of (old) star cluster populations 
and galaxies may eventually produce the present-day cluster sys- 
tems across the en tire galaxy mass range (for typical observables 
of interest, see e.g.lPeng et al.1120061 ; I Jordan et al] 2007; Pen g et al] 
120081" : iBurkert & Tremainell20 1 d) . 

A self-consistent analytic model for cluster populations that 
can be implemented in semi-analytic models of galaxy formation 
requires (1) a prescription for the formation of cluster populations 



25 Note that this prescription only applies to numerical resolutions at which 
the clustering of star formation is not resolved. At higher resolution, the 
influence of clustered star formation may begin to emerge naturally. 



and (2) a prescription for their dynamical evolution. Both compo- 
nents of the problem are still unsolved. While the theory of the 
CFE presented in this paper provides the fraction of star formation 
occurring in bound stellar clusters, the initial characteristics of the 
individual clusters within the population are yet to be understood in 
more detail. Their mass spec trum can be taken from observations 
jPortegies Zwart et al] l20l6T) and has been shown to arise from 
the hierarchical nature of the ISM l lElmegreen & Falgaronel! 1996), 
but the understanding of the corresp onding upper (e.g. Gieles et al. 
| 2006a ; lBastianll2008l; lLarserJl2009r) and lower (e.g. lLada & Ladal 
l2003l ; lMoeckel et alj2012T> cluster mass limits requi res more work , 
as does the distribution of initial cluster radii (e.g. lLarsenll2004l) . 
These quantities all need to be described with reasonable accu- 
racy to enable the inclusion of cluster formation and disruption in 
a single semi-analytic framework. The evolution of cluster popu- 
lations due to dynamics should be addressed by building on more 
detailed numerical modeling, both accounting fo r the internal clus- 
ter dynamics (e.g . iBaumgardt & Makinol 120031 ; iKruiissetj 120091 ; 
I Gieles et al] 201 lh and the disruptive influence of th e galactic en- 
vironment " dGieles et al j|2006bl ; lKruiissen et al j|201 ll) . While some 
of the above ingredients can only be roughly estimated at present, 
the recent progress of the field is encouraging. By addressing the 
above questions, a first inclusion of cluster populations in semi- 
analytic galaxy population modeling should become possible in the 
next couple of years. 



8 CONCLUSIONS 

We have described and applied an analytic theory of the cluster 
formation efficiency (CFE), i.e. the fraction of star formation oc- 
curring in bound stellar clusters. The conclusions of the paper are 
as follows. 

(i) The presented theory of cluster formation provides a frame- 
work in which the galaxy-scale environment is related to the prop- 
erties of star-forming regions. Bound stellar clusters naturally arise 
from the high-density end of the density spectrum of the interstel- 
lar medium (ISM). Due to short free-fall times, these high-density 
regions can achieve high star formation efficiencies (SFEs). This 
makes them insensitive to gas expulsion and enables them to form 
bound clusters. In regions of lower density, the SFE is lower and 
the collapse into spherically symmetric structure is not completed. 
As a result, such regions (or associations) remain gas-rich and are 
unbound upon gas expulsion. ( Sj2.3442.6t 

(ii) In the picture described here, gas expulsion does not af- 
fect centrally concentrated stellar clusters like in the classical sce- 
nario of 'infant mortality', but inhibits the merging of hierarchically 
structured stellar groups and associations. ( i)2.6| and Appendix IDT) 

(iii) Additionally, we have included the 'cruel cradle effect', 
which represents the tidal perturbation of star-forming regions by 
surrounding GMCs. This is shown to give a second-order decrease 
of the CFE in normal disc galaxies, and is most important in high- 
density galaxies. ( H2Ji 

(iv) The CFE is obtained by integrating the probability distribu- 
tion function (PDF) of the densities at which star formation pro- 
ceeds, and the corresponding part thereof that is both bound locally 
and sufficiently dense to survive the cruel cradle effect. (i )2.8| and 
Figure [TJ 

(v) Our theoretical framework allows the prediction of the CFE 
as a function of galaxy-scale observables: the gas surface density 
Eg, the angular velocity fl, and the Toomre Q parameter. For use 
in numerical simulations, a local formulation is also given, which 
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uses the gas volume density p\ oc , the gas velocity dispersion <7i oc , 
and the sounds speed c Si i oc . ( i]2.8l and i]7.3,3l l 

(vi) The model predicts that only some fraction (1-70%) of all 
star formation results in bound stellar clusters. This CFE increases 
with the gas surface density of the galaxy, because in denser galax- 
ies a larger fraction of the ISM is pushed into the high-density 
regime where bound cluster formation occurs. However, at surface 
densities E g <; 10 3 M© pc -2 the CFE is limited by the cruel 
cradle effect. For E g ^ 10 25 M pc~ 2 , the CFE deer eases with 
increasing angular velocity f2, while at lower densities it is gener- 
ally insensitive to Q. The CFE increases with Toomre Q, and hence 
with the stability of the galaxy disc. (S)3T] fl3~2l and COT l 

(vii) We find remarkable agreement between our theoretical 
model and the observed relation between the CFE and the star for- 
mation rate density Esfr in dwarf, spiral and starburst galaxies. 
Also when modeling individual galaxies, the model predictions are 
in excellent agreement with the observed CFEs. (§4$ 

(viii) The model is applied to investigate the spatial variation 
of the CFE within galaxies, using parameter sets that resemble the 
Milky Way, M3 1 , and M5 1 . At the low densities of the Milky Way 
and M31, the CFE simply follows the relatively flat surface density 
profile of the gas. However, for the exponential gas density profile 
of M51, the CFE increases towards the centre of the galaxy, but (at 
high densities) does so less than the gas density profile. (Sj5) 

(ix) We use samples of nearby and high-redshift galaxy discs 
and starbursts to predict the evolution of the CFE across a Hubble 
time of star formation. For galaxy discs, the CFE is found to de- 
crease from T = 30-70% to T = 5-10% between redshifts z ~ 2 
and z = 0, whereas starburst galaxies exhibit a moderate decrease 
of r = 20-60% to T = 10-50% over the same redshift range. 
High-redshift disc galaxies may thus have been the most efficient 
cluster factories in the history of the Universe, although presently 
the CFE is generally higher in starburst galaxies than in discs. We 
estimate that up to r un i v = 30-35% of all stars in the Universe 
was formed in bound stellar clusters. ([j6} 

(x) The predictions of our theoretical framework are well-suited 
for verification with Gaia or ALMA, and also for application to 
observational studies of cluster populations. We provide specific 
recipes for use of the model in theoretical and numerical work, 
such as numerical simulations of galaxies and their star cluster pop- 
ulations, the feedback efficiency in cosmological simulations, and 
semi-analytic galaxy models. ( i]7.2| and £|7.3b 

We have presented a new model that describes how populations of 
bound stellar clusters form, and that predicts which fraction of all 
star formation is constituted by these clusters (the cluster formation 
efficiency or CFE). The model enables the calculation of the CFE 
over a range of galactic environments, and shows that star clusters 
are no fundamental unit of star formation, but instead are a possible 
outcome. 
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APPENDIX A: ROUTINES FOR CALCULATING THE 
CLUSTER FORMATION EFFICIENCY IN FORTRAN OR 
IDE 

Fortran and IDL routines for calculating the CFE 
with the presented model are publicly available at 
http://www.mpa-garching.mpg.de/cfe For both languages, 
two versions of the model are provided. The first version of the 
routines adopts the global formulation that is followed throughout 
this paper, with as req uired input va riables the gas surface density 
of the galaxy E g , the l lToomrell 19641) Q parameter, and the angular 
velocity fi. Note that Q and Q, are optional. If they are left 
unspecified, the routines set Q — 1.5 and define Q, according 
to equation {9]l. Additional optional input parameters allow the 
user to adjust the star formation law or several of the parameters 
from Table [2] The second version of the routines adopts the local 
formulation from i)7.3.3l with as required input variables the local 
density pi oc , the local gas velocity dispersion <ti oc , and the sound 
speed of the cold ISM c Si i oc . 



APPENDIX B: THE FEEDBACK TIME-SCALE AND A 
GRADUAL TRUNCATION OF FEEDBACK 

This appendix considers the scenario in which the inflow of gas into 
a star-forming region is halted by supernova feedback. The energy 
injection rate into the cold ISM per unit stellar mass (f>ib due to 
supernovae can be estimated by writing it as 



At 



(Bl) 



where Xfb is the fraction of the supernova energy that couples 
to the cold ISM, <f> an is the number of supernovae per unit stel- 
lar mass, E Bn — 10 51 erg is the typical energy per supernova, 
and At = t 1 ^ — t Bn represents the time interval during which 
supernovae occur, with 4„ st the time at which the last super- 
nova explodes, and t sn the time of the first supernova as in H2.5\ 
This assumes that the supernova rate is approximately constant 
during At. The efficiency of the kinetic coupling of the super- 
nova energy to the ISM Vfb has been the subject of a wide range 
of studies (e.g. ISilkll 1 9971; lEfstathioul 120001: iNavarro & Stei nmetz 
200(1 lAbadi et alj|2003l ; iRobertson et alj|2005l : iDib et alj|2006l) ."lt 
is smaller than unity due to (1) rapid cooling in the cold ISM and 
(2) the porosity of the hierarchically structured ISM. In the liter- 
ature, values in the range Xfb = 0.01-0.25 have been reported, 
which illustrates that this quantity is the most uncertain of all vari- 
ables in equation l lB 1 1 . The remainder of this appendix assumes 
Xfb = 0.05, which is the logarithmic mean of the quoted range 
and is a commonly used value in numerical simulations. Physi- 
cally, this corresponds to a scenario in which about 10% of the su- 

[ ernova energy goes i nto the kinetic energy of an expanding shell 
Thornton et al.l 19981) and the total surface angle of the porous, 
cold ISM takes up about 50% of the spherical surface around the 
star-forming region. 

The number of supernovae per unit stellar mass is obtained by 
assuming a stellar initial mass fun ction (IMF) and integrating over 
the appropriate mass interval. For a l S alpeterl d 1 95 5h IMF dn /dm oc 
m~ 2 ' 35 in the mass range 0.1-100 M© and a minimum mass of 
8 Mq for a star to result in a supernova, this implies 
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Figure Bl. Left: Influence of <f>fo on the global cluster formation efficiency (CFE) as a function of the gas surface density of the disc E g , assuming Q = 1.5 
and writing f2 as a function of E g according to equation (9). The dotted, solid, and dashed lines show the relation for </>fb = {0.016, 0.16, 1.6} cm 2 s — 3 , 
respectively, hardly affecting the CFE despite varying by two orders of magnitude. Right: Same as the left-hand panel, but showing the influence on the CFE 
of a gradual truncation of star formation by feedback. The solid line shows the abrupt truncation as in £|2.5I while the dashed line represents a linear decrease 
of the specific star formation rate per free-fall time, obtained by substituting t cq — > t cq /2 in equation )2H . The CFEs are almost indistinguishable. 
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(B2) 



The time interval At during which the supernovae occur is set by 
t sn = 3 Myr (see f|23) and ffi* = 40 Myr (fo r the Padova stellar 
evolutionary isochrones, see Mar igo et alj2008l) . This range is very 
insensitive to the metallicity, as i^f* only varies by 0.05 dex in the 
metallicity range z = 0.0004-0.03. 

Substituting all variables into equation l lB 1 1 now provides 
0fb = 0.16 cm 2 s -3 = 3.2 x 10 32 erg s _1 M Q _1 , with a total 
uncertainty of ~ 1.5 dex. While this might sound substantial, its 
effect on the CFE is almost negligible, as is shown in the left-hand 
panel of Figure IbTI The curves represent the CFE as a function of 
gas surface density for different 0fb spanning two orders of mag- 
nitude, yet the variation of the CFE is always smaller than 0. 1 dex. 
The reason is that 0fb influences the star formation efficiency at 
low overdensities. It mainly controls the critical overdensity where 
the feedback-induced truncation of star formation becomes so inef- 
ficient that it is still ongoing at the time of observation. This critical 
point happens to lie in the overdensity range where the PDF is quite 
flat and where no bound stellar clusters are formed due to the cruel 
cradle effect. The overdensities where bound clusters do form is 
characterized by such powerful feedback that typcally i cq <C t sn 
regardless of the value of 4>{b, or by such short free-fall times that 
star formation is near-optimal, with e w e corc . In other words, the 
uncertainty on (f>fb may translate into a substantial shift of the criti- 
cal overdensity that separates ei nc from en, in equation J24b . but this 
hardly influences the integration of equation J26t and thus leaves 
the CFE unaffected. 

So far, the truncation of star formation by supernova feedback 
has been considered to be instantaneous once pressure equilibrium 
between the feedback energy and the cold ISM has been reached. 
As a first-order correction for this assumption, it is possible to adopt 
a linear decrease of the specific star formation rate per free-fall time 
between t sn and tft,. This is mathematically equivalent to writing 
t C q — > tcq/2 in equation ( 1211 ). Clearly, this reflects a much smaller 



change than the above variation of and indeed the right-hand 
panel of Figure IbTI shows that the CFE for a gradually decreasing 
star formation rate is nearly identical to the case of instantaneous 
truncation. 



APPENDIX C: THE TRUNCATION OF STAR 
FORMATION BY RADIATION PRESSURE 



In §2.5\ the local SFE is determined by allowing star formation to 
proceed until a certain point at which it is halted by feedback^ In 
the main body of this paper, this is assumed to happen due to super- 
nova feedback. However, recent work has shown that actually ra- 
diative feedb ack may be the responsible mechanism, especially on 
small scales ^Thompson et alJboOSl ; iMurrav et al]l20ld: iFall et al.l 
l201CtlD~ale et alj|2012h . To verify to what extent this can influence 
the CFE, we focus on the truncation of star formation by radiation 
pressure in this appendix. 

We require pressure equilibrium between radiation pressure 
P ra d and turbulent pressure within a GMC Pgmc- The radiation 
pressure reads 



(1 + ' 



>trT) 



(1 + 0trKR,EGMc] 



47r ^GMC' : 



= (1 + 0trKB.£GMc) 



Tc ' 



(CI) 



where r ~ «.rSgmc i s me optical depth, cj>t t = fu/r ~ 0.2 
dKrumholz & Thompsonll2012h is a constant that indicates the frac- 
tion of infrared radiation that is trapped at an optical depth r = 1, 
T is the radiative flux, c is the speed of light, ^ ~ 0.3 m 2 s -3 
is the light-to-mass ratio dThompson et alfeOQSl ). M* is the stellar 



26 If this takes longer than the age at which the cluster formation efficiency 
is determined, star formation is considered to be still ongoing. If star for- 
mation is completed before feedback can halt star formation, the local SFE 
is optimal (see £|2.5t . 
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mass, 7?gmc is the GMC radius, e rat j is the SFE when star forma- 
tion is limited b y radiative feedback, and kr is the Rosselan d-mean 
dust opacity (cf. lThompson et al . 2005: lMurrav et all2 010). which 
for T < 200 K is expressed as 



Mi 



(I) 



C"SB 



(C2) 



4(TSB 4<TSB 

where ~ 2.4 x 10 -5 m 2 kg -1 K~ 2 is a proportionality constant 
and ctsb is the Stefan-Boltzmann constant. The turbulent pressure 
in a GMC is taken to be 

Pgmc = — GE GM c- (C3) 

By demanding pressure equilibrium P ra <j = Pgmc, equa- 
tions can be solved for the SFE e ra d. This gives a critical 
SFE above which radiation pressure inhibits further star formation: 



trad 



2<TSB 



/ 1 + 27rcG'0 tr K 2 S^ MC _ i 

0"SB 



(C4) 



In the context of our model, e ra d sets the maximum SFE if it 
is smaller than the core efficiency e ra d < e C orc (see §2.5\ . For 
Sgmc = Egmc = 100 Mq pc -2 , we indeed obtain e ra( j ~ 0.1 < 
e corc . The resulting CFE is shown for the entire density range in 
Figure lCTl and is compared to the fiducial model of H2.5\ Radiative 
feedback strongly limits the SFE and CFE until a density E g — 
E G mc> above which the GMC density changes from being approx- 
imately constant to Egmc = E g . As shown by equation dC4b . the 
SFE then scales as e ra( j oc Egmc oc E g if the interstellar medium 
(ISM) is optically thin, and as e ra( j oc Eq^jq oc E" 1 if the ISM 
is optically thick. This makes the SFE (and hence the CFE) peak at 
Eg ~ 10 3 ' 5 Mq pc -2 . At higher densities, the cruel cradle effect 
and radiative feedback in an optically thick ISM both decrease the 
CFE. 

While the details of the CFE-E g relation are affected by as- 
suming that radiative feedback halts star formation, the overall be- 
havior is comparable to the fiducial model of this paper. We con- 
clude that while the effects should be kept in mind (see §4. 1 1 and 
97.11 . there is no immediate need to adjust the fiducial model. 



APPENDIX D: THE BOUND STELLAR FRACTION AND 
THE STAR FORMATION EFFICIENCY 



In iKruiissen et alj 1 2012a ), we used a minimum spanning tree 



(MST , see e.g. lZahnlll971hlMaschberger et alfeofctlKirk & Mversl 
l201ll) to identify subclusters in the Bonnell et al. ( 2008) simulation 
of turbulent fragmentation. This technique is commonly used in the 
analysis of structure in star-forming regions due to its advantage of 
not imposing any symmetry on the data set. The MST is the unique 
connection of all points such that there are no closed loops and 
the total length of all connecting lines is minimized. Subclusters 
are then identified by adopting a break distance dbreak, which is 
the maximum distance between connected points. Any connecting 
lines longer than dbreak are removed, which breaks the MST of an 
entire star-forming region up into individua l subclusters. The fidu- 
cial break distance of dbrcak = 0.035 pc in lKruiissen et alJd2012al) 
was chosen to result in subclusters that are comparable to those 
identified by eye, but a range of dbrcak = 0.020-0.100 pc was 
used to validate the results. Additionally, a minimum of 12 stars 
per subcluster was postulated. 
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Figure CI. Influence of radiative feedback on the global cluster forma- 
tion efficiency (CFE) as a function of the gas surface density of the disc 
E g , assuming Q = 1.5 and writing f2 as a function of S g according to 
equation {9). The solid line shows the fiducial model with supernova feed- 
back, while the result for radiative feedback only is given by the dotted 
line. The dashed line represents the case in which star formation is halted 
by whichever mechanism generates the required feedback energy first, and 
thus combines both feedback processes. 



In this appendix, the subcluster data from IKruiissen et al.l 

are revisited to determine how the local CFE is related 
to the local SFE in a hierarchically structured star-forming re- 
gion. To this end, the stellar-to-total mass ratios within the half- 
mass radii of the subclusters represent the SFE e, while the frac- 
tion of the stars that is bound to the subcluster when ignoring the 
gravitational potential of the gas is considered to reflect the local 
CFE 7(e). The result is shown for break distances of dbrcak = 
{0.020, 0.040, 0.075, 0.100} pc in the left panel ofFigurelrJTI and 
includes the subclusters at times t± = 0.442 Myr (briefly after the 
onset of star formation) and = 0.641 Myr (after one global free- 
fall time). The figure also shows data points for which the calcu- 
lation was applied to the whole box of the simulation - the largest 
(and thus most gas-rich) scale on which the CFE can be obtained. 
Figure iDTI shows that the subclusters follow a single trend of 



7(e) = e, 



(Dl) 



independently of the time, break distance or location. The scatter 
around the relation is modest, about 30%. A continuous relation 
such as the figure seems to suggest, implies that the fraction of the 
stellar mass that ends up being bound is given by the local stel- 
lar mass fraction itself. There is no obvious reason as to why the 
7(e) relation would be linear, let alone why 7 and e are simply ap- 
proximately equal. Qualitatively speaking, a high local stellar mass 
fraction translates into a higher gravitational attraction of stars to- 
wards the stellar loci, where the resulting increase of the stellar 
density may l ead to the runaway acc retio n and gas deplet i on tha t 
was found by IKruiissen et al. I J2012ah and lGirichidisetalJ J2012h . 
But again, the details are not clear a priori, and it is evident that 
additional study is needed to further elucidate the 7(e) relation in a 
hierarchical interstellar medium. 

The 7(e) relation that is found in simulations of the influence 
of gas expulsion on star cluster survival with spherically symmet- 
ric initial conditions, in which the gas is modelle d with an analytic 
background potential dGever & Burkertl l200ll ; iBoilv & K roupa 
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Figure Dl. Left: Local cluster formation efficiency 7 as a function of the local star formation efficiency e in the Bonnell et al. 12008) simulation. The symbols 
indicate the values for subclusters that were identified in the simulation with different minimum spanning tree break distances - plus signs, diamonds, triangles 
and dots correspond to dbrcak = {0.020, 0.040, 0.075, O.fOO} pc, respectively, with big dots indicating the entire simulation (see text). The dashed line 
shows the 1:1 relatio n and the dotted curve gives a parameterized 7(e) relation that is typically found in equilibrium-type, spherically symmetric simulations 
of gas expulsion (e.g. lBaumgardt & Kroupa 2007). Right: Influence of the choice of 7(e) relation on the global cluster formation efficiency T as a function of 
the gas surface density of the disc E g , assuming Q = 1.5 and writing C as a function of E g according to equation 



l2003ah iGoodwin & Bastianl 120061: iBaumgardt & Kroupal 120071) . 
differs from the relation that is obtained here. In such simulations, 
the central concentration of the gas implies a threshold SFE be- 
low which a stellar cluster is unable to remain bound after gas 
expulsion. Because these are equilibrium models, such a gas-rich 
star-forming region is then also unable to produce locally bound 
agglomerates of stars, implying that below the threshold SFE they 
predict 7 = 0. This implies that the transition from a zero CFE to 
almost 100% is more su dden, and is generally found t o occur in the 
30-50% SFE range (e.g. lBaumgardt & Kroupall2007l) . It should be 
noted that this is exact ly the range that is n ot well-resolved by the 
above analysis of the Bonnell et al. (2008) simulation. However, 
the behavior of 7(e) at SFEs of e > 0.5 differs as well. Due to the 
sudden transition from 7 = to 7 = 1 in spherically symmetric 
simulations, they show an elevated j(e > 0.5) with respect to the 
relation of equation JDU . 

From a physical perspective it is preferable to adopt the hi- 
erarchical 7(e) relation over the spherically symmetric one, but to 
assess the impact of this choice, the influence of the two differ- 
ent 7(e) relations on the CFE is shown in the right-hand panel of 
Figure IdTI For this purpose, the 7(e) relation for spherically sym- 
metric clusters is parame terized with a simple function that reflects 
the typical form found in Baumg ardt & Kroupal (2007): 



7(e) = max[0, 1-4(1 - e) 3 ], 



(D2) 



where in the model calculation e is replaced by e/e cola as in equa- 
tion l(25). The right-hand panel of Figure iDTI shows T as a function 
of the gas surface density E g . The difference between CFEs for the 
hierarchical and spherically symmetric 7(e) relations is typically 
^ 0.6 dex. However, the steep transition from 7 = to 7 = 1 at 
intermediate SFEs causes T in the spherically symmetric case to be 
two times more sensitive to the feedback time-scale than when us- 
ing the hierarchical 7(e) relation. It is therefore not only physically 
preferable to adopt a gradual 7(e) relation, but also for the purpose 
of minimizing the model uncertainties. 



APPENDIX E: THE EFFICIENCY AND STOCHASTICITY 
OF THE TIDAL DISRUPTION OF STAR-FORMING 
REGIONS 

To determine the contribution of the cruel cradle effect to the clus- 
ter formation efficiency, it is important to determine the fraction / 
of the energy injected by tidal shocks into star-forming regions that 
is used to unbind the regions. As explained in ^2.71 this fraction is 
smaller than unity due to the energy dissipation by the cold gas. If 
this turbulent energy decay is assumed to be exponential, it follows 
that the disruption is only possible if a sequence of tidal pertur- 
bations has the appropriate strengths and time spacing to cause a 
runaway effect. Since this is a highly stochastic process, / should 
be determined by modeling a sequence of tidal perturbations in a 
Monte-Carlo simulation. 

Each Monte-Carlo experiment carried out here uses three ran- 
domly drawn variables: the time interval between successive en- 
counters tdiSi the impact parameter b, and the relative velocity v. 
For a stochastic process that occurs on a characteristic time-scale 
£enc, the probability distribution function (PDF) of the time interval 
t between successive events is 



dp 

dtdiff 



-*diff/*e 



(El) 



where t cnc is taken from equation l !34t . The PDF for the impact 
parameter follows from geometric considerations as 

dp 2b 



db bL 



(E2) 



with 6 m ax defined as in §277] The relative velocity of each en- 
counter is drawn from a Maxwellian distribution of relative veloci- 
ties with a dispersion of \/2<7 g : 



(E3) 



— = e s . 

dv 

The example for which the results are shown below was set 
up for solar neighbourhood-type conditions, i.e. £ g = 12 Mq pc 2 , 
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Q = 1.5 and fi = 2.6 x 10~ 2 Myr -1 , but a broad parameter 
range (cf. Table [TJ was considered in other experiments to verify 
the universality of the results. For a given sequence of encounters, 
the energy after each encounter is calculated by writing 



E = E. 



I'll'. 



(AE) C 



(E4) 



where E plev is the energy after the previous perturbation, (A_E) t id 
is the energy change due to the tidal shock, and (A_E)di ss is the en- 
ergy loss since the last shock due to dissipation. For a given over- 
density x, the tidal disruption equations from Sj2j7]are used to cal- 
culate the injected energy per encounter as 



(AB)tid 



16G 2 (Jclose 4>sh MMgMC Rh 4>a.d 



(E5) 



3b 4 V 2 

in which M is the mass of the perturbed region, Rh its half-mass 
radius, and all other variables are defined as previously. For the 
example below, M = Mgmc and Rh = (3Mgmc /87rp g ) 1//3 . The 
density p g is set by x and the density of the ambient interstellar 
medium pism as giv en by equation (5) . The initial energy of the 
region Eq is set by a lPlummeg dl91 ll) profile as in equation l |29l >, 
and the energy that is dissipated between successive encounters is 
written as an exponential decay: 



(AE) dis 



\E V 



Eu\ (l 



-*dif£ /*diss"\ 



(E6) 



with tdiss given by equation l |32t . Note that many of the details 
of these equations (such as the mass and radius of the region) are 
actually irrelevant to the question at hand, which strictly involves 
the fractional energy loss due to dissipation and the stochasticity 
of the disruption process. Multiplication of equation l lE5t by a ran- 
dom factor may change the lifetimes of the modelled star-forming 
regions, but does not necessarily alter the fraction of the injected 
energy that is used to unbind the region. 

For each Monte-Carlo experiment, the region is subjected to 
tidal perturbations until either it is disrupted (i.e. E > 0) or a time 
t is reached. The fraction of the energy that was used to unbind the 
region then follows as 



/mc = 



E — En 



(A£) tid ,, 



(E7) 



where the subscript 'MC indicates the value is obtained from 
the Monte-Carlo experiment, (A_E) t id,tot = (A-E)tid* an d the 
summation is over all perturbations. For each experiment, the total 
disruption time i CC e,MC is recorded: 

E 



to 



E — Eo 



(E8) 



where t b s = min {idis, t} and tdis is the time of the tidal shock 
that causes E > 0. The definition of equation dE8t is appropriate 
for regions that are disrupted before time t (for which t Q b s = tdis) 
and corrects the disruption time-scale for any excess energy, but it 
also gives an estimate for regions with disruption times tdis > t (for 
which t bs = £). Without energy dissipation, the total disruption 
time would read 



t' 



t' 



So 



(E9) 



where t' ohs and (AE)' tid tot are set at the time of the shock that 
causes j (AE)' tld tot > \Eo\. The dissipational energy loss can 
then be parameterized by decreasing the total injected energy by 



a fixed factor fa t , i.e. writing t£ 



t 1 



: I /fit and choosing 



/At such that iccc.MC ~ iccc,MC- This formulation is useful when 
dealing with ensembles of star-forming regions as in the context of 




Figure El. Fraction /mc °f tidally injected energy that is used in unbind- 
ing the perturbed region instead of being dissipated away as a function of 
the time of observation t and the overdensity with respect to the ambient in- 
terstellar medium x. The grayscale indicates the median value of /mc ou t 
of a sample of 30,000 Monte-Carlo experiments per grid point. The thick 
solid line indicates where the median total disruption time due to the cruel 
cradle effect t CCB mc is equal to t, whereas the thick dotted line represents 
the loci where the median f/ ce MC (i.e. without energy dissipation) equals 
t, and the thick dashed line shows the same for t flt .»„ and f« t = 0.7. As 

the figure shows, the agreement between f- CC c,MC an d *ccc MC ' s reason " 
able, and the difference between them never exceeds 0.3 dex in x across the 
full parameter space probed in the performed Monte-Carlo experiments. 



this paper. We determine the constant fat from the results of the 
Monte-Carlo experiments below. 

The described Monte-Carlo experiment is executed 30,000 
times for each combination of the times 1 < t/Myr < 100 and 
overdensities 10° < x < 10 4 . The median value of / according to 
equation l(E7]l is shown in Figure IET1 as a function of t and x. The 
figure shows that the energy loss is typically largest for high over- 
densities and late times. This is not surprising, because high-density 
regions have short dissipation time-scales and take the longest to 
be tidally disrupted. Additionally, the energy loss should be low at 
young ages, because the energy simply has not had the time to be 
dissipated. 

Because the contribution of the cruel cradle effect to the CFE 
is formulated in equation J40t > as a threshold value a; ccc above 
which star-forming regions survive until time t, Figure lETI includes 
the line that satisfies median(t CCCj Mc) = t. For comparison, the 
line where the median of ^ccc.mc equals t is also shown, assuming 
/tit = 0.7. For this choice of fht, both lines never deviate by more 
than 0.3 dex, indicating that it provides a good representation for 
the fraction of the tidally injected energy that is used to unbind the 
region. As the discrepancy between the solid and dashed lines il- 
lustrates, this slightly underestimates the impact of the cruel cradle 
effect at low x and t, while slightly overestimating it at high x and 
t. This reflects the increase of the energy dissipation with x and 
t. However, the differences are so minor that / = 0.7 is adopted 
throughout this paper. 

The reason that the tidally injected energy is mainly used for 
unbinding a region (i.e. / > 0.5) if the disruption time t C cc,MC < t 
is partially caused by a selection effect. Due to the energy dissipa- 
tion, a sequence of tightly-spaced tidal perturbations is required to 



32 J. M. D. Kruijssen 



unbind a gas-rich region. Because the energy from such a sequence 
of perturbations is not dissipated, any region that is disrupted be- 
fore time t will have a high efficiency factor /. Conversely, regions 
that survive have not had a runaway sequence of tidal shocks and 
radiated away most of the energy, thus giving a low value of /. 
In a statistical ensemble of tidal shocks, the probability of having 
a sequence of encounters that successfully disrupts a region is set 
by its disruption time-scale. This probability is close to unity if 
tcce.MC ^ t, and hence / will be close to unity as well. 



REFERENCES 

Abadi M. G., Navarro J. R, Steinmetz M., Eke V. R., 2003, ApJ, 591, 499 
Adamo A., Ostlin G., Zackrisson E., 2011, MNRAS, 417, 1904 
Adamo A., Ostlin G., Zackrisson E., Hayes M., Cumming R. J., Micheva 

G., 2010, MNRAS, 407, 870 
Adams F. C, 2000, ApJ, 542, 964 

Adams F. C, Hollenbach D., Laughlin G., Gorti U., 2004, ApJ, 61 1, 360 
Aguilar L. A., White S. D. M., 1985, ApJ, 295, 374 
Allen L., Megeath S. T., Gutermuth R., Myers P. C., Wolk S., Adams F. C, 
Muzerolle J., Young E., Pipher J. L., 2007, Protostars and Planets V, 361 
Alves D. R., Nelson C. A., 2000, ApJ, 542, 789 

Arce H. G., Borkin M. A., Goodman A. A., Pineda J. E., Halle M. W., 

2010, ApJ, 715, 1170 
Assmann P., Fellhauer M., Kroupa P., Brims R. C., Smith R., 2011, 

MNRAS, 415, 1280 
Balbus S. A., 1988, ApJ, 324, 60 
Bastian N., 2008, MNRAS, 390, 759 

Bastian N., Adamo A., Gieles M., Silva- Villa E., Lamers H. J. G. L. M., 
Larsen S. S., Smith L. J., Konstantopoulos I. S., Zackrisson E., 2012, 
MNRAS, 419, 2606 

Bastian N., Covey K. R., Meyer M. R., 2010, ARA&A, 48, 339 

Bastian N., Gieles M., Lamers H. J. G. L. M., Scheepmaker R. A., De 
Grijs R., 2005, A&A, 431, 905 

Bastian N., Goodwin S. P., 2006, MNRAS, 369, L9 

Baumgardt H., Kroupa P., 2007, MNRAS, 380, 1589 

Baumgardt H., Makino J., 2003, MNRAS, 340, 227 

Bergvall N., Masegosa J., Ostlin G., Cernicharo J., 2000, A&A, 359, 41 

Bertoldi F, McKee C. R, 1992, ApJ, 395, 140 

Bigiel R, Leroy A., Walter R, Brinks E., de Blok W. J. G., Madore B., 

Thornley M. D., 2008, AJ, 136, 2846 
Binney J., Tremaine S., 1987, Galactic dynamics. Princeton, NJ, Princeton 

University Press, 1987, 747 pp. 
Blaauw A., 1964, ARA&A, 2, 213 

Bland-Hawthorn J., Krumholz M. R., Freeman K., 2010, ApJ, 713, 166 
Boily C. M., Kroupa P., 2003a, MNRAS, 338, 665 
— , 2003b, MNRAS, 338, 673 

Bolatto A. D., Leroy A. K, Rosolowsky E., Walter R, Blitz L., 2008, ApJ, 
686, 948 

Bonnell I. A., Clark P., Bate M. R., 2008, MNRAS, 389, 1556 

Bower R. G., Benson A. J., Malbon R., Helly J. C, Frenk C. S., Baugh 

C. M., Cole S., Lacey C. G., 2006, MNRAS, 370, 645 
Braun R., 1991, ApJ, 372, 54 

Braun R., Thilker D. A., Walterbos R. A. M., Corbelli E., 2009, ApJ, 695, 
937 

Bressert E., Bastian N, Evans C. J., Sana H, Henault-Brunet V., Good- 
win S. P., Parker R. J., Gieles M., Bestenlehner J. M., Vink J. S., Taylor 
W. D., Crowther P. A., Longmore S. N., Grafener G., Maiz Apellaniz J., 
de Koter A., Cantiello M., Kruijssen J. M. D., 2012a, A&A, 542, A49 

Bressert E., Bastian N., Gutermuth R., Megeath S. T., Allen L., Evans II 
N. J., Rebull L. M., Hatchell J., Johnstone D., Bourke T. L., Cieza L. A., 
Harvey P. M., Merin B., Ray T. P., Tothill N. F. H., 2010, MNRAS, 409, 
L54 

Bressert E., Ginsburg A., Bally J., Battersby C, Longmore S. N., Testi L., 

2012b, ApJ in press 
Bressert E., Testi L., Facchini A., Molinari S., Schisano E., Pezzuto S., 



Pestalozzi M., Sadavoy S., Andre P., Benedettini M., Di Francesco J., 
Elia D., Fallscheer C, Hennemann M., Hill T, Nguyen Luong Q., Motte 
R, Schneider N., Spinoglio L., Ward-Thompson D., White G. J., 2012c, 
A&A submitted 
Buckle J. V., Richer J. S., Davis C. J., 2012, MNRAS, 423, 1 127 
Bundy K., Ellis R. S., Conselice C. J., Taylor J. E., Cooper M. C, Willmer 
C. N. A., Weiner B. J., Coil A. L., Noeske K. G., Eisenhardt P. R. M., 
2006, ApJ, 651, 120 
Burkert A., Tremaine S., 2010, ApJ, 720, 516 

Cappellari M., McDermid R. M., Alatalo K., Blitz L., Bois M., Bour- 
naud R, Bureau M., Crocker A. R, Davies R. L., Davis T. A., de Zeeuw 
P. T, Due P.-A., Emsellem E., Khochfar S., Krajnovic D., Kuntschner 
H, Lablanche P.-Y., Morganti R., Naab T., Oosterloo T., Sarzi M., Scott 
N., Serra P., Weijmans A.-M., Young L. M., 2012, Nature, 484, 485 

Chabrier G., 2003, PASP, 115, 763 

Chapman S. C, Smail I., Windhorst R., Muxlow T, Ivison R. J., 2004, 
ApJ, 611,732 

Chemin L., Carignan C, Drouin N., Freeman K. C, 2006, AJ, 132, 2527 

Chemin L., Carignan C, Foster T., 2009, ApJ, 705, 1395 

Cimatti A., Daddi E., Renzini A., Cassata P., Vanzella E., Pozzetti L., Cris- 

tiani S., Fontana A., Rodighiero G., Mignoli M., Zamorani G., 2004, 

Nature, 430, 184 

Clarke C. J., Bonnell I. A., Hillenbrand L. A., 2000, Protostars and Planets 
IV, 151 

Cole S., Lacey C. G., Baugh C. M., Frenk C. S., 2000, MNRAS, 319, 168 
Combes F, Becquaert J.-R, 1997, A&A, 326, 554 

Cook D. O., Seth A. C, Dale D. A., Johnson L. C, Weisz D. R., Fouesneau 
M., Olsen K. A. G., Engelbracht C. W., Dalcanton J. J., 2012, ApJ, 751, 
100 

Cottaar M., Meyer M. R., Andersen M., Espinoza P., 2012, A&A, 539, 
A5 

da Silva R. L., Fumagalli M., Krumholz M., 2012, ApJ, 745, 145 

Daddi E., Bournaud R, Walter R, Dannerbauer H, Carilli C. L., Dickinson 
M., Elbaz D., Morrison G. E., Riechers D., Onodera M., Salmi R, Krips 
M., Stern D., 2010, ApJ, 713, 686 

Daddi E., Dickinson M., Morrison G., Chary R., Cimatti A., Elbaz D., 
Frayer D., Renzini A., Pope A., Alexander D. M., Bauer F. E., Giavalisco 
M., Huynh M., Kurk J., Mignoli M., 2007, ApJ, 670, 156 

Daddi E., Renzini A., Pirzkal N., Cimatti A., Malhotra S., Stiavelli M., Xu 
C, Pasquali A., Rhoads J. E., Brasa M., di Serego Alighieri S., Ferguson 
H. C, Koekemoer A. M., Moustakas L. A., Panagia N., Windhorst R. A., 
2005, ApJ, 626, 680 

Dalcanton J. J., Williams B. F., Lang D., Lauer T. R., Kalirai J. S., Seth 
A. C, Dolphin A., Rosenfield P., Weisz D. R., Bell E. R, Bianchi L. C, 
Boyer M. L., Caldwell N., Dong R, Dorman C. E., Gilbert K. M., Gi- 
rardi L., Gogarten S. M., Gordon K. D., Guhathakurta P., Hodge P. W, 
Holtzman J. A., Johnson L. C, Larsen S. S., Lewis A., Melbourne J. L., 
Olsen K. A. G., Rix H.-W., Rosema K., Saha A., Sarajedini A., Skillman 
E. D., Stanek K. Z., 2012, ApJS, 200, 18 

Dale J. E., Bonnell I. A., Clarke C. J., Bate M. R., 2005, MNRAS, 358, 
291 

Dale J. E., Ercolano B., Bonnell I. A., 2012, MNRAS, 424, 377 

Dalla Vecchia C, Schaye J., 2008, MNRAS, 387, 1431 

de Juan Ovelar M., Kruijssen J. M. D., Bressert E., Testi L., Bastian N., 

Canovas H., 2012, A&A submitted 
De Silva G. M., Freeman K. C, Asplund M., Bland-Hawthorn J., Bessell 

M. S., Collet R., 2007, AJ, 133, 1 161 
De Silva G. M., Sneden C, Paulson D. B., Asplund M., Bland-Hawthorn 

J., Bessell M. S., Freeman K. C, 2006, AJ, 131, 455 
de Wit W. J., Testi L., Palla R, Zinnecker H, 2005, A&A, 437, 247 
de Zeeuw P. T., Hoogerwerf R., de Bruijne J. H. J., Brown A. G. A., 

Blaauw A., 1999, AJ, 117, 354 
Dib S., Bell E., Burkert A., 2006, ApJ, 638, 797 

Dicaire I., Carignan C, Amram P., Marcelin M., Hlavacek-Larrondo J., 
de Denus-Baillargeon M.-M., Daigle O., Hernandez O., 2008, A J, 135, 
2038 

Dobbs C. L., Burkert A., Pringle J. E., 2011a, MNRAS, 417, 1318 
— , 2011b, MNRAS, 413, 2935 



The fraction of star formation in bound clusters 33 



Downes D., Solomon P. M., 1998, ApJ, 507, 615 
Dubois Y., Teyssier R., 2008, A&A, 477, 79 
Dukes D., Krumholz M. R., 2012, ApJ, 754, 56 
Efremov Y. N., Elmegreen B. G., 1998, MNRAS, 299, 588 
Efstathiou G., 2000, MNRAS, 317, 697 

Eggen O. J., Lynden-Bell D., Sandage A. R., 1962, ApJ, 136, 748 

Elmegreen B. G., 1983, MNRAS, 203, 1011 

— , 2000, ApJ, 530, 277 

— , 2002, ApJ, 577, 206 

— , 2007, ApJ, 668, 1064 

— , 2008, ApJ, 672, 1006 

—,2010, ApJ, 712, L184 

Elmegreen B. G., Efremov Y. N., 1997, ApJ, 480, 235 
Elmegreen B. G., Elmegreen D. M., 2001, AJ, 121, 1507 
— , 2005, ApJ, 627, 632 

Elmegreen B. G., Falgarone E., 1996, ApJ, 471, 816 

Elmegreen B. G., Hunter D. A., 2010, ApJ, 712, 604 

Elmegreen B. G., Lada C. J., 1977, ApJ, 214, 725 

Elmegreen B. G., Scalo J., 2004, ARA&A, 42, 211 

Elmegreen D. M., Elmegreen B. G, Marcus M. T., Shahinyan K., Yau A., 

Petersen M., 2009, ApJ, 701, 306 
Erb D. K., Steidel C. C, Shapley A. E., Pettini M., Reddy N. A., Adel- 

berger K. L., 2006, ApJ, 646, 107 
Evans II N. J., Dunham M. M., J0rgensen J. K., Enoch M. L., Merin B., 

van Dishoeck E. E, Alcala J. M., Myers P. C., Stapelfeldt K. R., Huard 

T. L., Allen L. E., Harvey P. M., van Kempen T., Blake G. A., Koerner 

D. W., Mundy L. G., Padgett D. L., Sargent A. I., 2009, ApJS, 181, 321 
Fall S. M., Krumholz M. R., Matzner C. D., 2010, ApJ, 710, L142 

Fall S. M., Zhang Q., 2001, ApJ, 561, 751 

Federrath C, Klessen R. S., Schmidt W., 2008, ApJ, 688, L79 

Federrath C, Roman-Duval J., Klessen R. S., Schmidt W., Mac Low M.- 
M.,2010, A&A, 512, A81 

Fedotov K, Gallagher S. C., Konstantopoulos I. S., Chandar R., Bastian 
N., Charlton J. C, Whitmore B., Trancho G., 201 1, AJ, 142, 42 

Ford H. A., Bregman J. N., 2012, ApJ submitted, ArXiv: 1205. 1066 

Forster Schreiber N. M., Genzel R., Bouche N., Cresci G., Davies R., 
Buschkamp P., Shapiro K., Tacconi L. J., Hicks E. K. S., Genel S., Shap- 
ley A. E., Erb D. K., Steidel C. C, Lutz D., Eisenhauer F, Gillessen S., 
Sternberg A., Renzini A., Cimatti A., Daddi E., Kurk J., Lilly S., Kong 
X., Lehnert M. D., Nesvadba N., Verma A., McCracken H., Arimoto N., 
Mignoli M., Onodera M., 2009, ApJ, 706, 1364 

Genzel R., Burkert A., Bouche N., Cresci G., Forster Schreiber N. M., 
Shapley A., Shapiro K., Tacconi L. J., Buschkamp P., Cimatti A., Daddi 

E. , Davies R., Eisenhauer F., Erb D. K., Genel S., Gerhard O., Hicks E., 
Lutz D., Naab T., Ott T., Rabien S., Renzini A., Steidel C. C, Sternberg 
A., Lilly S. J., 2008, ApJ, 687, 59 

Genzel R., Tacconi L. J., Gracia-Carpio J., Sternberg A., Cooper M. C, 
Shapiro K., Bolatto A., Bouche N., Bournaud E, Burkert A., Combes 

F. , Comerford J., Cox P., Davis M., Schreiber N. M. E, Garcia-Burillo 
S., Lutz D., Naab T, Neri R., Omont A., Shapley A., Weiner B., 2010, 
MNRAS, 407, 2091 

Geyer M. P., Burkert A., 2001, MNRAS, 323, 988 

Gieles M., 2010, in Astronomical Society of the Pacific Conference Se- 
ries, Vol. 423, Galaxy Wars: Stellar Populations and Star Formation in 
Interacting Galaxies, B. Smith, J. Higdon, S. Higdon, & N. Bastian, ed., 
p. 123 

Gieles M., Bastian N., Lamers H. J. G. L. M., Mout J. N., 2005, A&A, 
441, 949 

Gieles M., Heggie D. C, Zhao H., 201 1, MNRAS, 413, 2509 

Gieles M., Larsen S. S., Scheepmaker R. A., Bastian N., Haas M. R., 

Lamers H. J. G. L. M., 2006a, A&A, 446, L9 
Gieles M., Portegies Zwart S. F, 201 1, MNRAS, 410, L6 
Gieles M., Portegies Zwart S. F, Baumgardt H., Athanassoula E., Lamers 

H. J. G. L. M., Sipior M., Leenaarts J., 2006b, MNRAS, 371, 793 
Girichidis P., Federrath C, Banerjee R., Klessen R. S., 2012, MNRAS, 

420, 613 

Gnedin O. Y, Lee H. M., Ostriker J. P., 1999, ApJ, 522, 935 
Gnedin O. Y, Ostriker J. P., 1997, ApJ, 474, 223 



Goddard Q. E., Bastian N., Kennicutt R. C, 2010, MNRAS, 405, 857 

Goodwin S. P., 1997, MNRAS, 284, 785 

Goodwin S. P., Bastian N., 2006, MNRAS, 373, 752 

Gratton R., Sneden C, Carretta E., 2004, ARA&A, 42, 385 

Guo Q., White S., Boylan-Kolchin M., De Lucia G, Kauffmann G., Lem- 

son G., Li C, Springel V., Weinmann S., 201 1, MNRAS, 413, 101 
Gutermuth R. A., Megeath S. T, Myers P. C, Allen L. E., Pipher J. L., 

Fazio G. G., 2009, ApJS, 184, 18 
Gutermuth R. A., Pipher J. L., Megeath S. T, Myers P. C, Allen L. E., 

Allen T. S., 2011, ApJ, 739, 84 
Hansen C. E., Klein R. I., McKee C. E, Fisher R. T., 2012, ApJ, 747, 22 
Heiles C, Troland T. H., 2003, ApJ, 586, 1067 

Heyer M., Krawczyk C, Duval J., Jackson J. M., 2009, ApJ, 699, 1092 

Hill T., Motte E, Didelon P., White G. J., Marston A. P., Nguyen Luong 
Q., Bontemps S., Andre P., Schneider N., Hennemann M., Sauvage M., 
Di Francesco J., Minier V., Anderson L. D., Bernard J. P., Elia D., Griffin 
M. J., Li J. Z., Peretto N., Pezzuto S., Polychroni D., Roussel H., Rygl 
K. L. J., Schisano E., Sousbie T, Testi L., Thompson D. W., Zavagno A., 
2012, A&A, 542, A114 

Hills J. G., 1980, ApJ, 235, 986 

Hopkins A. M., Beacom J. F, 2006, ApJ, 651, 142 

Hopkins P. F, Lauer T. R., Cox T. J., Hernquist L., Kormendy J., 2009, 
ApJS, 181,486 

Hopkins P. F, Quataert E., Murray N., 2012, MNRAS, 421, 3488 
Huff E. M., Stahler S. W., 2006, ApJ, 644, 355 

Ilbert O., Salvato M., Le Floc'h E., Aussel H., Capak P., McCracken 

H. J., Mobasher B., Kartaltepe J., Scoville N., Sanders D. B., Arnouts 
S., Bundy K., Cassata P., Kneib J. -P., Koekemoer A., Le Fevre O., Lilly 
S., Surace J., Taniguchi Y, Tasca L., Thompson D., Tresse L., Zamojski 
M., Zamorani G, Zucca E., 2010, ApJ, 709, 644 

Johnson L. C, Seth A. C, Dalcanton J. J., Caldwell N., Fouesneau M., 
Gouliermis D. A., Hodge P. W., Larsen S. S., Olsen K. A. G., San Roman 

I. , Sarajedini A., Weisz D. R., Williams B. F, Beerman L. C, Bianchi L., 
Dolphin A. E., Girardi L., Guhathakurta P., Kalirai J., Lang D., Monach- 
esi A., Nanda S., Rix H.-W., Skillman E. D., 2012, ApJ, 752, 95 

Jordan A., McLaughlin D. E., Cote P., Ferrarese L., Peng E. W., Mei S., 

Villegas D., Merritt D., Tonry J. L., West M. J., 2007, ApJS, 171, 101 
J0rgensen J. K., Johnstone D., Kirk H., Myers P. C, Allen L. E., Shirley 

Y L., 2008, ApJ, 683, 822 
Joung M. K. R., Mac Low M.-M., 2006, ApJ, 653, 1266 
Kauffmann G., Heckman T. M., White S. D. M., Chariot S., Tremonti C, 

Peng E. W., Seibert M., Brinkmann J., Nichol R. C, SubbaRao M., York 

D., 2003, MNRAS, 341, 54 
Kauffmann G., White S. D. M., Guiderdoni B., 1993, MNRAS, 264, 201 
Kennicutt Jr. R. C, 1989, ApJ, 344, 685 
— , 1998, ApJ, 498, 541 
Kirk H., Myers P. C, 201 1, ApJ, 727, 64 

Konstantopoulos I. S., Bastian N., Smith L. J., Westmoquette M. S., Tran- 
cho G., Gallagher III J. S., 2009, ApJ, 701, 1015 

Korpi M. J., Brandenburg A., Shukurov A., Tuominen I., Nordlund A., 
1999, ApJ, 514, L99 

Krause M., Fierlinger K., Diehl R., Burkert A., Voss R., Ziegler U., 2012, 
A&A submitted, ArXiv: 1207.7231 

Kravtsov A. V., Gnedin O. Y, 2005, ApJ, 623, 650 

Kroupa P., 2001, MNRAS, 322, 231 

Kruijssen J. M. D., 2009, A&A, 507, 1409 

— , 2011, in Stellar Clusters & Associations: A RIA Workshop on Gaia, 
pp. 137-141 

Kruijssen J. M. D., Cooper A. P., 2012, MNRAS, 420, 340 

Kruijssen J. M. D., Maschberger T., Moeckel N., Clarke C. J., Bastian N., 

Bonnell I. A., 2012a, MNRAS, 419, 841 
Kruijssen J. M. D., Pelupessy F. I., Lamers H. J. G. L. M., Portegies Zwart 

S. E, Bastian N, Icke V., 2012b, MNRAS, 421, 1927 
Kruijssen J. M. D., Pelupessy F. I., Lamers H. J. G. L. M., Portegies Zwart 

S. E, Icke V., 2011, MNRAS, 414, 1339 
Krumholz M. R., Klein R. I., McKee C. F, 2012, ApJ, 754, 71 
Krumholz M. R., Klein R. I., McKee C. F, Offner S. S. R., Cunningham 

A. J., 2009, Science, 323, 754 



34 J. M. D. Kruijssen 



Kramholz M. R., McKee C. F., 2005, ApJ, 630, 250 
Krumholz M. R., Tan J. C, 2007, ApJ, 654, 304 

Kramholz M. R., Thompson T. A., 2012, ApJ submitted, 

ArXiv: 1203.2926 
Kundic T, Ostriker J. R, 1995, ApJ, 438, 702 

Lada C. J., 1987, in IAU Symposium, Vol. 115, Star Forming Regions, 

M. Peimbert & J. Jugaku, ed., pp. 1-17 
Lada C. J., Lada E. A., 2003, ARA&A, 41, 57 
Lada C. J., Margulis M., Dearborn D., 1984, ApJ, 285, 141 
Larsen S. S., 2004, A&A, 416, 537 
— 2009, A&A, 494, 539 

Larsen S. S., Brodie J. P., Huchra J. P., Forbes D. A., Grillmair C. J., 2001, 

A J, 121, 2974 
Larsen S. S., Richtler T, 2000, A&A, 354, 836 
Larson R. B., 1974, MNRAS, 169, 229 
— , 1981, MNRAS, 194, 809 

Lilly S. J., Eales S. A., Gear W. K. P., Hammer F, Le Fevre O., Crampton 

D. , Bond J. R., Dunne L., 1999, ApJ, 518, 641 

Longmore S. N., Bally J., Testi L., Purcell C. R., Walsh A. J., Bressert E., 
Pestalozzi M., Molinari S., Ott J., Cortese L., Battersby C, Murray N., 
Lee E., Kruijssen J. M. D., 2012a, submitted 

Longmore S. N., Rathborne J., Bastian N., Alves J., Ascenso J., Bally 
J., Testi L., Longmore A., Battersby C, Bressert E., Purcell C, Walsh 
A., Jackson J., Foster J., Molinari S., Meingast S., Amorim A., Lima J., 
Marques R., Moitinho A., Pinhao J., Rebordao J., Santos F. D., 2012b, 
ApJ, 746, 117 

Lundgren A. A., Olofsson H., Wiklind T, Rydbeck G., 2004, A&A, 422, 
865 

Mac Low M.-M., Ferrara A., 1999, ApJ, 513, 142 

Mac Low M.-M., Klessen R. S., 2004, Reviews of Modern Physics, 76, 
125 

Madau P., Ferguson H. C, Dickinson M. E., Giavalisco M., Steidel C. C, 

Frachter A., 1996, MNRAS, 283, 1388 
Madau P., Pozzetti L., Dickinson M., 1998, ApJ, 498, 106 
Majewski S. R., Nidever D. L., Smith V. V., Damke G. J., Kunkel W. E., 

Patterson R. J., Bizyaev D., Garcia Perez A. E., 2012, ApJ, 747, L37 
Marchesini D., van Dokkum P. G., Forster Schreiber N. M., Franx M., 

Labbe I., Wuyts S., 2009, ApJ, 701, 1765 
Marigo P., Girardi L., Bressan A., Groenewegen M. A. T., Silva L., 

Granato G. L., 2008, A&A, 482, 883 
Martin C. L., Kennicutt Jr. R. C, 2001, ApJ, 555, 301 
Maschberger T, Clarke C. J., Bonnell I. A., Kroupa P., 2010, MNRAS, 

404, 1061 

Matzner C. D., McKee C. F, 2000, ApJ, 545, 364 
McKee C. F, Ostriker E. C, 2007, ARA&A, 45, 565 
McKee C. F, Ostriker J. P., 1977, ApJ, 218, 148 
McKee C. F, Tan J. C, 2003, ApJ, 585, 850 
McKee C. F, Williams J. P., 1997, ApJ, 476, 144 

Megeath S. T., Allen L. E., Gutermuth R. A., Pipher J. L., Myers P. C, 
Calvet N., Hartmann L., Muzerolle J., Fazio G. G., 2004, ApJS, 154, 
367 

Moeckel N., Holland C., Clarke C. J., Bonnell I. A., 2012, 

MNRAS in press, ArXiv:1205.1677 
Molina F. Z., Glover S. C. O., Federrath C, Klessen R. S., 2012, 

MNRAS, 423, 2680 
Muratov A. L., Gnedin O. Y., 2010, ApJ, 718, 1266 
Murray N., Quataert E., Thompson T. A., 2010, ApJ, 709, 191 
Naab T., Johansson P. H., Ostriker J. P., Efstathiou G, 2007, ApJ, 658, 

710 

Nakamura F, Li Z.-Y., 2007, ApJ, 662, 395 
Navarro J. F, Steinmetz M., 2000, ApJ, 538, 477 

Nieten C, Neininger N, Guelin M., Ungerechts H., Lucas R., Berkhuijsen 

E. M., Beck R., Wielebinski R., 2006, A&A, 453, 459 

Offner S. S. R., Hansen C. E., Kramholz M. R., 2009, ApJ, 704, L124 
Olczak C, Pfalzner S., Spurzem R., 2006, ApJ, 642, 1 140 
Ossenkopf V., Mac Low M.-M., 2002, A&A, 390, 307 
Ostlin G., Amram P., Bergvall N., Masegosa J., Boulesteix J., Marquez I., 
2001, A&A, 374, 800 



Ostlin G., Amram P., Masegosa J., Bergvall N, Boulesteix J., 1999, 

A&AS, 137,419 
Padoan P., Nordlund A., 2002, ApJ, 576, 870 
—,2011, ApJ, 730, 40 

Padoan P., Nordlund A., Jones B. J. T, 1997, MNRAS, 288, 145 
Parker R. J., Goodwin S. P., 2007, MNRAS, 380, 1271 
Parker R. J., Meyer M. R., 2012, MNRAS accepted, ArXiv:1208.0335 
Parker R. J., Quanz S. P., 2012, MNRAS, 419, 2448 
Parmentier G, Goodwin S. P., Kroupa P., Baumgardt H., 2008, ApJ, 678, 
347 

Pelupessy F. I., Portegies Zwart S., 2012, MNRAS, 420, 1503 

Peng E. W., Cote P., Jordan A., Blakeslee J. P., Ferrarese L., Mei S., West 

M. J., Merritt D., Milosavljevic M., Tonry J. L., 2006, ApJ, 639, 838 
Peng E. W., Jordan A., Cote P., Takamiya M., West M. J., Blakeslee J. P., 

Chen C.-W., Ferrarese L., Mei S., Tonry J. L., West A. A., 2008, ApJ, 

681, 197 

Perez-Gonzalez P. G., Rieke G. H., Egami E., Alonso-Herrero A., Dole 
H., Papovich C, Blaylock M., Jones J., Rieke M., Rigby J., Barmby P., 
Fazio G. G., Huang J., Martin C, 2005, ApJ, 630, 82 

Pfalzner S., 2009, A&A, 498, L37 

Pflamm-Altenburg J., Weidner C, Kroupa P., 2007, ApJ, 671, 1550 
PlummerH. C, 1911, MNRAS, 71, 460 

Portegies Zwart S. F, McMillan S. L. W., Gieles M., 2010, ARA&A, 48, 
431 

Pozzetti L., Bolzonella M., Lamareille F, Zamorani G., Franzetti P., Le 
Fevre O., Iovino A., Temporin S., Ilbert O., Arnouts S., Chariot S., 
Brinchmann J., Zucca E., Tresse L., Scodeggio M., Guzzo L., Bottini D., 
Garilli B., Le Bran V., Maccagni D., Picat J. P., Scaramella R., Vettolani 
G., Zanichelli A., Adami C, Bardelli S., Cappi A., Ciliegi P., Contini 
T., Foucaud S., Gavignaud I., McCracken H. J., Marano B., Marinoni 
C, Mazure A., Meneux B., Merighi R., Paltani S., Pello R., Polio A., 
Radovich M., Bondi M., Bongiorno A., Cucciati O., de la Torre S., Gre- 
gorini L., Mellier Y, Merluzzi P., Vergani D., Walcher C. J., 2007, A&A, 
474, 443 

Prieto J. L., Gnedin O. Y, 2008, ApJ, 689, 919 
Rand R. J., 1993, ApJ, 410, 68 

Randich S., Sestito P., Primas F, Pallavicini R., Pasquini L., 2006, A&A, 
450, 557 

Reddy N. A., Steidel C. C, 2009, ApJ, 692, 778 

Rieke G. H., Cutri R. M., Black J. H., Kailey W. F, McAlary C. W., Lebof- 

sky M. J., Elston R., 1985, ApJ, 290, 1 16 
Robertson B., Bullock J. S., Font A. S., Johnston K. V., Hernquist L., 2005, 

ApJ, 632, 872 

Rochau B., Brandner W., Stolte A., Gennaro M., Gouliermis D., Da Rio 
N., Dzyurkevich N, Henning T., 2010, ApJ, 716, L90 

Ryder S. D., Staveley-Smith L., Malin D., Walsh W., 1995, AJ, 109, 1592 

Sakamoto K, Ho P. T. P., Peck A. B., 2006, ApJ, 644, 862 

Sales L. V., Navarro J. F, Schaye J., Dalla Vecchia C, Springel V., Booth 
C. M., 2010, MNRAS, 409, 1541 

Salpeter E. E., 1955, ApJ, 121, 161 

Scally A., Clarke C, 2001, MNRAS, 325, 449 

Scally A., Clarke C, McCaughrean M. J., 2005, MNRAS, 358, 742 

Schmidt M., 1959, ApJ, 129, 243 

Schuster K. F, Kramer C, Hitschfeld M., Garcia-Burillo S., Mookerjea 

B., 2007, A&A, 461, 143 
Scoville N. Z., Sargent A. I., Sanders D. B., Soifer B. T., 1991, ApJ, 366, 

L5 

Searle L., Zinn R., 1978, ApJ, 225, 357 

Shapiro K. L., Genzel R., Forster Schreiber N. M., 2010, MNRAS, 403, 
L36 

Silk J., 1997, ApJ, 481, 703 

Silva- Villa E., Larsen S. S., 2010, A&A, 516, A10+ 
—,2011, A&A, 529, A25+ 

Smith L. J., Bastian N., Konstantopoulos I. S., Gallagher III J. S., Gieles 
M., de Grijs R., Larsen S. S., O'Connell R. W., Westmoquette M. S., 
2007, ApJ, 667, L145 

Smith R., Fellhauer M., Goodwin S., Assmann P., 2011, MNRAS, 414, 
3036 



The fraction of star formation in bound clusters 35 



Solomon P. M., Rivolo A. R., Barrett J., Yahil A., 1987, ApJ, 319, 730 
Solomon P. M., Vanden Bout P. A., 2005, ARA&A, 43, 677 
Somerville R. S., Primack J. R., 1999, MNRAS, 310, 1087 
Somerville R. S., Primack J. R., Faber S. M., 2001, MNRAS, 320, 504 
Spitzer L., 1987, Dynamical evolution of globular clusters. Princeton, NJ, 

Princeton University Press, 1987, 191 p. 
Spitzer Jr. L., 1958, ApJ, 127, 17 
Springel V., Hernquist L., 2003, MNRAS, 339, 289 
Stanimirovic S., Staveley-Smith L., Jones P. A., 2004, ApJ, 604, 176 
Stil J. M., Israel F. P., 2002, A&A, 392, 473 
Strickland D. K., Stevens I. R., 1999, MNRAS, 306, 43 
Swaters R. A., Schoenmakers R. H. M., Sancisi R., van Albada T. S., 1999, 

MNRAS, 304, 330 
Tacconi L. J., Genzel R., Neri R., Cox P., Cooper M. C, Shapiro K., Bo- 
latto A., Bouche N., Bournaud F., Burkert A., Combes F, Comerford 
J., Davis M., Schreiber N. M. F, Garcia-Burillo S., Gracia-Carpio J., 
Lutz D., Naab T., Omont A., Shapley A., Sternberg A., Weiner B., 2010, 
Nature, 463, 781 

Tacconi L. J., Genzel R., Smail I., Neri R., Chapman S. C, Ivison R. J., 
Blain A., Cox P., Omont A., Bertoldi F., Greve T., Forster Schreiber 
N. M., Genel S., Lutz D., Swinbank A. M., Shapley A. E., Erb D. K., 
Cimatti A., Daddi E., Baker A. J., 2008, ApJ, 680, 246 

Tan J. C, Krumholz M. R., McKee C. F, 2006, ApJ, 641, L121 

Tasker E. J., Tan J. C, 2009, ApJ, 700, 358 

Tassis K., Christie D. A., Urban A., Pineda J. L., Mouschovias T. C, Yorke 

H. W., Mattel H., 2010, MNRAS, 408, 1089 
Testi L., Sargent A. I., Olmi L., Onello J. S., 2000, ApJ, 540, L53 
Thompson T. A., Quataert E., Murray N, 2005, ApJ, 630, 167 
Thornton K., Gaudlitz M., Janka H.-T, Steinmetz M., 1998, ApJ, 500, 95 
Tilley D. A., Pudritz R. E., 2007, MNRAS, 382, 73 
Toomre A., 1964, ApJ, 139, 1217 
Tutukov A. V., 1978, A&A, 70, 57 
van Dokkum P. G., Conroy C, 2010, Nature, 468, 940 
Vazquez-Semadeni E., 1994, ApJ, 423, 681 
Vesperini E., 2001, MNRAS, 322, 247 
Wada K., Norman C. A., 2007, ApJ, 660, 276 

Walsh W., Beck R., Thuma G., Weiss A., Wielebinski R., Dumke M., 

2002, A&A, 388, 7 
Wang K., Zhang Q., Wu Y, Li H.-b., Zhang H., 2012, ApJ, 745, L30 
Wang P., Li Z.-Y, Abel T., Nakamura F, 2010, ApJ, 709, 27 
Wang Q. D., Dong H., Lang C, 2006, MNRAS, 371, 38 
Weinberg M. D., 1994, AJ, 108, 1403 
White S. D. M., Frenk C. S., 1991, ApJ, 379, 52 
White S. D. M., Rees M. J., 1978, MNRAS, 183, 341 
Wilke K., Klaas U., Lemke D., Manila K., Stickel M., Haas M., 2004, 

A&A, 414, 69 
Williams J. P., McKee C. F, 1997, ApJ, 476, 166 

Wolfire M. G., McKee C. F, Hollenbach D., Tielens A. G. G. M., 2003, 
ApJ, 587, 278 

Wong T., Hughes A., Fukui Y, Kawamura A., Mizuno N., Ott J., Muller 
E., Pineda J. L., Welty D. E., Kim S., Mizuno Y, Murai M., Onishi T, 
2009, ApJ, 696, 370 

Wright E. L., 2006, PASP, 118, 1711 

Wiinsch R., Tenorio-Tagle G., Palous J., Silich S., 2008, ApJ, 683, 683 
Wylie-de Boer E., Freeman K., Williams M., 2010, AJ, 139, 636 
Yun M. S., Ho P. T. P., Lo K. Y, 1993, ApJ, 411, L17 
Zahn C. T., 1971, IEEE Transactions on Computers, 20, 68 

This paper has been typeset from a TpX/ KTpX file prepared by the 
author. 



